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A method  of  analyzing  the  modes  of  laser  resonators  con- 
taining homogeneous  media  is  presented  and  established.  This 
analysis  is  based  on  a derivation,  which  begins  with  Maxwell's 
equations  and  displays  the  required  assumptions,  of  a pair  of 
integral  equations  involving  the  tangential  fields  on  the  reson- 
ator mirrors.  This  pair  of  equations,  which  must  be  solved 
simultaneously,  is  specialized  to  apply  to  paraxial  resonators 
with  perfectly  conducting  mirrors.  The  result  of  the  special- 
ization is  a pair  of  integral,  eigenvalue  equations  for  the 
current  distributions  induced  on  the  resonator  mirrors. 

After  further  specializing  them  to  resonators  for  which 
the  spatial  dependence  of  the  modes  separates,  these  integral 
equations  are  solved  using  a straightforward  technique  based  on 
a variational  principle.  This  technique,  which  employs  a novel 
method  of  obtaining  modal  expansion  functions,  reduces  the  analysis 
to  a homogeneous  matrix  equation  that  is  solved  using  well-known 
numerical  methods. 

The  basic  theory  and  method  of  solution  presented  not  only 
produce  comparable  results  in  the  classical  cases  considered  by 
other  authors  in  this  field,  but  they  also  provide  a framework  for 
attacking  general  resonator  problems. 
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A THEORETICAL  ANALYSIS  OF  RESONATOR  MODES 
IN  THE  PRESENCE  OF  HOMOGENEOUS  MEDIA 

I . Introduction 

Since  the  discovery  of  the  laser  in  1961,  many  workers  have 
devoted  considerable  effort  to  analyzing  the  electromagnetic  fields 
associated  with  laser  resonators.  The  approaches  employed  in  these 
efforts  cover  the  spectrum  with  regard  to  derivation  of  basic  equations, 
essential  assumptions,  and  method  of  solution.  However,  despite  this 
wide  range  of  approaches  and  high  level  of  effort,  the  vast  majority  of 
these  approaches  has  not  been  derived  in  a completely  general  manner. 

In  addition,  all  of  these  approaches  leave  much  to  be  desired  with 
regard  to  ease  of  application,  depth  of  understanding  and  capability  of 
prediction  in  resonator  problems.  This  paper  is  directed  not  only 
towards  establishing  a sound,  general  theory  but  also  towards  improving 
the  existing  capabilities  in  these  three  aspects  of  the  analysis  of 
resonator  fields. 

The  improvement  in  the  ability  to  analyze  resonator  fields  is 
important  because  such  analyses  have  significant  impact  throughout  the 
development  and  use  of  laser  systems.  For  example,  although  the  first 
step  in  the  development  of  a new  laser  is  usually  the  finding  of  a gain 
medium,  the  next  step  is  devising  an  efficient  means  of  extracting  power 
from  the  medium  so  that  the  output  beam  has  certain  desirable  charac- 
teristics. These  desirable  characteristics  might  include  good  beam 
quality  (nearly  uniform  phase  and  amplitude  of  the  field),  low  losses, 
and  a large  volume  over  which  the  field  can  interact  with  the  active 
medium.  For  the  effective  design  of  resonators  with  these  or  other 


1 


I 

1 

extraction  characteristics,  one  needs  to  know  not  only  the  fields 
present  within  a given  resonator  but  also  how  they  are  affected  by 
changing  the  resonator  parameters. 

G.  Fox  and  T.  Li  (ref.  1,  pp.  453-488)  were  among  the  first  to 
recognize  the  importance  of  analyzing  the  fields  associated  with  lasers 
containing  homogeneous  media.  They  concentrated  their  efforts  on 
resonators  which  (with  regard  to  geometric  optics)  periodically  refocus 
paraxial  rays  so  that  they  always  remain  confined  to  the  resonator 
volume  (except  for  transmission  through  the  mirrors).  In  addition  to 
Fox  and  Li,  many  others  (refs.  2-10)  analyzed  the  fields  associated  with 
these  so-called  stable  resonators  by  calculating  their  normal  modes 
(which  could  be  used  to  represent  any  resonator  field).  These  normal 
modes  were,  by  definition,  the  eigensolutions  of  a certain  integral 
equation  that  was  derived  by  applying  the  theory  of  diffraction  to  the 
resonator  in  question. 

For  several  years,  it  was  felt  that  only  stable  resonators  could 
find  practical  application  in  the  laser  field.  However,  disadvantages 
such  as  small  mode  volume  and  poor  mode  discrimination  prompted  Siegman 
(ref.  11,  p.  278)  in  1965  to  propose  using  resonators  that  do  not 
confine  paraxial  rays  to  the  resonator  volume.  These  unstable  resona- 
tors, which  have  since  found  widespread  application,  are  characterized 
by  large  mode  volume,  diffraction  output  coupling,  good  transverse  mode 
discrimination,  and  totally  reflecting  optics.  This  last  characteristic 
is  especially  valuable  for  high  power  lasers  for  which  the  resonator 
mirrors  need  to  be  cooled.  | 


Just  as  laser  resonators  have  become  more  complex,  the  theory 
and  analytical  techniques  used  to  study  them  have  become  more  sophisti- 
cated. For  example,  the  initial  analyses  (refs.  11-16)  of  unstable 
resonators  were  based  on  geometric  optics.  However,  these  approaches, 
which  yield  average  values  of  the  losses  and  rough  estimates  of  the  mode 
distributions,  have  been  replaced  by  analyses  that  consider  the  diffrac- 
tive effects  introduced  by  the  finite  sizes  of  the  mirrors.  Although 
other  approaches  (refs.  17-19)  also  have  merit,  the  most  promising 
techniques  can  be  grouped  into  two  broad  categories;  the  Waveguide 
Analogy  (refs.  20-23)  and  the  Integral  Equation  Method  (refs.  6,  9,  10, 
24-28).  These  categories,  along  with  their  specific  deficiencies,  are 
described  in  chapter  II. 

Despite  the  many  positive  aspects  of  the  techniques  in  these  two 
categories,  no  method  exists  which  one  can  use  to  adequately  determine 
and  understand  the  characteristics  and  behavior  of  complex  resonator 
modes.  The  objective  of  this  work  is  to  develop  and  establish  such  a 
method  for  resonators  containing  a homogeneous  medium  and  two  perfectly 


conducting  mirrors.  The  purpose  of  this  report  is  to  present  this 
analytical  method  and  its  supporting  results. 

As  discussed  later  in  this  report,  this  analysis  provides  two 
significant  contributions  for  determining  and  understanding  laser 
resonator  modes.  The  first  is  a derivation  of  integral  equations  for 
the  tangential  fields  on  the  resonator  mirrors;  it  begins  with  Maxwell's 
equations  and  explicitly  displays  all  required  assumptions.  The  second 
contribution  is  a straightforward  technique  for  solving  these  equations 
for  a wide  variety  of  resonators  with  perfectly  conducting  mirrors. 

This  technique  includes  a novel  method  for  obtaining  expansion  functions 
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to  be  used  with  the  variational  principle  upon  which  the  technique  is 
based. 

This  report  is  organized  in  the  following  manner.  First, 
following  a presentation  of  important  background  material,  the  basic 
laser  resonator  problem  is  formulated.  The  formulation  is  followed  by 
derivations  of  the  integral  equations  for  the  current  induced  on  the 
mirrors  of  open  and  closed  laser  resonators.  Then,  in  chapter  IV,  the 
general  problem  of  solving  these  integral  equations  is  discussed,  and 
the  method  to  be  used  to  approximately  solve  these  equations  for  par- 
axial resonators  is  presented.  In  chapter  V,  the  theory  and  method  of 
soluticr:  of  the  previous  two  chapters  are  specialized  to  apply  to 
resonators  for  which  the  spatial  dependence  of  the  modes  can  be  sepa- 
rated to  yield  two  independent  governing  equations.  The  results  of  the 
specialization  are  presented  and  compared  to  existing  published  work  in 
the  following  chapter.  The  text  ends  with  chapter  VII,  which  presents 
specific  conclusions  and  recommendations.  The  eight  appendices  contain 
detailed  mathematical  derivations  and  calculations  in  support  of  this 
work. 
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II . Background  and  Theoretical  Preliminaries 

This  chapter  presents  and  then  formulates  the  basic  resonator 
problem.  To  that  end,  the  body  of  this  chapter  begins  with  a discussion 
of  important  background  material  and  a definition  of  a laser  resonator 
mode.  This  definition,  as  well  as  the  insight  gained  from  the  discus- 
sion which  follows,  forms  the  basis  for  the  formulation,  which  includes 
a discussion  of  the  analytical  approach  used  and  a presentation  of  the 
basic  equations  applied  in  the  following  chapters. 

Background 

Despite  the  large  number  of  configurations  in  use  today,  laser 
resonators  can  be  grouped  into  the  two  rather  broad  categories  of  open 
and  closed  resonators.  Since  the  analysis  presented  in  this  paper  is 
performed  by  first  considering  closed  resonators  and  then  extending  the 
results  to  open  resonators,  some  care  must  be  exercised  in  distinguish- 
ing between  the  two  resonator  types.  As  one  might  expect,  the  essential 
difference  between  the  two  types  involves  the  nature  of  the  resonator 
mirrors  and  the  surfaces  with  which  one  can  enclose  them. 

As  a first  step  in  determining  whether  a particular  resonator  is 
either  open  or  closed,  one  uses  one  or  more  closed  surfaces  to  enclose 
the  resonator  mirrors  such  that  only  the  mirrors  are  included  within  the 
enclosed  volume(s).  Thus,  for  each  case,  the  number  and  shape  of  the 
closed  surfaces  will  be  chosen  to  correspond  to  that  series  of  closed 
surfaces  which  most  nearly  conforms  to  the  mirror  shapes.  If  only  one 
closed  surface  is  used,  the  resonator  is  an  open  resonator  containing 
one  mirror  with  one  or  more  "holes"  in  the  mirror  surface.  If  two 
surfaces  are  used,  the  resonator  is  either  a closed  resonator  or  an  open 
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resonator  with  two  mirrors.  If  the  two  closed  surfaces  are  such  that 
one  surface  completely  encloses  the  second  surface,  the  resonator  is 
closed.  However,  if  one  surface  does  not  include  the  other,  the  reso- 
nator is  open.  Finally,  if  three  or  more  closed  surfaces  are  used  to 
enclose  the  resonator  mirrors,  the  resonator  is  an  open  resonator. 

To  clarify  these  ideas  concerning  open  and  closed  resonators, 
two  examples  are  shown  below.  First,  a planar  view  of  a closed  reso- 
nator is  shown  in  figure  1.  As  one  can  see  from  the  figure. 
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Figure  1.  A Planar  View  of  a Closed  Resonator 


the  surface  enclosing  the  unbroken  material  (mirror)  boundary,  which  may 
include  segments  of  several  materials,  consists  of  two  closed  surfaces 
and  $2*  As  indicated  in  the  previous  paragraph,  one  surface  (S^) 
completely  encloses  the  second  surface  ($2).  In  a similar  fashion,  a 
planar  view  of  an  open  resonator  containing  two  mirrors  is  shown  in 
figure  2.  From  this  figure,  one  can  see  that  neither  surface  encloses 
the  other. 
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Figure  2.  A Planar  View  of  an  Open  Resonator 

Open  Resonators  with  Spherical  Mirrors.  A significant  fraction 
of  the  open  resonators  in  use  today  employs  two  soherical  mirrors  of 
either  rectangular  or  circular  projection.  In  addition  to  the  material 
properties  of  the  host  and  lasing  media,  these  resonators  are  charac- 


terized by  the  following  mirror  parameters: 


1.  The  transverse  mirror  dimensions  (diameter  or 
width)  a and/or  b; 

2.  The  mirror  radii  of  curvature,  R;  and 

3.  The  mirror  separation,  L. 


These  last  two  parameters  play  an  important  role  in  the  cl  ass i- 


[ fi cation  of  laser  resonators  containing  spherical  mirrors.  For  that 


role,  the  mirror  separation  and  radii  of  curvature  are  combined  into  so- 
called  "o-oarameters"  as  shown  below. 


where  the  subscript  denotes  the  i^*^  mirror.  The  ranges  and/or  values  of 
these  g-parameters  for  several  resonators  are  shown  below  in  table  I. 


Table  I 

Resonator  g-parameters 


RESONATOR  TYPE 

g-PARAMETERS 

STABLE 

0 < q^g2  < 1 

UNSTABLE 

g^g2  < 0 OR  a^g2  > 1 

QUASI -STABLE 

g^g2  = 0 OR  9]P2  " ^ 

CONFOCAL* 

2g^g2  " 9] 

PLANE  PARALLEL 

' 92  = 1 

*A  confocal  resonator  is  a resonator  for  which  the 
foci  of  the  two  resonator  mirrors  are  colocated. 


The  resonator  Fresnel  number  N is  another  parameter  that  is 
often  used  to  characterize  these  spherical  resonators.  Actually,  a 
Fresnel  number,  which  is  a quantity  often  reserved  for  mirrors  of 
circular  projection,  is  designated  for  each  resonator  mirror  and  is 
defined  by 

where  a is  the  mirror  radius  and  x is  the  wavelength  of  the  radiation 
within  the  resonator.  To  adapt  this  quantity  to  mirrors  of  rectangular 
projection,  one  uses  Eq.  (2)  to  define  two  Fresnel  numbers  for  each 
mirror.  However,  for  this  case,  the  quantity  "a"  denotes  each  of  the 
transverse  mirror  dimensions  (length  and  width)  instead  of  the  mirror 
radius. 
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Even  thouph  spherical  mirror  resonator  configurations  of  the 
types  shown  in  table  I are  widely  used,  they  do  have  several  short- 
comings. Unfortunately,  these  shortcomings  are  not  limited  to  the  small 
mode  volume  and  poor  mode  discrimination  exhibited  by  stable  resonators. 
For  instance,  for  many  of  the  applications  involving  unstable  resonators 
with  mirrors  of  circular  projection,  the  output  beam  has  an  annular 
shape.  This  annular  shape  not  only  makes  the  beam  difficult  to  use,  but 
it  can  also  prevent  a nontrivial  (several  percent  or  more)  portion  of 
the  energy  from  reaching  the  far-field  central  spot.  Also,  spherical 
unstable  resonators  do  not  efficiently  accommodate  new  lasers  having 
central  obscurations  (to  the  radiation  within  the  resonator)  caused  by 
the  use  of  a radial  gas  flow  or  a radial  electron  beam.  These  and  other 
deficiencies  have  led  to  the  study  and  development  of  some  new  and 
rather  exotic  resonator  configurations. 

Toroidal  Resonators.  Many  of  these  new  resonator  configurations 
involve  resonators  which  include  at  least  one  toroidal  mirror.  As  the 
name  implies,  the  mirrors  falling  into  this  category  simply  correspond 
to  different  portions  or  cross  sections  of  a toroid;  however,  there  are 
two  kinds  of  toroidal  mirrors.  The  first  kind,  which  is  rarely  used,  is 
simply  characterized  by  two  radii  of  curvature.  The  second  kind  is  a 
mirror  that  is  characterized  by  a particular  surface  of  revolution;  that 
is,  a plane  which  passes  through  the  resonator  axis  and  intersects  the 
mirror  yields  two  arcs.  These  two  arcs,  which  have  radii  of  curvature  R 
and  which  may  or  may  not  be  connected,  have  centers  of  curvature  which 
are  displaced  from  the  resonator  axis  by  a distance  pg. 
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As  with  spherical  mirror  resonators,  one  may  use  a pair  of  q- 
parameters  to  characterize  or  classify  toroidal  resonators.  However, 
for  these  resonators,  one  modifies  Eq.  (1)  by  replacing  the  radius  of 
curvature  of  a spherical  mirror  with  the  radius  of  curvature  of  the 
individual  arcs  discussed  above.  With  this  modification,  one  can  use 
table  I to  classify  toroidal  resonators  with  regard  to  stability,  etc. 
Although  this  method  of  classifying  toroidal  resonators  is  somewhat 
artificial,  it  does  provide  a useful  framework  for  the  analysis. 

Four  resonators  using  this  second  kind  of  toroidal  mirror  are 
depicted  in  figure  3.  Actually,  each  of  the  sketches  in  this  figure 
shows  the  intersection  of  a toroidal  resonator  with  a plane  passing 
through  the  resonator  axis  (denoted  AA'  in  the  figure).  To  obtain  a 
three  dimensional  view  of  each  resonator,  it  is  necessary  to  revolve 
each  sketch  through  an  angle  of  180°  about  the  resonator  axis.  In 
figure  3(d),  the  lines  with  the  tic-marks  designate  the  boundaries  of  a 
central  obscuration  to  the  radiation  within  the  resonator.  The  lines 
and/or  arcs  in  each  of  the  other  three  sketches  depict  a wavefront  as  it 
passes  through  ard  out  of  that  resonator. 

Past  Work.  Existina  analytical  techniques  are  deficient  in  that 
they  cannot  adequately  treat  many  spherical  or  toroidal  resonator  con- 
figurations. This  is  especially  true  of  the  Waveguide  Analogy 
(refs.  20-23)  in  which  the  resonator  is  treated  as  if  it  were  a wave- 
guide section  (witn  the  resonator  mirrors  as  the  guiding  surface)  which 
is  coupled  to  the  surrounding  space.  The  auide  is  taken  to  be  operat- 
ing near  cut-off,  and  the  modes  are  expanded  in  terms  of  the  fields 
which  would  be  present  if  the  waveguide  section  were  infinite.  The 
expansion  coefficients  are  calculated  by  first  determining  how  the 
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Fiqure  3.  Toroidal  Resonators 
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infinite  guide  modes  are  reflected  by  the  interfaces  present  between  the 
guide  section  and  its  surroundings,  and  then  requiring  that  the  resulting 
fields  be  self-reproducing.  Of  the  two  approaches  mentioned  earlier,  it 
is  this  Waveguide  Analogy  which  lends  the  greatest  insight  into  the 
physical  behavior  of  the  modes.  However,  since  its  application  not  only 
requires  solving  rather  involved  transcendental  equations  in  the  complex 
plane,  but  also  calculating  reflection  coefficients  at  all  interfaces 
(between  the  guide  and  its  surroundings),  it  would  be  very  difficult  to 
apply  to  resonators  with  output  coupling  apertures  or  toroidal  mirrors. 

Despite  these  limitations,  the  Waveguide  Analogy  has  been  used 
to  excellent  advantage  by  Vainshtein  (ref.  21,  p.  711)  to  obtain  closed 
form  expressions  for  the  modes  of  finite,  plane  parallel  resonators.  As 
will  be  seen  on  page  63,  these  expressions  play  an  important  role  in 
the  techniques  applied  in  this  paper. 

The  second  and  more  widely  used  analytical  approach  is  based  on 
solving  certain  integral  equations  for  components  of  the  electromagnetic 
field  (precisely  which  components  depends  on  the  problem)  on  the  resona- 
tor mirrors.  The  popularity  of  this  approach  is  largely  a result  of  the 
availability  of  mathematical  methods  for  approximately  solving  integral 
equations;  it  is  also  due  to  the  fact  that  once  the  basic  forms  have 
been  derived,  the  equations  can  be  easily  modified  to  apply  to  a wide 
variety  of  complex  resonators. 

Several  different  authors  (refs.  1,  29-31)  have  derived  integral 
equations  which  are  widely  considered  to  be  the  basic  forms  used  in  this 
approach.  However,  each  of  these  derivations  and,  in  some  cases,  tne 
equations  themselves  are  deficient  in  one  or  more  aspects  of  the  prob- 
lem. One  major  deficiency  of  many  of  the  derivations  is  the  tendency  to 
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make  important  assumptions  at  the  beginning  of  the  analysis  concerning 
the  nature  of  the  electromagnetic  field.  This  tendency  not  only  limits 
the  applicability  of  the  entire  analysis,  but  it  can  inhibit  under- 
standing the  behavior  of  the  field  by  obscuring  the  precise  implications 
of  the  various  assumptions.  For  example,  many  derivations  are  based  on 
the  implicit  assumption  that  the  mirror  material  properties  do  not 
affect  the  mode  distributions. 

This  assumption,  which  may  not  always  be  valid,  obscures  the 
fact  that  in  determining  the  mode  distributions  one  must  consider  the 
effect  of  both  the  electric  and  magnetic  fields  on  the  mirror  surfaces. 
For  some  cases,  either  the  electric  or  magnetic  field  at  the  mirror 
boundary  may  actually  have  a negligible  effect  on  the  mode  distribu- 
tions. However,  one  should  show  that  this  is  true  in  each  case  rather 
than  assuming  it  is  true  in  general.  A second  example  of  these  ob- 
scuring assumptions  is  the  assumption  that  the  field  (or  current)  on  one 
resonator  mirror  can  be  expressed  entirely  in  terms  of  the  field  (or 
current)  on  the  second  resonator  mirror.  This  assumption  inherently 
obscures  the  fact  that,  in  some  cases,  one  must  include  the  effect  of 
the  current  on  both  resonator  mirrors. 

A second  deficiency  applies  to  derivations  that  employ  assump- 
tions which  are  inconsistent  with  known  mathematical  theorems.  One  such 
assumption  is  that  a scalar  field,  along  with  its  normal  derivative, 
vanishes  identically  on  a finite  surface  element.  Finally,  many  deri- 
vations are  deficient  in  that  they  do  not  properly  include  the  damping 
of  the  modes  of  open  resonators  and  its  effect  on  the  resulting  mode 
distributions.  For  some  cases,  the  damping  has  a negligible  effect  on 
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the  actual  distributions;  however,  in  other  cases,  its  effect  can  be 
significant.  By  including  the  effects  of  damping  at  the  beginning  of 
the  analysis,  one  can  show  that  it  is  not  always  possible  to  formulate 
the  laser  resonator  problem  in  terms  of  a linear  eigenvalue  equation. 

Unfortunately,  the  bulk  of  the  numerical  work  performed,  especially  for 
complex  resonators,  has  been  based  on  derivations  containing  this  last  | 

deficiency. 

As  indicated  in  the  introduction,  the  remainder  of  this  paper  is 
directed  toward  developing  and  establishing  an  analytical  approach  which 

,] 

will  overcome  many  of  the  deficiencies  of  the  existing  theories. 

Laser  Resonator  Modes 

In  a wide  range  of  problems  in  electromagnetic  theory,  it  is 
standard  procedure  to  work  in  terms  of  fields  (modes)  which  depend  only 
on  the  characteristics  of  the  material  bodies  and  surrounding  media. 

This  useful  procedure  was  adapted  to  this  analysis  by  defining  a reso- 
nator mode  to  be  a member  of  that  class  of  linearly  independent,  source 
free,  electromagnetic  fields  which  satisfy  the  boundary  conditions 
imposed  by  the  resonator. 

To  gain  some  insight  into  the  nature  of  these  modes,  consider 
the  following  experiment.  A radiation  source  is  placed  in  the  vicinity 
of  a resonator.  Prior  to  time  t=0,  the  source  is  turned  off,  and  only 
the  null  field  is  present.  Then  at  t=0,  the  source  is  pulsed  and 
immediately  turned  off.  After  sufficient  time  has  passed  for  the  wavefront 
to  reach  and  interact  with  the  resonator,  the  resulting  field,  which  is 
composed  entirely  of  the  modes  of  the  resonator,  is  expressed  in  terms 
of  the  Fourier  decomposition  shown  in  Eq.  (3). 
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To  analyze  the  temporal  behavior  of  that  field,  the  function  ^(F,u)  is 
analytically  continued  into  the  complex  oj-plane,  and  the  resulting  form 
is  examined  using  contour  integration  in  conjunction  with  the  theory  of 
residues.  That  examination  reveals  that  since  the  field  is  nonvanishing 
after  time  t=0,  the  integrand  must  have  singularities  in  the  upper  half 
of  the  complex  a}-plane.  The  nature  of  these  singularities,  which  occur 
for  u).  = 0)1  + io)'.'  (a)'.'>0),  can  be  used  to  determine  the  behavior  of 

J J J J 

the  field  through  the  relation 

^(r,t)  = Ziri  23  Lj(r,u),t)  (4) 

j 

where  ri(F,o),t)  is  the  residue  of  the  integrand  at  o)..  Although  the 
J w 

exact  form  depends  on  the  details  of  the  source  and  the  resonator,  each 
residue  can  be  written  in  the  form 


+io)lt  -a)'.'t 

ri(r,w,t)  = e e 


(i>) 


where  F-Cr,*!).)  is  the  residue  of  ^(r,o))  at  oj  = o)-.  The  result  obtained 

J J J 

by  substituting  Eq.  (5)  into  Eq.  (4)  expresses  the  field  as  a combina- 
tion of  the  modes  of  the  resonator,  where  each  mode  has  a time  depen- 
+io)lt  -lol't 

dence  of  the  form  e 3 e J . In  this  last  expression,  o)l  is  the 

J 

frequency  and  a)'l  is  the  decay  constant  of  the  j^*^  mode. 

J 

Another  important  aspect  of  the  nature  of  these  modes  involves 
the  type  of  energy  flow  associated  with  each  mode.  For  this  discussion, 
the  two  types  of  energy  flow  of  importance  are  (1)  energy  flow  from  the 
resonator  to  the  surroundings  and  (2)  energy  flow  from  the  surroundings 
to  the  resonator.  Of  course,  in  the  typical  laser  application. 


lb 


radiation  emanates  from  within  the  resonator,  and  the  net  radiative 
energy  flows  from  the  resonator  to  the  surroundings. 

For  closed  resonators,  such  as  the  one  shown  below  in  figure  4, 
distinguishing  between  the  two  types  of  energy  flow  is  relatively 
straightforward.  For  these  resonators,  which  require  two  closed  sur- 
faces to  enclose  the  unbroken  material  (mirror)  boundary,  one  distin- 
guishes between  the  two  types  of  flow  by  considering  energy  flow  across 
the  unbroken  material  boundary. 


Figure  4.  A Closed  Resonator 


However,  for  resonators  not  characterized  by  an  unbroken  material 
boundary,  these  ideas  need  to  be  made  more  precise.  To  do  that,  con- 
sider the  two  closed  surfaces  involving  the  open  resonator  shown  in 
figure  5.  The  first  surface,  which  I call  the  resonator  hull,  is 
designated  H in  the  figure.  This  surface  consists  of  the  front  surfaces 
of  the  two  mirrors  and  the  family  of  straight  lines  connecting  the 
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mirrors  such  that  H encloses  the  maximum  possible  volume.  The  second 
surface,  which  I call  the  resonator  shell,  is  designated  S in  the 
figure.  The  resonator  shell  consists  of  the  front  surfaces  of  the  two 
mirrors  and  the  family  of  straight  lines  enclosing  the  maximum  volume 
such  that  each  straight  line  begins  and  ends  at  a mirror  edge. 


Figure  5.  Surfaces  Defining  the  Inside  and 
Outside  of  an  Open  Resonator 


If  the  resonator  hull  encloses  a volume  greater  than  that 
enclosed  by  the  resonator  shell  (which  might  correspond  to  the  presence 
of  a lip  on  one  of  the  mirrors),  the  ideas  related  to  energy  flow  are 
not  easily  understood,  and  such  a resonator  will  not  be  treated  in  this 
analysis.  However,  in  the  event  that  the  resonator  hull  and  shell 
coincide,  the  important  energy  flow  ideas  can  be  precisely  formulated. 

To  do  this,  the  volume  contained  within  and  on  the  closed  surface  H, 
designated  region  a,  is  considered  to  be  inside  the  resonator.  Obviously 
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then,  the  volume  outside  H,  designated  region  b,  is  outside  the  reso- 
nator. For  this  analysis,  only  those  fields  for  which  the  net  energy 
flows  from  region  a to  region  b will  be  designated  as  laser  resonator 
modes. 


Formulation  of  the  Problem 

General  Approach.  Based  on  the  material  covered  in  the  previous 
section,  the  modes  of  a laser  resonator  can  be  found  by  determining  the 
linearly  independent  members  of  that  class  of  electromagnetic  fields, 
with  time  dependence  e^^“  ^ e”‘^  which  satisfy  the  boundary  conditions 
imposed  by  the  resonator  to  produce  energy  flowing  from  the  resonator  to 
the  surroundings.  Figure  6 depicts  a planar  view  of  this  basic  problem, 
which  is  analyzed  in  the  remainder  of  this  paper. 


Figure  6.  The  Laser  Resonator  Problem 
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As  indicated  in  the  figure,  a rectangular  system  of  coordinates 
(only  the  x-z  plane  is  shown  in  the  figure)  is  established  in  the  vicinity 
of  the  mirrors  with  the  z-axis  chosen  as  the  optic  axis  of  the  resonator. 
In  addition  to  the  two  mirrors,  the  resonator  contains  a homogeneous 
medium  with  parameters  u,  e,  and  a.  The  position  vector  of  a point  on 
the  mirror  is  denoted  r. , and  the  distance  between  two  arbitrary 

J 

points  on  the  two  mirrors  is  R.j2.  Finally,  the  mirror  separation, 
measured  along  the  optic  axis,  is  denoted  L.  Figure  6 a-lso  applies  to 
closed  resonators,  except  in  that  case,  there  is  one  mirror  rather  than 
two  separate  ones. 

The  following  approach  is  used  to  determine  the  fields  and  other 
characteristics  associated  with  the  laser  resonator  depicted  in  figure  6. 
First,  the  basic  equations  of  electromagnetic  theory  are  used  to  derive 
a pair  of  coupled  integral  equations  for  the  electric  and  magnetic 
fields  within  the  resonator.  These  equations  express  these  fields 
within  the  resonator  volume  in  terms  of  the  electric  and  magnetic  fields 
tangential  to  the  mirror  surfaces.  This  pair  of  equations  is  then 
specialized  to  apply  to  resonators  with  perfectly  conducting  mirrors. 

The  result  of  this  specialization  is  a pair  of  equations  which  relate 
the  electric  and  magnetic  fields  within  the  resonator  to  the  currents 
induced  on  the  perfectly  conducting  mirrors.  Then,  using  the  boundary 
conditions  for  perfect  conductors  and  letting  the  field  point  approach  a 
point  on  the  mirror  surfaces,  an  integral  equation  is  obtained  for  the 
current  distributions  on  the  resonator  mirrors.  Once  this  equation  has 
been  solved,  these  current  distributions  can  be  used  in  the  equations 
that  relate  the  electric  and  magnetic  fields  to  the  currents  induced  on 
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the  resonator  mirrors.  In  this  manner,  one  may  specify  the  electro- 
magnetic fields  associated  with  each  mode  throughout  the  resonator 
volume.  The  fields  obtained  using  this  approach  are  unique  (Harrington, 
ref.  32,  p.  102). 

Basic  Equations.  Since  the  modes  are  electromagnetic  in  nature, 
they  satisfy  Maxwell's  equations,  which  are  written  below  in  mks  units 
for  uniform  media  (Stratton,  ref.  33,  p.  464). 


- V X F(r,t)  = + uR’(r,t)  + fT(r,t) 

(6) 

• 

V X iT(7,t)  = o’(F,t)  + T(F,t) 

(7) 

V . lT(F,t)  = \/y 

(8) 

7 • 0(F,t)  = p 

(9) 

where  C = EqF  + V.  In  these  equations,  V is  the  polarization  of  the 
medium,  M is  a fictitious  magnetic  current,  and  q^^  is  a fictitious 
magnetic  charge.  These  equations  can  be  manipulated  to  yield  the 
following  equations  of  continuity. 

V . J(r,t)  = -p(r,t)  (10) 

V . lT(r,t)  = -\(r,t)  (11) 

In  addition  to  Eqs.  (6)  through  (11),  the  following  two 
equations  govern  the  behavior  of  the  fields  in  two  different  regions 
separated  by  a surface  electric  or  magnetic  current  (Harrington, 
ref.  32,  p.  34). 
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n X [H^  - H2]  = (12) 

[E^  - £3]  X n = (13) 

where  n,  the  unit  normal  to  the  surface,  points  into  region  #1. 

To  apply  these  equations  to  the  resonator  problem,  the  following 
conditions  concerning  the  medium  and  fields  within  the  resonator  are 
assumed  to  hold: 

1.  The  homogeneous  medium  is  assumed  to  be  linear  and  isotropic 

so  that  D"  = (n^  + where  n is  the  index  of  refraction 

of  the  host  material,  and  x is  the  susceptibility  of  the 
lasing  material . 

2.  Within  the  medium,  no  magnetic  sources  or  electric  charge 

exists.  Then  p = = 0 and  M = 0. 

3.  Within  the  medium,  J = aF. 

4.  For  each  mode,  the  fields  have  a time  dependence  of  the  form 

+icjit 

e where  w.  = w.  + iu-";  oj.'S  0. 

J J J J 

With  these  assumptions,  the  field  equations  within  the  resona- 
tor become 


V X Fj(r)  = - iojjvHj(f) 

(14) 

V X Fj(r)  = ia)je(a)j)F.(r) 

(15) 

V • H.(r)  =0 

\j 

(16) 

V . E.(r)  = 0 

(17) 
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where  e(o)j 


j)  = Bq 


10 

“j"0j 


. By  taking  the  curl  of  Eq.  (14), 


substituting  Eq.  (15)  into  Eq.  (14),  and  simplifying  the  result,  the 
wave  equation  for  T.  can  be  obtained  as 

J 


v2rj.(?)  + kj2  Fj(F)  = 0 


(18) 


where  = a)j2ye((j^).  a similar  procedure  yields  an  equation  identical 
to  Eq.  (18)  with  E".  replaced  by  H" . . 

J J 

These  wave  equations  form  the  basis  for  the  derivation  of  the 
modal  integral  equations  in  the  following  chapter. 
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TII . Derivation  of  the  Integral  Equations 


The  purpose  of  this  chapter  is  to  derive  the  basic  integral 
equations  which  govern  the  behavior  of  the  modal  fields  of  laser  reso- 
nators. The  text  begins  with  derivations  of  these  integral  equations 
for  two  classes  of  closed  resonators.  The  derivations  are  then  modified 
to  apply  to  open  resonators.  Following  these  modifications,  the  integral 
equations  for  open  resonators  are  specialized  to  apply  to  resonators 
with  two  perfectly  conducting  mirrors.  The  specialization  is  then 
manipulated  to  yield  the  integral  equations  tc  be  solved  for  the  modal 
currents  induced  on  the  resonator  mirrors. 


Integral  Equations  for  Closed  Resonators 


The  two  classes  of  closed  resonators  to  be  considered  in  this 
section  are  closed  resonators  with  either  homogeneous  or  segmented 
boundaries.  As  one  might  expect,  a homogeneous  boundary  consists  of  a 
single  homogeneous  medium.  On  the  other  hand,  a segmented  boundary 
contains  at  least  two  materials  with  different  electromagnetic  prop- 


erties. 


For  closed  resonators  containing  either  homogeneous  or  segmented 


boundaries,  it  is  assumed  that  the  boundary  is  smooth.  This  assumption, 
which  also  applies  to  open  resonators,  implies  that  (Taylor,  ref.  30, 
pp.  360,  371) 

1.  The  boundary  (surface)  does  not  intersect  itself, 

2.  The  boundary  (surface)  has  a tangent  plane  at  each  point 
whose  direction  varies  continuously  as  the  point  moves  along 


the  boundary. 
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It  is  further  assumed  that  any  discontinuities  in  the  fields  are  such 
that,  at  the  points  of  discontinuity,  the  changes  in  the  fields  are 
finite.  These  field  discontinuities  could  occur  at  the  intersection 
of  the  different  materials  in  a closed  resonator  with  a segmented 
boundary  or  at  the  edge  of  a mirror  in  an  open  resonator.  It  is  also 
assumed  that  the  resonator  fields  are  differentiable  at  all  points  in 
space  with  the  exception  of  these  points  of  discontinuity  along  the 
resonator  boundary. 

Closed  Resonators  with  Homogeneous  Boundaries.  As  indicated  at 
the  end  of  the  previous  chapter,  the  wave  equations  for  E".  and  H.  form 
the  basis  for  the  integral  equation  derivations  contained  in  this 
chapter.  Those  two  vave  equations  are  summarized  by  a single  equation 

v'2  KAP)  + k.2  S'.(r')  = 0 (19) 

J sj  U 

where  A-  = E.  or  H”.  as  desired  and  k.^  = a).2ye(co.).  Equation  (19)  can 

J 0 J yj  yj 

be  solved  using  Green's  function  techniques  by  first  solving  the 
equation 

7'2p,  + k.2  G = - ! 6(7:7^  (20) 

J 

where  5(r-r' ) is  the  Dirac  delta  function  (Collin,  ref.  34,  p.  565)  and 
^ ^y®y  ®z^z  dyadic  in  rectangular  coordinates. 

One  solution  to  this  equation  is  a Green's  function  corresponding  to 
outgoing  waves  (see  appendix  A)  of  the  form  G = T(}>(R)  where 
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*(R)  = i 


-ik.R 
e ^ 


V(n^  + xj)-^ 


with  R = |r-r‘ 


To  apply  this  Green's  function  to  a closed  resonator  with  a 
homogeneous  boundary,  the  scalar  products  of  S with  Eq.  (19)  and  A.(r) 

J 

with  Eq.  (20)  are  taken,  the  two  resulting  equations  are  subtracted,  and 
the  difference  is  integrated  over  the  volume  inside  the  resonator  to 
obtain  Eq.  (23). 


A.(r)  = I V'2A.(r')  • P,  - A.(r')  • v'2  G dt'  (23) 
J I J J 


Applying  the  identity  (Collin,  ref.  34,  o.  60), 


V • X (v'x  G)  + (v‘  X S"j)  X G + Sj.  (V  • G)  - (V  • ^j) 


= A.  • v'2G  - v'2T.  • G 
J J 


in  conjunction  with  the  divergence  theorem  applied  to  dyadics  (Collin, 
ref.  34,  p.  569),  the  volume  integral  is  converted  to  the  surface  inte- 
gral shown  below 


^j(r)  (n  X Kj) 


7'xG+(nxv'xA.)  *G 

J 


+ (n  • . G - (n  . G)v  - S'.  dS' 

«J  J 
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where  n points  into  the  resonator  volume,  and  S'  is  the  surface  just 
inside  the  resonator  boundary.  Then  substituting  iTj  and  for  and 
using  Eqs.  (14)  and  (15)  to  eliminate  v x F.  and  7 x H. , the  following 

J J 

pair  of  equations  is  obtained. 


Ej(f) 


i 


(n  X F.)  • V'  X S - iu^ii  (n  x FT.)  • § 

J J s/ 


+ (n  • F. )v'  • 6 

J 


Tj(?)  =/  (S  X Hj 
S' 


dS' 


(26) 


) • 7'  X 5 + ioj-e  (h  X F. ) • S 

J J 


+ (n  • H.)7'  • 5 

V 


dS' 


(27) 


Equations  (26)  and  (27)  are  the  integral  equations  which  govern 
the  behavior  of  the  modal  fields  within  closed  resonators  with  homoge- 
neous boundaries.  However,  before  applying  these  equations  throughout 
the  resonator  volume,  they  must  be  specialized  to  apply  to  points  on  the 
surface  S',  and  the  specializations  must  be  solved  simultaneously  for 

(n  X F,. ) and  (n  x H. ).  When  the  surface  fields  have  been  determined, 

J J 

Eqs.  (26)  and  (27)  can  be  used  to  compute  the  fields  within  the  resonator 
volume. 


To  compute  the  fields  outside  the  resonator,  one  must  specify 
how  the  fields  behave  as  the  surface  S'  is  crossed.  This  specification, 
which  simply  amounts  to  describing  the  particular  resonator  boundary  or 
interface,  will  yield  known  values  for  the  field  over  the  closed  surface 
just  outside  the  resonator  boundary.  This  knowledae  of  the  tangential 
fields  over  the  closed  surface  uniquely  specifies  the  fields  throughout 
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any  bounded  volume  outside  the  resonator  (Harrington,  ref.  32, 
p.  102). 

It  is  interesting  to  note  that  Eqs.  (26)  and  (27)  can  be  applied 
directly  to  a closed  resonator  with  a perfectly  conducting  boundary 
(n  X F = 0).  No  modification  of  the  equations  at  any  point  is  required 
as  the  functions  described  by  the  closed  surface  integrals  change 
discontinuously  as  the  surface  S'  is  crossed  (from  region  a to  region  b 
in  figure  5).  These  discontinuous  changes  are  such  that  the  fields 
outside  S'  vanish  identically  (Stratton,  ref.  33,  p.  468),  thus  corre- 
sponding to  the  fields  of  a closed  conducting  boundary. 

Closed  Resonators  with  Segmented  Boundaries.  In  the  above 
derivation,  the  divergence  theorem  was  applied  to  a closed  surface 
separating  the  resonator  from  its  surroundings.  That  application  is 
completely  justified  as  the  fields  and  their  derivatives  are  continuous 
throughout  the  volume  of  a closed,  homogeneous  resonator  with  a smooth, 
homogeneous  boundary.  However,  if  the  resonator  boundary  is  segmented, 
discontinuities  in  the  fields  and  their  derivatives  may  occur  at  the 
intersection  of  the  different  materials  in  the  boundary  surface.  As  a 
result,  any  derivation  applied  to  a resonator  with  a segmented  boundary 
must  allow  for  the  presence  of  discontinuities  in  the  fields  at  points 
along  the  resonator  boundary. 

As  discussed  in  Stratton  (ref.  33,  p.  468),  the  presence  of 
such  discontinuities  can  be  reconciled  with  the  field  equations  only  if 
one  assumes  the  existence  of  charges  or  currents  at  the  points  of  dis- 
continuity. These  sources  produce  fields  which,  when  added  to  the  fields 
in  Eqs.  (26)  and  (27),  yield  net  fields  that  satisfy  Maxwell's  equations 
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(Baker  and  Copson,  ref.  35,  pp.  114-117;  appendix  B of  this  paper). 
These  net  fields  are  given  by  Eqs.  (28)  and  (29)  below. 


m J 


(n  X Ejj^)  • V X R + (n  X 7'  X Ej^^)  • 5 


+ (n  . EjjV  . R dS; 


-7^  L /’■ 
110. e „ J 


jm  m 


J m 


Hj(r)  = 2 y*  (n  X Hjj^)  • V x (5  + (n  x 7'  x H^^^)  • ^ 


+ (n  • F.  ) 7'  • 5 I dS* 
jm  1 m 


Ep,  • 

m m 


where  S'  is  the  surface  of  the  m continuous  seament  of  the  boundary, 
m 

is  the  contour  enclosina  S',  and  and  fL  are  the  fields  just 
m m m m 

inside  C . 
m 

Thus,  for  resonators  with  segmented  boundaries,  it  is  Eqs.  (28) 
and  (29)  which  must  be  specialized  to  apply  to  the  tangential  fields  on 
the  surface  S'.  Once  these  specializations  are  obtained,  the  procedure 
for  computing  the  fields  at  other  points  in  space  is  the  same  as  the  one 
described  above  for  closed  resonators  with  homogeneous  boundaries. 
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However,  for  either  case,  the  final  equations  may  be  simplified 
somewhat.  This  simplification  results  from  takinq  the  curl  of  the 
expressions  for  Tj  and  Rj,  applying  the  fact  that  v'  • 5 = 7'<(»  = - 7<(>, 
and  using  Eqs.  (14)  and  (15).  The  final  result  is  the  pair  of  equations 
shown  below. 


ffj(F). 


V xi5x^+7x5 


(n  X Hj) 


dS' 


(30) 


■ / pT  (S  x 

S'  I'-j' 


) • 7'  X X V + 7 X ^ • (n  X F.)}  dS'  (31) 


where  the  surface  integrals  are  evaluated  over  the  closed  boundary  with 
any  points  of  discontinuity  removed.  (See  appendix  B.) 


Integral  Equations  for  Open  Resonators 

For  open  resonators,  these  integral  equations  can  be  simplified 
further  by  eliminating  that  portion  of  the  closed  surface  integral  which 
corresponds  to  the  open  surface  of  the  resonator  hull.  To  do  that,  the 
open  resonator  is  considered  to  be  a special  case  of  a closed  resonator 
with  a segmented  boundary.  Then,  Eqs.  (28)  and  (29)  are  applied  to 

1.  The  resonator  hull,  and 

2.  The  closed  surface  consisting  of  the  backs  of  the  mirrors 
and  the  open  portion  of  the  hull. 

The  equations  for  F and  FT  in  these  two  cases  are  then  subtracted  to  yield 
the  following  pair  of  equations  for  the  fields. 
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The  next  step  is  to  specialize  these  equations  to  apply  to  laser 
resonator  modes  at  points  inside  the  resonator  hull.  For  these  fields, 
the  effect  of  the  fields  on  the  backs  of  the  mirrors  (S^'  and  S^')  can 
be  dropped.  This  is  justified  because  S^'  and  S^'  are  outside  the  reso- 
nator, and  laser  resonator  modes  are  fields  for  which  the  vast  majority 
of  the  energy  flows  from  the  resonator  to  the  surroundings. 

It  should  be  apparent  that  this  procedure  (of  neglecting  S^'  and 
S^')  will  yield  accurate  results  only  to  the  degree  that  laser  resonators 
actually  produce  fields  such  that  the  energy  flow  from  the  surroundings 
is  negligible.  For  some  or  perhaps  all  resonators,  there  may  be  some 
modes  for  which  n2  energy  flows  into  the  resonator.  However,  for  a 
resonator  operating  at  a low  frequency,  it  is  easy  to  imagine  that  the 
fields  on  the  backs  of  the  mirrors  could  affect  the  field  inside  the  hull. 
Conversely,  for  resonators  containing  mirrors  much  larger  than  the  wave- 
length of  operation  (which  corresponds  to  the  conditions  for  the  physical 
optics  approximation),  the  effect  of  the  fields  on  S^'  and  S^'  should  be 
negligible.  As  a result,  they  will  be  neglected  for  the  remainder  of 
this  analysis,  and  the  sums  in  Eqs.  (32)  and  (33)  will  be  evaluated 
from  m=l  to  m=2. 

As  with  the  closed  resonator  case,  one  takes  the  curl  of  the  two 
equations,  applies  7 x Vif  = 0,  and  uses  Eqs.  (14)  and  (15)  to  eliminate 
V X and  v x ITj.  The  resulting  equations,  which  are  the  chief  result 
of  this  section,  are  shown  below. 


Hj(f) 


Ts  f\^  X X ^ + 7 X ^-(nxHj^) 

m=l  I ^ 

m 


(34) 


Ll4j 
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(35) 


F,(F) 


^ (nxH.  )-V  X X V + V X ^i-(nxE.  ) dS' 
(DjE  jm'  ' jm'  ( m 


Open  Resonators  with  Perfectly  Conducting  Mirrors 

For  resonators  with  perfectly  conducting  mirrors,  the  fields 
on  the  mirror  surfaces  satisfy  the  boundary  conditions 

n X E.(r')  = 0 (36) 

J 


n X Hj(r')  = 


(37) 


Substituting  these  two  equations  into  Eqs.  (34)  and  (35),  the  basic 
equations  for  this  case  take  the  form  shown  below. 


2 


m 


Now  if  Eqs.  (37)  and  (38)  can  be  manipulated  to  yield  forms  which  can  be 
solved  for  the  then  Eqs.  (38)  and  (39)  can  be  used  to  compute  the 
fields  throughout  the  resonator  volume.  The  procedure  for  obtaining 
equations  for  the  is  the  subject  of  the  next  several  pages,  while 
the  methods  used  for  solving  these  equations  are  discussed  in  Chapter  IV. 


J 


The  Derivation.  To  begin  the  procedure,  the  cross  product 


1 


between  Eq.  (38)  and  an  arbitrary  unit  vector  u is  computed.  Then  the 
surface  current  is  expressed  in  terms  of  its  rectangular  components, 
and  the  resulting  equation  is  projected  onto  the  x,  y,  axes  to 

obtain  the  following  three  equations  relating  the  rectangular  components 
of  R".  and  T. , 

J J 


• [u  X H(r)] 


m=l  . L 


a (J  d)  - J A ) 
y'  ym^x  xm^y' 


m 


"a(d  A “J  a)I  dS 


2'  xm"z  zm"x'  m 


3y  • [u  X H(r)] 


i^-f[ 

m=l  -4 . L 


a(J  A "**3  a) 


(40) 


a (J  A - J A ) dS' 


x'“ym-"x  xm^y'  m 


• [u  X H(r)]  = 


m=l  L 


a (J  A " J A ) 
X xm^z  zm^x' 


(41) 


- a (J  A - J A ) 
y'  zm^y  ym^z' 


dS' 

J 


where  the  subscript  j has  been  dropped  and  a„  = etc. 

X dX 


(42) 


From  these  equations,  which  aooly  to  resonators  containing 
mirrors  of  arbitrary  curvature,  it  is  evident  that  the  rectangular 
components  of  the  current  induced  on  the  mirrors  and  the  magnetic  field 
within  the  resonator  are  related  in  a very  comolicated  manner.  However, 


: j' 

I h 
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of  the  open  resonators  which  find  practical  application,  a large  fraction 
employs  mirrors  for  which 


— « 1 (43) 

|Rl 

where  a is  the  radius  and  R is  the  radius  of  curvature  of  the 

mirror.  For  these  mirrors,  the  surface  normal  is  nearly  parallel  to  the 

optic  axis,  and  the  component  of  current  parallel  to  that  axis  is  very 

small.  Since  the  effect  of  these  deviations  /rom  currents  which  are 

» * 

transverse  and  mirror  normals  which  are  parallel  to  the  optic  axis  is 
likely  to  be  small,  their  effect  will  be  neglected  for  the  remainder  of 
the  analysis. 

It  is  worth  noting  that  for  plane  parallel  resonators  there  is 

no  deviation  from  transverse  currents  or  longitudinal  mirror  normals. 

Although  this  fact  does  not  prove  anything  for  curved  mirrors  for  which 
a 

-jYj"  « 1,  it  does  support  the  idea  that,  for  these  mirrors,  the  longitu- 

dinal  component  of  the  current  and  the  transverse  components  of  the  mirror 
normals  are  negligible.  Then  settina  = 0 and  letting  u = ±a^in  Eqs. 

(40) through  (42),  the  following  pair  of  equations  is  obtained. 

• [li,  X H(F)]  f ds;  (44) 

m=l 

m 

2 

Sy  • C±a,  X H(F)]  = 7 f ds; 

m=l  s- 


(45) 


/ 


These  equations  are  identical  in  form  and  the  x and  y components  of  the 
current  are  not  coupled.  Then  without  loss  of  generality,  only  the 
x-component  of  the  current  will  be  treated  for  the  remainder  of  this 
paper. 

To  convert  Eq.  (44)  into  integral  equations  in  which  the  only 
unknowns  are  the  currents  on  the  mirrors,  the  field  point,  with  position 
vector  r,  will  be  allowed  to  approach  an  arbitrary  point  on  each  mirror 
surface  in  the  limit  as  the  distance  between  the  field  point  and  the  mirror 
goes  to  zero.  In  each  case,  the  approach  will  be  made  from  within  the 
resonator.  As  the  field  point  approaches  first  one  mirror  and  then  the 


other,  the  appropriate  sign  will  be  chosen  so  that  ±a^  corresponds  to 


the  unit  normal  at  the  point  of  interest  on  the  surface. 

Then  letting  the  field  point  approach  a point  with  position 
vector  r^  on  mirror  #1,  and  choosing  the  upper  sign  in  Eq.  (44),  Eq.  (46) 
is  obtained. 


-/ 


'^x2^'"2  ds; 


(46) 


r=r. 


where  Eq.  (37)  and  the  relation  = a^  • V4>  with  = a^  have  been 
employed.  Using  the  expression  for  <()  given  in  Eq.  (21),  the  normal 
derivative  of  (t  is  written  in  the  form  shown  below. 


an. 


(1  + ikR.|2)ni  • ^ 


(Al) 


(47) 
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I 

This  expression  is  used  in  Eq.  (46)  to  rewrite  the  limit  in  the 

1 form 


where 


u(r')  = e‘‘*'^’^12  (;  + ikR^^)  (4gy 

Using  Eq.  (50)  to  evaluate  this  limit  (Stakgold,  ref.  36,  vol . 2,  p.  119), 


and  substituting  the  result  into  Eq.  (46),  the  following  equation  re- 
lating the  currents  on  the  mirrors  is  obtained. 


S'  r=ro 


f I + ■ ’"pi 

S'  |_  l'^2  *^2' 


(|)(r2|r^)cos(r2  - r^,  n2)dS^  (52 


Classification  of  the  Modes.  Based  on  the  analysis  presented  on 
the  past  few  paqes,  the  problem  of  determining  the  modes  for  laser  reso- 
nators with  perfectly  conducting  mirrors  has  been  reduced  to  solving  the 
following  pair  of  integral  equations: 


where 


q)l  XJI 


'^xl  “ '^n'^xl  '^12'^x2 

(53 

“^xZ  " ^21^x1  hz^y^Z 

(54 

(55 

r,  f 1 ■'k|r  - rM 

K J = 2(-l)^/  J (r')  3_  ^(r'lr  ) 

qq  xq  ^ ’ J xq^  gM  i..  _ ..n  q'  q' 


xq'  q' 


r - r ' 


cos(r.  - r',  n„)dS' 
'■  q q’  q'  q 


The  solutions  of  this  pair  of  equations,  which  are  laser  resonator  modes, 
can  be  divided  into  the  three  classes  defined  below. 


1.  A Class  I mode  is  defined  to  be  a mode  for  which  the  self- 
induction  terms  in  Eqs.  (53)  and  (54)  are  neqliqible  in  com- 
parison to  the  mutual -induction  terms;  that  is 


Ik  J |«|J  |. 

I qq  xq'  ' xq' 

2.  A Class  II  mode  is  defined  to  be  a mode  for  which  the  self- 
induction  terms  are  proportional  to  the  mutual -induction  terms; 
that  is, 

qq  xq  q q*-  x)i 

3.  A Class  III  mode  is  any  mode  not  falling  into  Classes  I or  II. 
Since  the  modes  falling  into  this  last  class  are  beyond  the 
scope  of  this  paper,  they  will  not  be  considered  further. 

For  the  modes  falling  into  Class  II,  Eqs.  (53)  and  (54)  take  the 
form  shown  below. 


Jxi  = (T  + Y])  Kt2'^x2 


(57) 


'x2  = ^ ^2^  ^21^x1 

Substituting  Eg.  (58)  into  Eq.  (57)  and  vice  versa,  one  obtains  the 
familiar  eigenvalue  problems  of  the  form. 


(58) 


- 0 *Y,U1  (59) 

= (1  ♦ YjjO  * r,)  (60) 

which  apply  for  all  values  of  k such  that  k'>0. 

It  should  be  apparent  that  not  all  of  the  solutions  to  Eqs.  (59) 
and  (60)  are  solutions  to  the  more  general  set,  Eqs.  (53)  and  (54).  Thus, 
not  all  solutions  of  Eqs.  (59)  and  (60)  are  laser  resonator  modes.  For 
such  solutions  to  be  Class  II  modes,  they  must  also  satisfy  the  consistency 
conditions. 
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2*^x2 


'^22'^x2  " ^2'^21'^xl 


for  an  points  on  mirrors  #7  and  #2.  However,  due  to  the  difficulty  in- 
volved in  demonstrating  these  consistency  conditions,  this  potentially 
important  class  of  modes  will  not  be  considered  further  in  this  paper. 
For  Class  I modes,  Eqs.  (53)  and  (54)  take  the  form. 


■^xl  ■ "l2'’x2 


■ K2,J^,  (64) 

As  with  Class  II  modes,  one  can  obtain  an  eigenvalue  formulation  by  sub- 
stituting for  *^x2  and  vice  versa.  The  result  is  the  following  pair  of 
equations. 


'^xl  = ^12^21'^xl 


'^x2  " '^2l'^12'^x2 


However,  not  all  solutions  of  Eqs.  (65)  and  (66)  are  solutions  to 
Eqs.  (53)  and  (54).  Thus,  for  solutions  of  Eqs.  (65)  and  (66)  to  be 
laser  resonator  modes,  they  must  satisfy  the  consistency  conditions 


'^22'^x2l*'l'^x2l 


at  all  points  on  the  resonator  mirrors. 

In  appendix  C,  calculations  are  performed  to  estimate  the  condi- 
tions under  which  the  ineoualities  in  Eqs.  (67)  and  (68)  will  apply. 
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For  mirrors  such  that  a «|R|  and  a < j |l  estimate  (which 

was  obtained  by  specializing  Eqs.  (67)  and  (68)  to  apply  to  azimuthally 
symmetric  modes  at  the  center  of  circular  mirrors)  corresponds  to  the 


conditions  shown  below. 


‘“I 


1 * K fx 

2TRT®  4 


/l  + (J)2  + |-  a «1 


where  6g  = - l'.  Other  conditions,  which  correspond  to  cases  where 


a > j y given  in  appendix  C. 

Not  surprisingly,  these  estimates  of  the  consistency  conditions 
yield  rather  complex  results  that  depend  on  X,  a,  jR],  L,  and  the  values 
of  k"  for  each  candidate  mode.  For  ranges  of  these  parameters  where  these 
conditions  are  marginal,  the  rigorous  conditions  of  Eqs.  (67)  and  (68) 
should  be  applied. 

Then  assuming  that  the  inequalities  in  Eqs.  (69)  and  (70)  hold, 
the  equations  for  Class  I modes  take  the  form 


~ikR2i 


COSQ21  dSj 


'^xl^'l^  = 


'^x2^’"2^  R^2  ^°^“12  ^^^2 
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where  it  has  been  assumed  that  |k|»l/R^-j  and  is  the  angle  between 
R^j  and  the  optic  axis.  Finally,  substitutina  Eq.  (72)  into  Eq.  (71) 
and  vice  versa,  one  obtains  the  following  pair  of  symmetric  integral 
equations  for  the  currents  on  the  mirrors. 


'^xl^’^l^  " brj  / '^xl^l 


-i  k( R^  2 R2] ) 


f 

■/,  ^12  *^21 


cosa-|2  cos32^  d$2  dS-j 


'^x2^'"2^  " 


Cl  Cl  ^ * 


-ik(R2^  + R^2^ 


COSa2i  COSa^2 


This  last  pair  of  equations  bears  considerable  resemblance  to  the 
integral  equations  that  are  normally  used  (Fox  and  Li,  ref.  1,  p.  454)  to 
analyze  the  modes  of  laser  resonators.  However,  three  important  dif- 
ferences are  listed  below. 

1.  Eqs.  (73)  and  (74)  are  not  in  the  form  of  a linear  eigenvalue 
problem.  In  this  paper,  the  eigenvalue  k is  not  only  a multi- 
plier, but  it  is  also  included  in  the  integrand. 

2.  The  obliquity  factor,  cosa^2»  paper  is  different 

from  the  factor,  1/2  (1  + cosa^2)»  which  is  normally  used 
(Fox  and  Li,  ref.  1,  p.  454). 
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3.  Eqs.  (73)  and  (74)  apply  to  the  current  distribution  in- 
duced on  perfectly  conducting  mirrors,  while  the  equations 
normally  used  are  widely  considered  to  apply  to  the  field 
(electric  or  magnetic)  on  resonator  mirrors  of  any 
material . 

These  first  two  differences,  which  correspond  to  differences  in  the  form 
of  the  equations,  can  have  a significant  effect  on  the  mode  distributions 
and  losses.  However,  for  some  of  the  modes  of  paraxial  resonators,  the 
effects  of  these  two  differences  in  form  are  negligible.  Although  the 
third  difference  will  have  no  effect  on  the  forms  of  the  solutions  ob- 
tained, it  does  affect  the  physical  interpretation  of  the  terms  in  the 
equations  as  well  as  the  conditions  under  which  the  equations  can  be 
applied.  Finally,  it  is  important  to  note  that,  (1)  not  all  solutions  of 
Eqs.  (73)  and  (74)  are  resonator  modes;  and  (2)  these  modes  are  not  the 
only  modes  of  the  resonator.  Other  possible  modes  are  Class  II  and  III 
laser  resonator  modes  and  the  entire  set  of  incoming  wave  modes  (which 
have  not  been  addressed  in  this  paper). 

In  any  case,  to  determine  the  Class  I modes  of  a laser  resonator 
with  perfectly  conducting  mirrors,  a procedure  must  be  developed  for 
solving  Eqs.  (73)  and  (74).  That  procedure,  which  must  specify 
simultaneously  both  the  real  and  imaginary  parts  of  k as  well  as  the 
current  distributions  on  the  two  mirrors,  is  the  subject  of  the  next 
chapter. 
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solving  the  integral  equations  derived  in  the  previous  chapter  and  to 
present  the  method,  based  on  a variational  principle,  to  be  used  to 
obtain  approximate  solutions  to  these  equations  for  open  resonators 
satisfying  the  paraxial  approximation.  To  accomplish  these  objectives, 
this  chapter  begins  with  a discussion  of  the  general  form  of  the  inte- 
gral equations  and  the  relationships  between  the  solutions  to  these 
equations  and  the  fundamental  resonator  parameters.  The  paraxial 
approximation  is  then  discussed  and  applied  to  the  general  equations  to 
obtain  an  eigenvalue  formulation  for  paraxial  resonators.  Following  the 
formulation,  the  variational  method  for  obtaining  approximate  solutions 
to  these  equations  is  presented  and  discussed. 

The  Basic  Problem 


As  indicated  at  the  end  of  the  preceding  chapter,  the  first  step 


in  determining  the  modal  properties  of  open  resonators  is  solving  Eqs. 
(73)  and  (74).  To  solve  these  eouations,  the  quantities  , J^2* 
k must  be  specified  simultaneously  for  each  mode.  Once  k is  known,  the 
resonator  mode  and  material  parameters  may  be  determined  using  Eq.  (22), 
which  is  repeated  without  subscripts  below. 


k = + 7 + x)  - — (22) 

Writing  the  square  root  in  series  form,  assuming  that  ^Ix  - <<  U 

and  writing  k,  u,  and  x "in  terms  of  their  real  (')  and  imaginary  ('') 
parts,  the  following  relationships  are  obtained. 
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Using  these  relations  with  the  computed  values  of  k'  and  k"  (which  have 
been  obtained  as  part  of  the  solution  of  the  integral  equations)  and  the 
given  values  of  four  of  the  parameters  u',  u",  x'»  x"»  n,  and  a,  the 
remaining  two  quantitites  can  be  determined. 

The  main  difficulty  with  the  general  approach  just  described 
lies  in  obtaining  the  solutions  to  Eqs.  (73)  and  (74).  This  is  espe- 
cially true  when  information  concerning  higher  loss  modes  is  desired. 

One  possible  approach  to  obtaining  such  solutions  begins  with  an  itera- 
tion approach  similar  to  that  used  by  other  workers  (refs.  1,  26)  in  the 
field.  By  applying  that  procedure  to  Eqs.  (73)  and  (74),  one  can  deter- 
mine the  current  distributions  and  the  value  of  k for  the  lowest  loss 
mode.  To  extend  that  procedure  to  the  next  lowest  loss  mode,  one  must 
insure  that  each  successive  approximation  for  that  mode  is  orthogonal  to 
the  solution  for  the  lowest  loss  mode.  Similarly,  the  approximation  for 
each  higher  loss  mode  must  be  made  orthogonal  to  all  of  the  lower  loss 
modes  that  have  been  obtained. 

This  iteration/orthogonal ization  procedure  is,  in  general,  very 
difficult  and  time  consuming  to  apply.  However,  as  shown  in  the  next 
section,  a considerably  simpler  procedure  is  available  for  open  resonators. 


Formulation  of  the  Eigenvalue  Problem  for  Paraxial  Resonators 


For  many  resonators  containing  mirrors  for  which  a/|R|<<l, 


the  mirror  separation  L is  much  greater  than  the  mirror  radius  a.  Such 
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resonators  satisfy  the  paraxial  approximation.  Some  of  the  modes  of 

these  paraxial  resonators  satisfy  the  more  stringent  condition  that, 

for  the  mode,  the  quantity  e^*^j*^12  is  accurately  approximated  by 
+k'-L 

e J for  any  two  points  on  the  mirror  surfaces.  For  these  paraxial 
modes,  the  coupled  integral  equations  take  the  form  below. 


-ik'(R^2  '^21^c, 


dS^  [ dS ■ 


-ik  (^21  ^12^ 


dS'  dS^ 


For  resonators  operating  at  or  near  optical  wavelengths,  the 
relations  k'»k"  and  a)'»a)"  hold  for  all  but  the  lossiest  modes.  Using 
these  relations,  Eq.  (77)  is  rewritten  in  the  form 


where 


■ -ik'(R^2  ^21^  ... 

e dS2 


Using  this  symmetric  kernel,  which  is  not  a function  of  k",  Eq.  (79)  is 
now  cast  in  the  form  of  a linear  eigenvalue  problem  with  an  external 
parameter  k'. 
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(81) 


'i 

Then  subject  to  the  constraint  of  finding  the  proper  value  of  k'  such 
that  k"  = ^ In  Y.  where  y is  real  and  positive,  the  problem  of  deter- 
mining the  modes  has  been  reduced  to  a linear  eigenvalue  problem. 

When  the  precise  frequency  spectrum  is  not  required,  this 
constraint  on  y can  be  relaxed  without  significantly  affecting  the 
current  distributions  or  the  damping  rates  for  the  modes  of  interest  in 
this  analysis  (relatively  low  loss  modes  at  optical  or  infra**ed  wave- 
lengths). This  relaxation  is  a result  of  the  relatively  high  mode 
density  which  is  characteristic  of  even  lossy  resonators  at  the  wave- 
lengths of  interest.  This  high  mode  density  limits  the  difference 
between  an  arbitrarily  chosen  value  of  k'  and  one  which  will  produce  a 
real,  positive  value  of  y.  For  example,  the  axial  mode  spacing  of  a 
plane  parallel  resonator  in  free  space  with  mirror  separation  L is 

Af  = C/2L  (82) 

a 

As  the  maximum  shift  off  a resonance  is  one-half  of  this  spacing,  the 
maximum  shift  or  error  in  the  value  of  k'  (at  threshold  and  line  center) 
is  ir/2L.  This  maximum  shift  corresponds  to  a relative  change  in  the 
resonator  Fresnel  number,  N = a^/xL,  of 

AN  = a2/4L2  (83) 


where  a is  the  mirror  radius. 


For  paraxial  resonators  of  all  but  the  largest  magnification, 
this  small  change  in  the  Fresnel  number  will  have  negligible  effect  on 
the  modal  parameters.  As  a result,  the  requirement  that  y must  be  real 
and  positive  will  be  dropped.  This  reduces  the  analysis  to  solving  an 
eigenvalue  problem  (with  complex  eigenvalue)  for  any  value  of  k'  corres- 
ponding to  optical  or  infrared  wavelengths.  As  complex  values  of  y will 
now  be  allowed,  the  expression  for  k"  must  now  be  written 

k"  = ^lnH  (84) 


For  certain  cases,  the  problem  can  be  simplified  further.  Those 
cases  correspond  to  problems  for  which  the  spatial  dependence  of  the 
kernel,  ) » can  be  separated  or  the  dependence  on  one  coordinate 

can  be  solved  in  closed  form.  In  the  first  case,  which  corresponds  to 
resonators  with  rectargular  mirrors,  the  separation  of  variables  reduces 
to  solving  two  eigenvalue  problems  as  shown  below. 


u(xi)  = y]J u(x])U|^,(x^  Ixpdx] 

(85) 

v(yi)  = 

(86) 

with  ^ = lYlYgl’  '^xl^’^1^  ^ u(x^)v(y^),  and  K,^,(r^  Irp  = 

U,^(xi  |xpVk.(yi  |yj). 

In  the  second  case,  v/hich  corresponds  to  mirrors  of  circular 
cross  section  (where  the  azimuthal  deoendence  can  be  solved  exactly), 
a single  eigenvalue  equation  involving  only  the  radial  coordinate 
results.  The  corresponding  equation  has  the  form 
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+2k'*L 

with  e = IyI  as  before. 


The  Variational  Method  of  Solution 


To  see  how  a variational  approach  can  be  applied  to  these  eigen- 
value problems,  it  is  useful  to  determine  the  variational  principle  (or 
problem)  satisfied  by  the  eigensolutions  of  the  type  under  consideration. 
To  determine  that  principle,  one  first  considers  eigensolutions  y and 


u(x)  which  satisfy 


Yu(x)  = J u(x')K(xIx')dx' 


where  K{x)x')  is  symmetric  with  respect  to  x and  x',  and  y = y"^  • 
then  multiplies  Eg.  (88)  by  an  arbitrary  function  w(x)  and  intearates 
the  result  over  x to  obtain  Eg.  (89). 

Y J*u(x)w(x)dx  = J*w(x)  J u(x' )K(x|x' )dx'dx  (89) 


Next,  one  varies  the  guantities  y»w(x),  and  u(x),  uses  the  fact  that 
K(x|x')  is  symmetric,  and  recombines  terms  to  obtain  the  result  shown 
in  Eg.  (90). 

J 6w(x)  |yu(x)  - J K(xlx' )u(x' )dx'|  dx  + J 6u(x)  |yw(x) 

XX  X 

x' )K(x|x' )dx' I dx  = - (Sy y*u(x)w(x)dx  (90) 


I 


Letting  w(x)  = u(x)  in  this  equation,  it  can  be  seen  that  if  Ju^(x)dx.  t 0, 
then  the  eigensolutions  7 and  u(x)  satisfy  the  condition 


6y  = 0 


where 


y*u(x)  yK(xlx' )u(x' )dx'dx 
— _ X x' 

y*u^(x)dx 


Thus  the  eigenfunctions  u(x)  are  those  functions  which  make  7 

stationary  for  all  variations  about  the  functions  u(x).  This  implies 

that  if  one  were  to  substitute  all  possible  trial  functions  for  u(x) 

in  Eq.  (92),  by  writing  u(x)  in  terms  of  an  infinite  number  of  parameters 

{A„},  7 would  have  stationary  values  for  those  values  of  the  {A„>  which 
n n 

would  yield  the  eigenfunctions  u(x).  One  cannot  in  practice  generate 
all  possible  trial  functions,  but  one  can  generate  all  possible 
variations  within  a given  class  of  functions  and  then  require  that 
67  = 0.  To  the  degree  that  the  chosen  class  has  the  capability  to 
represent  the  modes,  one  can  obtain  accurate  approximations  to  these 
modes.  One  well-known  technique  for  obtaining  such  stationary  approxi- 
mations is  the  Rayleigh-Ritz  procedure  (Morse  and  Feshback,  ref.  37, 
vol . 2,  p.  1118)  described  in  the  following  paragraphs. 

The  Rayleiah-Ritz  Procedure.  To  apply  the  Rayleigh-Ritz  pro- 
cedure to  the  eigenvalue  problem  shown  in  Eq.  (88),  one  begins  with 
the  variational  principle  just  described.  One  then  expands  u(x)  in 
terms  of  a set  of  known  functions  (lii^. (x)}. 
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u(x)  = LA.<|;.(x) 


T 


(93) 


and  substitutes  the  result  into  Eq.  (92)  to  obtain  the  followinq 
expression. 


i 


)'l'n,(x)dx 


2 f ij)^.(x)  rK(x|x')i(j^(x')dx'dx  (94) 

i ,m  J J , 


The  next  step  in  the  procedure  is  to  find  those  values  of  the 


independent  set  {A^}  which  make  y stationary.  It  can  be  shown  that  those 


values  satisfy  (Morse  and  Feshback,  ref.  37,  vol . 2,  p.  1119) 


_ n 
JK,  - ° 


(95) 


for  each  value  of  i.  Applying  that  condition,  the  followinq  set  of 
equations  is  obtained  for  i = 1,2,  — ,n, 


m m 


(96) 


where 


’ml 


(97) 


and 


/ K(x|x'  )ijj^. (x'  )dx'dx 

\J  w * 


(98) 
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This  set  of  homogeneous  equations  corresponds  to  the  matrix 


problem 


YfA  = KA  (99) 

where  y and  K are  square  matrices  with  elements  f . and  K . as  defined 

mn  mi 

above  and  A is  a column  matrix  with  elements  A^ . 

In  select! no  the  variational  method  to  solve  the  linear  eigen- 
value problem,  consideration  was  oiven  to  the  fact  that  the  eigenfunctions 
can  be  normalized  so  that  they  obey  the  orthogonality  condition  (Siegman 
and  Miller,  ref.  28,  p.  2730), 

/u„(x)u^(x)<lx  . (100) 

X 

Then  using  the  Rayleigh-Ritz  procedure,  stationary  approximations  to  the 

eigenvalues  7,  which  yield  approximations  to  the  orthogonal  eigenfunctions, 

can  be  found  if  the  matrix  elements  v . and  K . can  be  evaluated.  These 

mi  mi 

approximations  can  also  be  made  to  obey  Eq.  (100)  by  using  the  procedure 
discussed  in  appendix  R to  solve  the  matrix  problem. 

The  General  Form  of  the  Matrix  Elements.  To  evaluate  the  matrix 
elements,  it  will  be  necessary  to  know  the  specific  form  for  the  ip^-(x^) 
and  the  kernel  K(x^lxj).  Although  these  precise  forms  are  covered  in  the 
following  chapter,  the  basic  equations  are  summarized  below. 

Referring  to  Eq.  (80),  the  kernel  K(x.||xp  can  be  written  in  the 
general  form 

K(Xi|xp=y*  |x2)K2i(x^  |x2)dx2  (101) 

X2 
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where  the  integration  over  X2  is  taken  over  the  space  corresponding 
to  the  second  mirror.  Substituting  Eg.  (101)  into  Eg.  (98)  yields 


r 


c 

\ . 

t- 


To  simplify  the  notation,  the  function  ^^(^2)  is  defined  according  to 
Eg.  (103)  below. 


~ ^ I^12(^"|I^2^  '^m^^l  ^^^1  (103) 

^1 

Using  this  definition,  the  matrix  elements  can  be  compactly  summarized 
as  shown  below. 

./dx^Ajx^)  Ai(x2>  (104) 

*2 

■'ml  ■/  '‘"l  ♦m<’‘l>  ♦l<’‘l>  <’°5> 

^1 

As  previously  indicated,  the  next  step  in  the  procedure  for 
obtaining  the  modes  of  open  resonators  is  the  application  of  these  general 
expressions  to  the  particular  resonators  of  interest.  This  specializa- 
tion is  the  subject  of  the  next  chapter. 

However,  before  proceeding  to  that  chapter,  it  may  be  helpful 
to  summarize  the  types  of  resonators  and  modes  for  which  the  basic 
theory  and  method  of  solution  just  outlined  are  applicable.  First, 
the  basic  eguations  themselves  (Egs.  (79)  and  (80)  are  applicable  only 
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for  resonators  with  perfectly  conducting  mirrors  where  the  effect  of  any 
gain  is  essentially  constant  for  all  points  on  the  resonator  mirrors. 

In  addition,  these  equations  are  applicable  only  for  resonators  for 
which  the  oscillation  wavelength  is  small  in  comparison  to  all  resonator 
dimensions  and  where  the  mirror  diameters  are  small  in  comparison  to 
the  mirror  separation  and  radii  of  curvature.  Finally,  these  equations 
are  applicable  only  for  cases  in  which  the  self-induction  integrals  in 
Eqs.  (53)  and  (54)  are  negligible,  and  the  method  of  solution  should  be 
employed  only  for  cases  that  do  not  require  knowledge  of  the  precise 
frequency  spectrum. 


. Application  of  the  VariationaT  Method 

The  purpose  of  this  chapter  Is  to  specialize  the  basic  theory 
and  method  of  solution  to  open  resonators  of  rectangular  or  circular 
projection  for  which  the  spatial  dependence  separates.  These  speciali- 
zations Include  obtaining  the  specific  form  of  the  Integral  equation  as 
well  as  the  expansion  functions  for  each  geometry. 

Since  It  Is  to  be  applied  to  several  cases,  the  general  pro- 
cedure for  obtaining  these  specializations  Is  the  first  subject  of 
discussion  In  this  chapter.  This  procedure  Is  then  applied  to  resona- 
tors with  spherical  mirrors  of  circular  projection;  and  to  Illustrate 
Its  various  aspects,  the  calculations  are  covered  In  some  detail. 

Finally,  the  chapter  ends  with  a similar  section  for  rectangular  mirror 
resonators;  however,  the  details  for  this  case  are  covered  In  appendix  D. 

The  General  Procedure 

Specific  Forms  of  the  Integral  Equations.  The  procedure  for 
obtaining  the  integral  equation  Is  usually  relatively  simple.  To  apply 
It,  one  begins  with  Eqs.  (71)  and  (72)  as  the  basic  forms.  One  then 
makes  the  following  assumptions: 

1.  The  mirror  radii  are  much  smaller  than  the  mirror 
radii  of  curvature, 

2.  The  paraxial  approximation  holds, 

+k"R,2  +k"L 

3.  e lae  across  the  mirror  surfaces, 

4.  The  precise  frequency  spectrum  Is  not  of  interest, 

5.  k'»k". 
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Using  these  assumptions  and  substituting  Eq.  (72)  into  Eq.  (71)  and  vice 

versa,  one  obtains  Eqs.  (80)  and  (81)  from  chapter  IV.  Then,  the 

distance  between  two  arbitrary  points  on  the  mirrors,  which  is  expressed 

in  terms  of  the  resonator  constants  and  transverse  mirror  coordinates, 

is  written  in  series  form  using  a binomial  expansion.  In  this  expansion, 

a sufficient  number  of  terms  is  kept  to  accurately  approximate  the 
-ik'R 

exponential  e across  the  mirror  surfaces.  Finally,  where  possible, 
the  spatial  variation  of  the  equation  is  separated  to  yield  one  or  more 

equations  of  the  type  shown  in  Eqs.  (85)  through  (87). 

Selection  of  Expansion  Functions. 

Desirable  Characteristics--A1thouqh  the  precise  forms  are 
somewhat  arbitrary,  the  expansion  functions  used  to  calculate  the  reso- 
nator modes  should  have  two  desirable  characteristics.  First,  one 
should  be  able  to  establish  the  precise  form  of  the  functions  from  the 
resonator  parameters  with  relative  ease.  Second,  the  functions  should 
exhibit  (to  whatever  degree  is  practical)  the  behavior  expected  from  the 
actual  modes  so  that 

1.  Trends  can  be  identified  with  respect  to  changes  in 
parameters , 

2.  Physical  insight  can  be  gained  into  the  behavior  of 
the  modes, 

3.  The  modes  can  be  adequately  represented  using  a rela- 
tively small  number  of  expansion  functions. 

Excluding  the  actual  modes  themselves,  three  categories 

that  exhibit  some  of  these  desired  characteristics  are  listed  below; 

1.  Approximate  solutions  of  the  derived  resonator 
integral  equation, 

2.  Exact  solutions  of  the  derived  integral  equation  for 
a resonator  similar  to  the  one  being  considered. 
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3.  Approximate  solutions  of  the  derived  integral  equation 
for  a resonator  similar  to  the  one  being  considered. 

Unfortunately,  to  obtain  the  functions  corresponding  to  the  first  two 

cases,  one  must  solve  a problem  which  is  at  least  as  difficult  as  the 

original.  As  a result,  the  expansion  functions  will  be  chosen  ta  be 

approximate  solutions  of  an  integral  equation  for  a resonator  similar  to 

the  one  being  considered.  The  "similar"  resonator  chosen  for  each  case 

is  a symmetric  resonator,  which  consists  of  two  mirrors  such  that 

1.  The  mirrors  are  identical  to  the  one  on  which  the 
current  is  being  analyzed,  and 

2.  The  mirror  separation  is  the  same  as  for  the  original 
problem. 

Thus,  to  determine  the  expansion  functions,  approximate 
solutions  of  integral  equations  of  the  following  form  must  be  obtained, 

= Y K(jr|)T')  djT'  (106) 

x' 


K(x|x') 


(107) 


where  the  subscript  x denoting  the  x-component  of  the  current  has  been 
dropped,  and  N()T|)r')  is  symmetric  with  respect  to  the  spatial  coordinates 
X and  x' . As  indicated  before,  when  the  spatial  dependence  separates, 

Eq.  (106)  reduces  to  one-dimensional  equations  with  symmetric  kernels. 

Possible  Approaches  for  Obtaining  the  Approximate  Solutions-- 
To  discuss  the  procedure  for  approximately  solving  these  symmetric  inte- 
gral equations,  Eq.  (106)  is  shown  below  in  operator  notation. 

J(>T)  = yK  J(x)  (108) 
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For  the  first  approach,  one  assumes  that  an  operator  M can 

A A 

be  found  such  that  M commutes  with  K.  Then,  by  definition,  Eq.  (109) 
holds, 

M K J(x)  = K M J(3T)  (109) 

A A 

and  the  operators  K and  M have  simultaneous  eigenfunctions.  This  implies 
that  if  the  eigenfunctions  of  the  operator  M can  be  found,  the  eigen- 
functions of  K can  be  expressed  as  a linear  combination  of  them.  Thus, 

if  functions  v (x)  can  be  found  such  that 
s 

M v^()D  * -s2Vg(7)  (110) 

then  the  eigenfunctions  0(7)  of  K can  be  expressed  in  the  form 

The  particular  linear  combination  to  be  used  in  each  case  can  be  deter- 
mined using  the  variational  procedure  discussed  in  chapter  IV  to  specify 

the  A .. 
ys 

Unfortunately,  for  all  but  the  simplest  cases,  finding  an 

A A A 

operator  M,  and  its  eigenfunctions,  so  that  M commutes  with  K is  very 
difficult.  However,  in  some  cases  it  may  be  possible  to  find  an 

A 

A A 

operator  M such  that  the  commutator  between  M and  K is  relatively  small, 
especially  for  low  loss  modes.  For  such  cases,  Eq.  (109)  must  be  modi- 
fied to  include  this  relatively  small  difference,  which  is  represented 
by  RJ(7)  in  each  of  the  following  equations, 

M i<  J(>T)  = K M J(x)  - R J(x)  (112) 

M J()D  = yK  M J(7)  - yR  J()T)  (113) 
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where  Eq.  (113)  was  obtained  by  combining  Eq.  (108)  with  Eq.  (112).  For 

A A 

this  case,  M and  K do  not  have  simultaneous  eigenfunctions,  but  if 
|R  J(x)|  is  indeed  small  in  comparison  to  |K  M J(30|.  the  eigenfunctions 
of  M may  still  form  an  excellent  expansion  set  for  the  eigenfunctions  of 
K. 

Unfortunately,  it  may  be  quite  difficult  to  actually  show 
that  Eq.  (114)  holds. 

|R  J(>0|  « |K  M J()T)|  (114) 

However,  even  if  one  cannot  find  operators  M and  R such  that  Eq.  (114) 
can  be  shown  to  hold,  it  may  still  be  possible  to  find  forms  for  M and 
R such  that  the  eigenfunctions  of  M form  an  excellent  expansion  set  for 
the  eigenfunctions  of  K.  With  that  in  mind,  the  following  rather  simple 
procedure  is  used  to  obtain  the  approximate  solutions  to  Eq.  (108). 

The  Chosen  Procedure— First,  using  the  operations  discussed 
in  the  following  paragraphs,  Eq.  (108)  will  be  manipulated  to  yield  an 
equation  of  the  form  shown  in  Eq.  (113).  The  particular  manipulations 
used  will  be  selected  to  yield  an  operator  M such  that  one  can  solve  the 
eigenvalue  problem  shown  in  Eq.  (110).  In  addition,  the  manipulations 
will  be  chosen  so  R J()T)  contains  only  terms  which  either 

1.  Depend  on  J(jT)  and  its  derivatives  only  at  the  mirror 
edges,  or 

2.  Involve  products  of  the  coordinates  of  both  resonator 
mirrors . 

When  Eq.  (110)  has  been  solved,  the  functions  v^(x)  will  be  chosen  as  the 
approximate  solutions  of  the  symmetric  integral  equation  shown  in 
Eq.  (108).  For  mirrors  which  are  conic  sections,  this  selection  will 
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yield  expansion  functions  that  are  obtained  by  neglecting  the  finite 
sizes  of  the  resonator  mirrors. 

The  operations  used  to  obtain  Eq.  (113)  from  Eq.  (108)  are 
similar  to  those  used  by  Bergstein  (ref.  2,  p.  497)  and  others 
(refs.  22,  23).  These  manipulations  are  outlined  below  for  resonators 
such  that 

1.  The  mirrors  are  conic  sections,  and 

2.  The  integral  equation  separates. 

The  approach  used  for  extending  this  procedure  to  toroidal  mirrors  is 
covered  in  appendix  D. 

One  begins  by  separating  Eq.  (106)  to  yield  one  or  more  of 
the  equations  of  the  form  of  Eqs.  (85)  through  (87),  which  are  repre- 
sented by  Eq.  (115) , 

u(t)  ~ y J u(c')  K(?l!;')  dc'  (115) 

where  ? = x,  y,  or  p.  One  then  computes  the  second  derivative  of 
Eq.  (115)  to  obtain  Eq.  (116). 

= Y f uU')  015) 

d?2  dc2 

Next,  one  integrates  the  term  K( ?]!;')  ^ by  parts  twice  to  yield 

dc'2 


y f K(cU')  dg  ■ = yfuUn  d^'  + yR'uU')  (117) 

dc’2  J,  d?'2 
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where  the  expression 


R'u(i;') 


K(cU') 

dc' 


u(c') 


dK(gk') 

d?' 


018) 


is  evaluated  at  the  corresponding  limits  of  integration.  One  then  sub-  ] 
tracts  Eq.  (117)  from  Eq.  (116)  to  obtain  the  following  rather  compli-  | 
cated  relation.  ] 


dfuM.  ^rK(5|c')  - rR'u(t') 

d?2  J d?'2 

C' 


+ yJ  u(c') 
c' 


d?2 


dc'2 


K(ck')  dc' 


(119) 


For  symmetric,  separable  kernels  of  the  type  shown  in  Eq.  (107), 
the  derivatives  of  K(c|c')  in  Eq.  (119)  can  be  written  as  the  sum  of 
three  terms  as  shown  below. 


K(tU')  (120) 


where 

1.  The  function  h(g,;c')  contains  all  terms  involving  products 
of  c and  c'  (the  coordinates  of  both  resonator  mirrors),* 

2.  The  function  f(g,c)  contains  all  terms  involving  c but 
not  c' , and 

3.  The  function  f(g,;')  contains  all  terms  involvina  c' 
but  not  ?. 

*For  cases  where  the  mirrors  are  conic  sections,  h(g,cc')  vanishes 
(Bergstein,  ref.  2,  p.  497). 


d?2 


d2 


dc 


i2 


K(ck')  = 


- f(g.c)  + f(g,t')  + h(g,cc') 


I 

! 

i 

I 
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Then,  substitutinq  Eq.  (120)  into  Eq.  (119),  and  using  Eq.  (115)  to 


show  that 


yj  f(g,c)  K(c|c')  u(c')dc'  = ^(q.c)  u(?) 


Eq.  (119)  can  be  rewritten  in  the  form  shown  below. 


f (a,c)  u(c)  = y/'k(5U')  + f(g.c')  u(c')d!;' 

•i,  (dc'2 


+ y/  K(cU')  h(q,ci:-)u(c')dc’  - yR'u(;').  (122 

c' 

Examining  the  last  two  terms  in  Eq.  (122),  one  can  see  that 

1.  R'u(c')  depends  on  u(c')  and  its  derivatives  only 
at  the  limits  of  integration,  and 

2.  The  integral  involves  products  of  both  resonator 
coordinates. 

Then,  subject  to  the  condition  that  u(c')  must  be  selected  so  that 
R'u(c')  depends  on  u(c')  and  its  derivatives  only  at  the  mirror  edges, 
one  identifies  Ru(c')  by  comparing  Eq.  (122)  with  Eq.  (113).  The  result 


Ru(c')  = R'u(c')  - y K(cU')  h(a,cc')u(c')dc' 

i' 


Substituting  this  form  into  Eq.  (122),  one  obtains  the  following  result. 
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+ i"(g.?) 


u(!;)  = -yRu(c’) 


f K(clc') 

+ f(q.c') 

f 

u(5')d?'  (124) 


Specializing  this  equation  to  resonators  containing  mirrors  which  are 
conic  sections,  and  comparing  the  result  with  Eg.  (113),  the  operator 
M is  identified. 


•'  h2 

M = + f(g,c) 

dt^ 


(125) 


and  the  eigenvalue  problem  to  be  solved  has  the  form  shown  below. 


: j 

■ 

■I. 

: 

r i 


+ f(q,4)  u(?)  = -s2u(?)  (126) 

d?2 

The  eigenfunctions  obtained  by  solving  Eg.  (126)  will  be 
too  general  to  be  applied  to  the  resonator  problem  without  restrictions. 
Three  appropriate  restrictions  are  listed  below. 

1.  For  rectangular  mirror  resonators,  the  modes  (and  the 
expansion  functions)  are  either  even  or  odd  functions. 

2.  The  expansion  functions  must  be  finite  at  all  points 
'on  the  mirrors. 

3.  The  expansion  functions  must  be  consistent  with  other 
known  behavior  for  the  modes. 

For  unstable  resonators,  the  last  condition  corresponds  to 

the  requirement  that  as  q->l , the  eigenfunctions  of  M must  reduce  to  the 
known  forms  for  plane  parallel  resonators.  These  known  forms,  which 
were  obtained  by  Vainshtein  (ref.  21,  p.  711)  using  the  Waveguide 


i 
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Analogy  and  which  are  to  be  used  as  expansion  functions  for  plane 
parallel  resonators,  are  listed  below. 

For  rectangular  mirrors  of  width  2a  and  seoaration  2L, 
the  expansion  functions  for  modes  of  odd  symmetry  are 

f (x)  = sin  ^ ^;n=  2,4,6--  (127) 

" 2a  {1  + (1  + i)f} 

and  the  expansion  functions  for  modes  of  even  symmetry  are 


COS 


nirx 


2a  {1  + (1 


;n  = 1 ,3,5  — 


(128) 


5(7)  

In  these  expressions,  3 = w 0.824,  M = /8ttN,  and  ^ is  the  Riemann 

zeta  function  (Erdelyi,  ref.  41,  vol . 2,  p.  32,  Eg.  2). 

For  circular  mirrors  of  radius  a and  separation  2L,  the 
expansion  functions  are  given  by  • 


^ h 'f*  h 

where  v is  the  m^  root  of  the  n order  Bessel 
nm 

kind,  Jp(x),  and  3 and  M are  given  above. 


(129) 

function  of  the  first 


Resonators  with  Spherical  Mirrors  of  Circular  Projection 

To  aoply  this  procedure  to  a circular  mirror  resonator  with 
mirrors  of  radii  a-|  and  beoins  with  Eos.  (71)  and  (72).  The 

appropriate  forms  of  these  equations,  which  express  the  current  on 
mirror  #1  in  terns  of  the  current  on  mirror  ^2  and  vice  versa,  are 
shown  below. 
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Jx2(p2’®2^  " ■ ^ J J 


R2^  cosa2iP^dp^de^ 


0 0 


2i:  a. 


i k r r“^  e * 

^xi^Pi*®])  ^ ~^J  J ^x2^'^2’®2^  ~~^  cosa^2P2‘^P2'^®2 


0 0 


where  R^2  "'s  given  by  Eq.  (132) 


R-| 2 ~ ^2  ~ ^*^lP2  ^®2  ~ ®1  ^ "*”  ^^2  ~ 


(132) 


To  obtain  the  expansion  functions  for  the  current  on  mirror  , 
one  specializes  Eq.  (131)  to  apply  to  a symmetric  resonator  containing 
two  mirrors  identical  to  mirror  #1  with  the  same  mirror  separation 
present  in  the  original  problem.  Since  the  resonator  being  considered 
(to  obtain  the  expansion  functions)  is  symmetric,  ai  = a2>  R-|  = R2, 


Jx2(P2»®2^  = e ' Jx1^^2’®2^ 


where  n'  = uq,  and  q is  an  integer.*  Thus,  the  equation  for  the 
current  on  mirror  #1  of  this  symmetric  resonator  can  be  written  in 
the  following  form. 


(133) 


Oxi(pi .^i) 


^7T  d-j  -1KK-|2 

= -^e  /*  ^®2 /*  '^x1^^2’®2^~r7T  cosa.|2p2dp2  (134) 

•'0 

*0ne  can  see  this  by  choosing  the  origin  midway  between  the  mirrors 
and  realizing  that  the  fields  must  be  even  or  odd  with  respect  to  z. 
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To  simplify  the  procedure  for  obtaining  these  expansion 
functions,  the  following  approximations  are  made.* 


(136) 


k"  « k' 


(137) 


a^  « L 


(138) 


Under  these  cond.tions,  cosa.|2«1  and 

^-ik'R^2  _ik'L  '’^U  '^^1^  ^2^^  '^^'^1^2  ‘^°^(®2  " ®1  ^ 


where  g = 1 - and  the  first  two  terips  of  the  binomial  expansion  have 
been  used  to  approximate  R^2*  Finally,  it  is  assumed  that 


k-,2  +k"L 
e « e 


(140) 


across  the  entire  mirror  surface.  With  these  approximations,  Eg.  (134) 
can  be  manipulated  to  yield 


where  Eos.  (135)  throuah  (138)  are  not  valid,  usina 
r.-r*'’ijci~atlons  will  simply  require  the  use  of  a relatively  large 
•*~anslon  'unctions  to  represent  the  modes. 
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+-j-p,02  cos(92  - 6,) 
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041) 


where  y = -e 


+in'-ikL 


Next,  one  assumes  solutions  of  the  form. 


Jj^l(p.e)  = u^(p)  e 


and  applies  the  identity  (Erdelyi,  ref.  41,  vol.  2,  p.  7), 


(142) 


1"  2,J^(2)  = J 5-Mxcos*+1n* 


0 43) 


where  Jp(z)  is  a Bessel  function  of  the  first  kind.  The  result  is  the 
following  integral  equation  for  the  radial  mode  function. 


• n + 1 
1 Y 


(144) 


where  ^ P]  » 5]  = ^2'  ^1  * 

To  complete  this  step  in  the  procedure,  one  makes  the  sub- 
stitutions u^(C|)  = ^n(CT)/  and  Yp  = yi"^  ^ \ to  obtain  the  foV 
lowing  integral  equation. 
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v„(5i)K„(S,|5i)  <15) 


(145) 


'■n<5l> 


'"i 


where  the  kernel,  K^(C-||c-j).  which  is  given  by 


K„(5,U))  " ■JS,S)  e 


5)4) 


(146) 


is  symmetric  with  respect  to  and  c-J . 

The  next  step  in  the  procedure  is  to  obtain  the  operator  M, 
using  the  manipulations  described  in  the  section  entitled  "The  Chosen 
Procedure."  These  manipulations  involve  the  following  pair  of  eguations 


dCi^ 


d2K^(5lUj) 


dcf 


d?i‘ 


(147) 


/*a 


dV(Ci') 

■'  — ” ' ■ dc-l 
2 • 


d^K  (?^Up 

J — ^ — ! — !-  d£ ' + 


1 j-  I 2 


dC]  + R v^(c])  (148) 


where 


dv  (cn  dK  (qUi) 

Rv„(5))  »k„(5,l5))-^-  V5))— 


(1^9) 


Equation  (147)  was  obtained  by  operating  on  Eq.  (145)  with 

while  Eq.  (148)  was  obtained  by  integrating  the  term  K^(5iU]) 
— g^T2 — by  parts  twice. 


Continuing  to  follow  the  section  entitled  "The  Chosen 
Procedure,"  one  subtracts  Eq.  (148)  from  Eq.  (147),  computes  the 
derivatives  of  K^(5i[5-j)>  and  uses  Eq.  (145)  to  obtain  the  result 
shown  below. 


1 ,1 


H. 


- -Y„Rv„(5j)  K„(5,l5j) 

0 
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Referring  to  Eq.  (124),  the  following  conditions  may  be  shown  to  apply 

1.  f(g.5)  = (g^  - "I)  ^ (1  ■ 

2.  h(g,C^Cp  = 0,  and 

3.  Unless  the  term  in  braces  in  Eq.  (149)  vanishes  at  C-j  = 0, 
Rv^(Cp  will  depend  on  v^(c^)  and  its  derivative  at  the 
origin  as  well  as  the  mirror  edge. 

Thus,  to  meet  all  the  conditions  described  on  page  58,  the  v^(s) 


must  be  selected  so  that  at  = 0, 


(151) 


—jq  3c]  ^ 


For  the  actual  demonstration  that  Eq.  (151)  holds  for  the  selected 

forms  of  the  v„(Ci)9  the  reader  is  referred  to  a more  detailed  version 
n r 

of  this  discussion  in  appendix  D.  For  this  presentation,  it  is  simply 
assumed  that  Eq.  (151)  holds. 

Therefore,  subject  to  this  assumption  and  in  accordance  with 
Eq.  (125),  the  operator  M is  chosen  to  be 

M = (g"  - 1)  ^ (J--  n^)  (152) 


Thus,  the  eigenvalue  problem  of  Eq.  (110)  has  been  reduced  to  solving 
the  following  differential  equation  and  applyino  the  restrictions 
presented  on  paoe  62. 


d\(c) 

~W~ 


(q2  . 1)  ^2  + S2 


v,(c)  = 0 


(153) 


For  stable  resonators,  the  most  convenient  form  of  this  equa- 


tion is 


d\(z) 

dz2 


1 + (2VJL1I  + 

‘ Z Z2 


w^(z)  = 0 


(154) 


This  equation,  which  was  obtained  usina  the  substitutions 
w^(c)  = v^(c)  yr.  z = ot'c^.  s2  = - (2v  + 1)a',  and  a'  = \/l  - g^,  is 
Whittaker's  differential  equation  (Whittaker  and  Watson,  ref.  39, 
p.  337).  It  has  solutions  of  the  form. 
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(155) 


w (z)  = A W (z)  + B W _(-z) 
K.P  ic,p'  -<,p 


i 

t 


where  the  p(z)  are  Whittaker  functions  of  the  second  kind, 

K = - ^ (2v  +1),  and  P “ With  z given  above,  these  functions 
are  characterized  by  the  following: 

1.  Singularities  at  the  origin,  except  for  the  special 

case  where  k = - m - ")  with  m = 0,1 ,2, , and 

2.  Exponential  decay  (for  W „(^))  exponential  growth 
(for  W_^  p("Z))  with  increasing  values  of  the  radial 
coordinate  p. 

Thus,  to  make  the  solutions  given  by  Eq.  (155)  consistent  with 
the  known  behavior  for  the  modes  of  stable  resonators  (i.e.,  the  fields 
do  not  increase  exponentially  with  increasing  mirror  radius),  one  sets 
B = 0 and  k = - m - ( ^ ) with  m = 0,1,2, — . Using  these  conditions 

and  well  known  relationships  involving  Whittaker  functions  (see 
appendix  H),  one  can  show  that  the  expansion  functions  for  stable 
resonators  have  the  form 

('56) 

where  (a'c^)  is  a generalized  Laguerre  Polynomial  (Erdelyi,  ref.  41, 
vol . 1 , p.  268) . 

For  unstable  resonators,  it  is  advantageous  not  only  to  con- 
vert Eq.  (153)  back  to  an  equation  for  Up(c).  but  also  to  make  the  sub- 
stitutions z = nc,T  = - (s/2fl)^,  and  ^ v/o^  - l'  to  obtain  the 
following  equation. 


+ 4z2  - i!-  - 4t  u (z)  = 0 (157) 

-,2 


This  equation  has  solutions  of  the  form 


= r ^iT,n  ^ SW.T,n 
2 2 


(158) 


where  M (z)  and  W (z)  are  Whittaker  functions  of  the  first  and  second 
kind,  respectively.  As  indicated  in  the  discussion  for  stable  resona- 
tors, the  function  W (z)  has  a singularity  at  the  orinin  unless 

1C  yp 

ix  = - (m  + g— ^) with  m = 0,1,2,  — . However,  as  solutions  with 

these  values  of  t do  not  reduce  to  the  functions  shown  in  Eq.  (129),  they 


will  not  be  considered  further.  Then,  in  line  with  the  fact  that 


ix  - (m  + 


n + 1 


■) , the  coefficient  B in  Eq.  (158)  must  be  set  equal  to 


zero.  The  result  is 


2 


(159) 


In  appendix  D,  it  is  shown  that. 


(160) 


where  J^(s5)  is  a Bessel  function  of  the  first  kind  and  C is  a constant. 
Thus,  if  s is  chosen  so  that 


I 


then  the  functions  ^(z)  defined  in  Eq.  (160)  will  reduce  to  the 
functions  specified  in  Eq.  (129).  Accordingly,  the  expansion  functions 
for  unstable  resonators  with  circular  mirrors  are  aiven  by  Eq.  (159) 


nm 

Q 

M +■  (1  + i)b 


(162) 


where  6 = 0.824,  M = y 8ttN  ',  and  N is  the  resonator  Fresnel  number. 

For  resonators  containing  mirrors  with  centered  coupling 

apertures,  there  are  n£  difficulties  involving  singularities  at  the 

origin.  As  a result,  the  functions  W.  ^(iz^),  with  t given  by 

2 

Eq.  (162),  must  be  included  in  the  expansion  set. 

Finally,  to  determine  the  expansion  functions  to  be  used  for 
resonators  with  toroidal  mirrors,  the  reader  is  referred  to  appendix  D. 


Resonators  with  Spherical  Mirrors  of  Rectangular  Projection.  If  the 
same  procedure  is  applied  to  resonators  with  mirrors  of  rectangular 
projection,  the  followino  expansion  functions  are  obtained, 


^ *^"7’  7’  ^ functions 


u(5)  = 


e 5$  (-^,  ^^|^)for  odd  functions  (163) 


where  d^  = 2 v/g^  - 1,  c 


X , and  4i(a,Y,x)  is  a confluent  hyper- 


geometric function  of  the  first  kind  (see  aopendix  H).  For  stable 
resonators,  the  parameter  v must  be  zero  or  a positive  integer.  For 
unstable  resonators. 
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(164) 


mm:  .-.li  iiiii' 


= . i . i(l)2 

2 


with  s = 


mir 


^ ^ . For  even  modes,  m = 1,3,5,  — ; and  for 

odd  modes,  m = 2,4,6, — . 

To  support  the  basic  theory  and  method  of  solution  presented 
in  this  paper,  these  expansion  functions  have  been  used  in  conjunction 
with  the  Rayleigh-Ritz  procedure  to  calculate  the  modes  of  several 
laser  resonators.  The  results  of  these  calculations  are  presented 
and  discussed  in  the  following  chapter. 


This  chapter  has  two  important  objectives.  The  first  is  to 
present  the  results  of  the  numerical  computations  performed  by  applying 
the  basic  theory  and  method  of  solution  discussed  in  the  previous 
chapters.  The  second  is  to  support  the  validity  of  that  theory  and 
method  by  comparing  the  results  with  existing  published  work. 

To  accomplish  these  objectives,  the  chapter  is  organized  in  the 
following  manner.  First,  the  computational  procedure  used  to  obtain 
these  results  is  outlined  and  discussed.  Then,  the  results  obtained  for 
circular  resonators  are  presented  in  some  detail,  followed  by  a summary 
of  the  results  for  rectangular  mirror  resonators,  which  are  detailed  in 
appendix  F.  This  chapter  ends  with  an  evaluation  of  these  calculations 
and  the  associated  procedure.  However,  before  discussing  the  computa- 
tional procedure,  the  manner  in  which  the  term  "mode"  is  used  throughout 
this  chapter  is  briefly  discussed. 

The  reader  may  recall  from  chapter  III  that  laser  resonator 
modes  were  divided  into  three  classes.  The  modes  in  Classes  II  and  III 
were  considered  beyond  the  scope  of  this  paper  and  were  eliminated  from 
further  consideration.  In  addition,  the  following  procedure  was  estab- 
lished for  determining  the  Class  I modes  of  a laser  resonator. 

1.  Assume  the  self-induction  terms  in  the  coupled  integral 
equations  (Eqs.  (53)  and  (54))  are  negligible. 

2.  Solve  the  resulting  integral  equations  (Eqs.  (65)  and  (66)). 

3.  Show  that  the  solutions  are  consistent  with  the  assumption 
that  the  self-induction  terms  are  negligible. 

Since  the  calculations  in  this  chapter  were  performed  only  to 
check  the  basic  theory  and  method  of  solution  presented  earlier,  the 
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third  step  of  this  procedure  has  not  been  accomplished.  Despite  that 
fact,  to  remain  consistent  with  commonly  used  terminology,  the  term 
"mode"  is  used  to  designate  all  solutions  to  Eqs.  (65)  and  (66)  (or 
specialized  forms  of  these  equations)  regardless  of  whether  these 
solutions  produce  self-induction  terms  that  are  negligible. 

Computational  Procedure 

This  section  is  divided  into  three  parts.  The  first  discusses 
an  ideal  procedure  for  using  the  Rayleigh-Ritz  technique  to  determine 
laser  resonator  modes.  The  second  oart  covers  the  computational  diffi- 
culties encountered  as  a result  of  the  limitations  of  the  computer 
programs  used  by  the  author  for  these  calculations.  This  ideal  pro- 
cedure and  these  limitations  combine  to  yield  the  actual  procedure, 
which  is  discussed  in  part  three. 

Ideal  Procedure.  The  steps  of  an  ideal  procedure  for  applying 
the  Rayleigh-Ritz  technique  to  laser  resonator  problems  are  listed 
below. 

1.  Select  a number  of  expansion  functions  for  the  first 
stage  of  the  calculation. 

2.  Using  that  set  of  expansion  functions,  compute  the  eigen- 
values and  corresponding  mode  distributions  (magnitude  and 
phase) . 

3.  Increase  the  number  of  expansion  functions  and  repeat  steps 
1 and  2. 

4.  Compare  the  results  obtained  using  the  two  sets  of  functions 
to  determine  if  the  desired  accuracy  has  been  achieved. 

5.  If  the  desired  accuracy  has  been  achieved,  the  procedure  is 
terminated.  If  not,  steps  3 and  4 are  repeated. 

The  convergence  criterion  used  in  step  4 depends  on  the  data 
desired  from  the  calculation.  If  only  eigenvalues  are  required,  the 
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details  of  the  mode  distributions  need  not  be  considered  since  the 
eigenvalues  are  relatively  insensitive  to  small  perturbations  in  the 
mode  distributions.  Thus,  convergence  with  respect  to  eigenvalue  data 
may  be  achieved  with  a few  expansion  functions.  On  the  other  hand,  many 
expansion  functions  may  be  required  (at  least  10)  to  determine  the 
details  of  the  modes  of  large  magnification  resonators. 

Limitations  of  Existing  Programs.  Experience  obtained  from 
many  applications  of  this  procedure  indicates  that,  for  resonators 
departing  significantly  from  g = 1.0  (g>1.25),  accurate  calculations 
with  at  least  lo  expansion  functions  are  required  to  determine 

1.  The  eigenvalues  of  several  modes,  or 

2.  Any  mode  distributions. 

As  shown  in  chapter  V,  the  expansion  functions  selected  for 
all  resonators  considered  in  this  paper  are  directly  related  to  the 
confluent  hypergeometric  functions  of  the  first  ($(a,Y,z))  and  second 
i'v{a,r,z))  kind.  In  this  analysis,  the  parameter  a has  a strong 
dependence  on  the  number  of  functions  used  (see  Eqs.  (162)  and  (163) 
in  chapter  V).  That  is  to  say,  as  the  number  of  expansion  functions 
increases,  the  magnitude  of  a for  the  last  function(s)  becomes  large. 
This  large  magnitude  makes  it  difficult  to  compute  these  functions 
accurately  for  certain  ranges  of  the  spatial  variable  z.  These  ranges 
correspond  to  intermediate  values  of  z,  where  z is  too  small  to  apply 
asymptotic  series  expansions  and  where  round-o'ff  error  difficulties 
are  encountered  with  the  usual  Taylor  series  expansions. 

As  a result  of  these  limitations,  which  were  inherent  in  the 
computer  routines  readily  available  to  the  author,  it  was  not  always 
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possible  to  perform  accurate  computations  with  the  10  expansion 
functions  needed.  Due  to  these  inaccuracies,  the  ideal  procedure 
was  modified  as  discussed  in  the  followinq  paragraphs. 

Actual  Procedure.  The  procedure  actually  employed  in  this 
paper  depended  on  the  particular  resonator  being  studied.  For  reso- 
nators with  small  equivalent  Fresnel  number,  N^,  the  ideal  procedure 
was  used.  Typically,  however,  only  the  first  two  steps  of  that  pro- 
cedure were  performed  since  the  results  obtained  for  the  lowest  loss 
modes  at  that  stage  of  the  calculation  were  in  close  agreement  with 
those  obtained  by  other  authors.  For  these  cases,  only  the  results 
for  the  first  three  or  four  modes  are  reported.  The  eigenvalues  for 
higher  loss  modes,  which  are  not  reoorted,  are  considered  as  not  having 
converged. 

For  resonators  corresponding  to  the  regions  of  difficulty 
mentioned  on  the  previous  page,  this  procedure  was  modified  further. 

For  these  cases,  a series  of  eigenvalue  computations  was  performed 
using  a relatively  small  number  (usually  six  or  seven)  of  expansion 
functions  for  each  case.  For  these  resonators  only,  the  series 
reported  corresponds  to  the  one  yielding  the  best  agreement  in  the 
eigenvalue  of  the  lowest  order  mode.  For  example,  the  data  reported 
later  for  g = 2.6  fall  into  this  category. 

The  rationale  behind  selectina  the  series  yielding  the  best 
agreement  points  up  the  main  computational  difficulty  encountered  in 
this  paper.  To  understand  the  difficulty,  it  is  helpful  to  recall 
that,  with  any  expansion  procedure,  one  expects  the  accuracy  of 
calculations  to  increase  as  the  number  of  expansion  functions  is 
increased.  However,  due  to  the  difficulties  with  the  higher  order 
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expansion  functions  (those  with  large  |a|  discussed  on  the  previous  page), 
that  accuracy  could  actually  decrease  as  the  number  of  functions  is 
increased.  Thus,  with  the  tools  available,  there  was  an  optimum  number 
of  expansion  functions  for  each  resonator  in  this  category.  Tt  was 
assumed  that  the  optimum  number  corresponded  to  the  one  yielding  the 
best  agreement  (with  existing  published  work)  for  the  lowest  loss  mode. 

The  results  obtained  using  this  modified  procedure  are  presented 
in  the  following  two  sections. 

Presentation  of  Circular  Mirror  Results 

Two  types  of  results  are  presented  in  this  section.  First, 
data  related  to  the  eigenvalues  in  the  integral  equation  are  presented 
for  a wide  range  of  symmetric  resonators  with  circular  mirrors.  Except 
when  mentioned,  these  data  are  presented  in  the  same  format  used  by 
other  authors.  Second,  plots  are  presented  of  the  relative  magnitude 
and  phase  of  the  current  induced  on  one  mirror  of  a plane  parallel 
resonator  with  M = 10.  For  brevity,  only  plots  for  the  first  two 
azimuthally  symmetric  modes  are  included  in  this  section.  For  other 
mode  distributions  involving  this  resonator,  see  the  first  section  of 
appendix  F. 

Eigenvalue  Data.  As  indicated  earlier,  the  eigenvalue  data 
presented  in  this  section  cover  a broad  range  of  laser  resonators.  The 
results  for  the  first  case,  a olane  parallel  resonator  with  N = 10,  are 
summarized  in  tables  II  and  III.  The  data  used  for  comparison  were 
taken  from  Fox  and  Li  (ref.  25,  p.  465,  table  I),  and  the  indicated 
percentage  power  loss  corresponds  to  a single  pass  through  the 
resonator. 
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Table  II 


Percentage  Power  Loss  for  Circular  Mirror 
Resonator  with  N =10  and  g = 1.0 


DOUGHTY 

FOX  & LI 

VAINSHTEIN 

DESIGNATION 

1 

0.848 

0.84 

0.82 

TEMoo 

w 

\ ! 
r 

' 

4.447 

4.45 

4.31 

TEMio 

11.198 

10.6 

10.6 

TEM20 

1 

19.579 

18.7 

19.7 

TEM30 

2.019 

2.02 

2.08 

TEMoi 

6.857 

6.66 

6.96 

TEMii 

14.452 

14.4 

14.6 

TEM21 

: 

24.737 

23.3 

25.1 

TEM31 

r 

Table  III 

i 

Phase 

Shift  for  Circular  Mirror  Resonator  with 

N = 10  and 

g = 1.0  (in  degrees) 

DOUGHTY 

FOX  & LI 

VAINSHTEIN 

DESIGNATION 

- i 

2.36 

2.36 

2.75 

TEMqo 

[• 

12.5 

12.4 

11.9 

TEMio 

30.7 

30.7 

29.2 

TEM20 

57.2 

57.0 

54.2 

TEM30 

6.025 

6.03 

5.72 

TEMoi 

1 20.15 

20.1 

19.2 

TEMn 

42.35 

42.2 

40.3 

TEM21 

73.0 

73.0 

69.2 

TEM31 

i 
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Comparison  of  the  results  in  these  tables  reveals  that  the 
maximum  difference  in  the  percentage  power  loss  is  less  than  6 percent; 
and  in  most  cases,  the  results  obtained  in  this  paper  fall  between  those 
obtained  by  Fox  and  Li  and  L.  Vainshtein.  Also,  the  values  computed 
for  the  phase  shift  are  in  almost  perfect  agreement  with  those  computed 
by  Fox  and  Li. 

The  next  case  considered  was  a circular  mirror  resonator  with 
g = 1.1.  The  results  for  the  first  four  azimuthal ly  symmetric  modes  for 


Ng  = 1,2,  and  3 are  presented  in  tables  IV,  V,  and  VI.  For  those  cases. 


the  eigenvalue  phase,  arg  y,  corresponds  to  the  round  trip  eigenvalue 
into  which  has  been  absorbed  the  terms  e i . A brief  examina- 

tion of  these  results  reveals  that  the  maximum  difference  in  the 
magnitude  of  the  one-way  eigenvalue  is  less  than  3.5  percent,  and  the 
maximum  phase  difference  is  less  than  11  degrees.  Comparison  data 
were  taken  from  Siegman  and  Miller  (ref.  28,  p.  2734,  fig.  5). 


Table  IV 

Eigenvalue  Data  for  Circular  Mirror 


Resonator  with  n = 0,  g = 1.1 , and  = 1 


DOUGHTY 


SIEGMAN  & MILLER 


IyI 

ARG  Y 

>Th 

ARG 

0.816 

143.8° 

0.82 

140. 

0.778 

-134.2° 

0.78 

-140. 

0.544 

+ 2.9° 

0.54 

+ 10. 

0.308 

-136.9° 

0.30 

-140. 

n 
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Table  V 

Eigenvalue  Data 
with  n = 

for  Circular  Mirror  Resonator 
0,  g = 1.1,  and  N = 2 

DOUGHTY 


SIEGMAM  & MILLER 


ARG 


1.76 

1 

1.78 

-1 

1.67 

+ 

1.60 

- 

Table  VI 

Eigenvalue  Data  for  Circular  Mirror  Resonator 
with  n = 0,  q = 1.1,  and  = 3 


DOUGHTY 


SIEGMAN  & MILLER 


>Th 


■ 

.68 

+ 

.62 

The  third  case  considered  was  a syiroietric,  circular  mirror  resonator 
with  g = 1.25.  The  magnitudes  of  the  eigenvalues  (one-way)  for  the  first 
four  modes  for  N = 4 and  N = 8 are  presented  in  tables  VII  and  VIII. 
Examination  of  these  results  reveals  two  things.  First,  agreement 
between  the  results  for  the  first  two  modes  with  N = 4 and  for  the 
first  mode  with  N = 8 is  excellent.  Second,  agreement  for  other 
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modes  of  these  resonators  is  generally  poor.  The  reason  for  this 
disagreement  is  discussed  later  in  this  chapter.  Comparison  results 
for  this  case  were  also  taken  from  Siegman  and  Miller  (ref.  28,  p.  2733, 
fig.  3). 


Table  VII 


Eigenvalue  Magnitudes  for  Circular  Mirror 
Resonator  with  n = 0,  g = 1 .25,  N = 4 


DOUGHTY  >/TyT 

SIEGMAN  & MILLER 

0.596 

0.59 

0.560 

0.56 

0.559 

0.46 

0.429 

0.42 

Table  VIII 

Eigenvalue  Magnitudes  for 
Resonator  with  n = 0,  g 

Circular  Mirror 
= 1.25,  N = 8 

DOUGHTY  IyI 

SIEGMAN  & MILLER 

0.631 

0.63 

0.468 

0.58 

0.425 

0.46 

0.420 

0.45 

i 
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Finally,  the  results  for  a range  of  equivalent  Fresnel  numbers 
for  a symmetric,  circular  mirror  resonator  with  g = 2.6  are  summarized 
in  table  IX.  As  before,  the  eigenvalue  phase,  arg  y.  corresponds  to 
the  round  trip  eigenvalue  for  which  the  spatially  independent  terms  have 
been  absorbed.  In  addition  to  reporting  the  basic  results,  the  number 
of  expansion  functions  used  to  obtain  these  data  is  shown  in  the  far 
right  column  of  this  table. 


Table  IX 

Eigenvalue  Data  for  Circular  Mirror  Resonator 
with  n * 0,  q = 2.6  for  Several  Values  of 


DOUGHTY 

VlT  _AF 


SIEGMAN  & MILLER 
>Tm~  arg  y 


Number 


0.235  - 4.4“ 

0.202  -179.° 


0.255 

0.225 

0.125 

0.333 

0.212 

0.173 

0.278 

0.149 

0.284 

0.129 


108.6° 
- 94.8° 
127.4° 


■135.6“ 


•115.1° 


0.285 


+100.' 

-180.' 


+100.' 


+100.' 


-no.' 


I 


AMPLITUDE 


As  indicated  in  the  discussion  of  the  computational  procedure, 
the  results  obtained  using  only  six  or  seven  expansion  functions  may  be 
marginal,  especially  for  the  higher  loss  modes.  That  expectation  is 
certainly  borne  out  by  the  results  summarized  in  this  table.  However, 
despite  that,  the  maximum  difference  for  the  lowest  loss  mode  is  less 
than  16  percent  (Siegman  and  Miller,  ref.  28,  p.  2732,  fig.  2),  and 
the  difference  is  typically  less  than  10  percent.  As  can  be  seen,  the 
magnitude  of  the  difference  for  higher  loss  modes  is  rather  sporadic, 
but  generally  high.  These  results  are  discussed  further  following  the 
presentation  of  the  circular  mirror  mode  plots  and  the  summary  of  the 
data  for  rectangular  mirror  resonators. 

Mode  Plots.  This  portion  of  the  chapter  contains  plots  of  the 


o 


Figure  8.  Relative  Magnitude  and  Phase  Distributions 
for  the  Lowest  Loss,  n = 0 Mode 


Figure  9.  Relative  Magnitude  and- Phase  Distributions 
for  the  Second  Lowest  Loss,  n = 0 Mode 


relative  magnitude  and  phase  for  the  twa  lowest  loss,  azimuthally 
symmetric  modes  for  a symmetric  resonator  with  N = 10.  The  results 
for  the  lowest  loss  mode,  denoted  by  +,  are  presented  in  figure  8, 
while  the  results  for  the  next  mode,  also  denoted  by  -t-,  are  given 
in  figure  9.  The  distributions  to  which  these  data  are  compared, 
denoted  by  *,  were  taken  from  Fox  and  Li  (ref.  25,  p.  464,  fig.  2). 

As  can  be  seen,  there  is  excellent  agreement  between  these 
results  and  those  obtained  by  Fox  and  Li.  Similar  agreement  was 
obtained  for  several  higher  loss  modes  of  this  resonator  (including 
several  modes  which  are  not  azimuthally  symmetric).  These  results 
are  presented  in  appendix  F. 
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Summary  of  Rectangular  Mirror  Results 

As  was  the  case  for  circular  mirrors,  the  computations 
performed  for  rectangular  mirror  resonators  fall  into  two  categories. 
These  two  categories  are  eigenvalue  computations  and  mode  plots  con- 
sisting of  the  relative  distribution  of  the  intensity  (magnitude 
squared)  and  phase  for  the  current  induced  on  one  of  the  resonator 
mirrors. 

The  eigenvalue  results  correspond  to  one  of  the  cases  listed 

below. 

1.  Plane  parallel  resonators  with  N = 10  or  N = 8/tt. 

2.  Power  loss  computations  for  odd  symmetric  modes  with 
g = 1.2. 

3.  Eigenvalue  magnitude  computations  for  even  symmetric 
modes  with  g = 1.8. 

All  of  these  computations  were  performed  for  a single  rectangular 
component  (strip  case)  of  a symmetric  resonator  so  that  the  eigen- 
value corresponded  to  a single  pass  (one-way)  through  the  resonator. 
These  results,  which  are  presented  along  with  comparative  data  in 
appendix  F,  are  summarized  below. 

As  expected,  the  plane  parallel  resonator  results  for  the  two 
Fresnel  numbers  shown  were  excellent.  For  these  two  cases,  the 
relative  difference  in  the  percentage  power  loss  for  any  of  the  first 
three  even  symmetric  modes  was  less  than  2.5  percent.  Similarly,  the 
results  for  the  g = 1.2  unstable  resonator  for  equivalent  Fresnel 
numbers  of  1,  1.5,  and  2 were  quite  good.  In  each  of  these  cases,  the 
difference  in  the  percentage  power  loss  for  the  first  two  even  symmet- 
ric modes  was  less  than  2 percent.  In  any  case,  the  maximum  difference 
was  7 percent.  Finally,  computations  were  performed  for  the  magnitudes 
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of  the  eigenvalues  for  at  least  the  first  three  modes  of  the  g = 1.8 
resonator  mentioned  above.  For  that  series,  nine  values  for  the 
equivalent  Fresnel  number  were  selected  so  that  2.5  - - 20/tt.  For 

this  series,  the  difference  was  typically  4 percent  or  less,  although 
in  one  isolated  case,  it  did  reach  7 percent. 

The  relative  intensity  and  phase  distributions  determined  for 
rectangular  mirror  resonators  also  corresponded  to  symmetric  resonators 
with  g = 1.8.  These  computations  were  limited  to  determining  the 
distributions  for  the  first  two  even  symmetric  modes  for  values  of  Ng 
such  that  2irNg  = 18,  34,  36,  and  40.  With  the  oossible  exception  of 
the  phase  distribution  of  the  second  mode  for  2irNg  = 34,  these  distribu- 
tions are  certainly  not  in  strong  disagreement  with  those  obtained  by 
Sanderson  and  Streifer  (ref.  42,  figs.  13-21).  In  fact,  in  most  cases, 
these  distributions  can  be  said  to  be  in  qualitative  agreement  with  those 
obtained  by  these  authors.  By  qualitative  agreffnent,  I mean  that 
although  the  relative  values  for  the  peaks  and  troughs  may  differ  some- 
what, the  basic  nature  and  trends  for  the  distributions  are  the  same. 

The  reason  for  any  deviations  in  these  results  is  covered  in 
the  following  section. 

Discussion 

The  purpose  of  this  section  is  to  evaluate  the  computational 
procedure  employed  for  this  analysis,  as  well  as  the  results  obtained 
using  that  procedure.  The  procedure  itself  will  be  considered  first. 

From  earlier  discussions,  it  is  apparent  that  the  existing 
procedure  is  inadequate  for  determining  the  modal  characteristics  of 
a wide  range  of  laser  resonators.  However,  the  shortcomings  associated 
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with  that  procedure  result  from  the  limited  range  of  applicability 
of  the  author's  computer  routines  and  are  not  due  to  any  fundamental 
deficiency  inherent  in  the  procedure  itself.  In  addition,  the  purpose 
of  these  computations  is  somewhat  limited  in  that  they  were  performed 
only  to  support  the  basic  theory  and  method  of  solution  presented  in 
the  earlier  chapters.  When  evaluated  in  terms  of  this  limited  purpose, 
the  computational  procedure  is  certainly  adequate. 

To  aid  in  evaluating  the  results  obtained  using  this  procedure, 
the  results  for  which  excellent  agreement  (difference  generally  5 
percent  or  less)  was  achieved  are  listed  below. 

1.  Mode  distributions  for  plane  parallel  resonators  with 
circular  mirrors. 

2.  Eigenvalue  data  for  plane  parallel  resonators  with  both 
circular  and  rectangular  mirrors. 

3.  Eigenvalue  magnitudes  and  phases  for  circular  mirror 
resonators  with  g = 1.1. 

4.  Lowest  loss  eigenvalues  for  circular  mirror  resonators 
with  g = 1.2. 

5.  Eigenvalues  for  odd  modes  for  rectangular  mirror 
resonators  with  g = 1.2. 

6.  Eigenvalues  for  even  modes  for  a wide  range  of  rectangular 
mirror  resonators  with  g = 1.8  (one  point  had  a difference 
of  7 percent). 

Listed  below  are  the  results  which  are  considered  as  fair  or 
marginal . 

1.  Eigenvalues  for  higher  order  modes  for  a circular  mirror 
resonator  with  g = 1 .2  and  N = 8. 

2.  Eigenvalues  for  several  modes  for  circular  mirror 
resonators  with  g = 2.6. 

3.  Mode  distributions  for  rectangular  resonators  with  g = 1.8. 


All  of  the  cases  falling  into  this  second  list  correspond  to  instances 
where  computational  difficulties  were  expected.  Difficulties  were 
expected  for  cases  1 and  2 as  they  correspond  to  eigenvalue  computa- 
tions (mostly  for  higher  loss  modes)  for  which  the  equivalent  Fresnel 
number  is  in  the  intermediate  to  high  range.  Difficulties  were 
expected  for  case  3 because  the  actual  mode  distributions  are  espe- 
cially sensitive  to  errors  in  the  values  computed  for  the  expansion 
functions.  However,  it  is  worth  mentioning  that  the  basic  nature  or 
character  of  the  modes  is  correctly  described,  even  though  some 
differences  are  present  in  these  plots. 

From  this  discussion,  one  can  see  that  in  regions  where  no 
computational  difficulties  were  expected,  the  results  obtained  in  this 
analysis  are  in  excellent  agreement  with  those  obtained  by  other 
workers.  In  addition,  in  every  case  where  agreement  was  fair  or  poor, 
that  lack  of  agreement  can  be  correctly  attributed  to  the  computer 
routines  employed  and  not  to  the  basic  theory  or  method  of  solution. 
Thus,  the  numerical  results  discussed  in  this  chapter  and  the  analytic 
results  covered  in  appendix  E strongly  support  the  validity  of  the 
basic  theory  and  method  of  solution  for  symmetric  resonators  of  various 
geometries,  Fresnel  numbers,  and  magnifications.  More  importantly, 
these  same  results  should  inspire  considerable  confidence  in  the 
applicability  of  this  theory  and  method  of  solution  for  other  more 
complex  systems. 


Summary,  Conclusions,  and  Recommendations 


The  basic  theory  needed  to  analyze  the  modes  of  complex  laser 
resonators  containing  homogeneous  media  has  been  presented  in  this 
report.  That  theory,  which  is  firmly  based  on  the  fundamental  equations 
and  principles  of  electromagnetic  theory,  culminates  in  a pair  of 
integral  equations  for  the  electromagnetic  field  within  the  resonator. 

It  is  this  pair  of  equations  which  must  be  solved  simultaneously  to 
determine  the  modes  of  an  arbitrary  laser  resonator.  Unfortunately, 
these  equations  are  very  complex,  and  solving  them  for  the  general  case 
is  far  beyond  the  scope  of  this  effort. 

As  a result  of  this  complexity,  the  integral  equations  were 
gradually  simplified  by  making  various  assumptions  concerning  the  nature 
of  the  resonators  being  analyzed.  For  example,  it  was  assumed  that  the 
resonators  consist  of  two  perfectly  conducting  mirrors  and  a homogeneous 
medium.  This  assumption  resulted  in  a pair  of  equations  for  the  currents 
induced  on  the  resonator  mirrors.  It  was  further  assumed  that  the 
effect  of  any  gain  or  amplification  (resulting  from  polarization  of  the 
lasing  medium)  is  essentially  constant  between  any  two  points  on  dif- 
ferent mirrors.  In  addition,  the  oscillation  wavelength  was  assumed  to 
be  small  in  comparison  to  all  resonator  dimensions,  and  the  mirror 
diameters  were  assumed  to  be  small  in  comparison  to  the  mirror  separa- 


tion and  radii  of  curvature.  Finally,  it  was  assumed  that  the  self- 


After  making  these  assumptions,  the  resulting  equations  for  the 
modal  currents  were  solved  using  a technique  that  combines  a variational 
principle  with  the  novel  expansion  functions  discussed  in  chapter  V.  By 
comparing  these  solutions  to  existing  published  work,  the  validity  of 
the  basic  theory  and  method  of  solution  presented  here  has  been  verified 
for  a wide  range  of  symmetric  resonators.  However,  in  addition  to 
directly  verifying  the  analysis  for  these  symmetric  resonators,  these 
same  calculations  and  comparisons  strongly  support  the  assertion  that 
the  theory  and  method  of  solution  are  also  applicable  for  more  complex 
systems. 

In  addition  to  providing  the  basis  for  actually  computing  the 
mode  losses  and  distributions,  the  derivation  of  the  basic  equations 
sheds  considerable  light  on  laser  resonators  and  their  physical  char- 
acteristics and  fields.  For  example,  the  derivation  clearly  shows  that 
the  modes  of  open  resonators  are  damped  with  time.  For  many  cases,  such 
as  the  lower  loss  modes  of  paraxial  resonators,  this  damping  has  a 
negligible  effect  on  the  actual  mode  distributions.  However,  in  other 
cases,  its  effect  can  be  significant.  In  fact,  by  properly  including 
the  effect  of  uamping  at  the  beginning  of  the  analysis,  it  has  been 
shown  that  it  is  not  always  possible  to  formulate  the  laser  resonator 
problem  in  terms  of  a linear  eigenvalue  problem.  In  addition,  by 
properly  including  the  temporal  behavior  of  the  modes,  this  derivation 
has  eliminated  the  need  to  define  laser  resonator  modes  in  terms  of  some 
mythical  round  trip  through  the  resonator.  In  fact,  from  this  analysis, 
the  modes  are  seen  to  be  fields  for  which  the  relative  distributions 
do  not  change  at  any  time. 
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This  analysis  also  points  out  the  need  to  distinguish  between 
the  two  types  of  energy  flow  associated  with  resonant  systems.  More- 
over, by  using  continuity  arguments,  it  shows  that  for  energy  flowing 
from  the  resonator  to  the  surroundings,  the  open  portion  of  the  reso- 
nator hull  does  not  enter  into  the  calculations  for  the  fields  within 
the  resonator. 

In  addition  to  the  points  mentioned  in  the  preceding  paragraphs, 
this  analysis  clearly  shows  that  one  cannot  always  analyze  the  modes  of 
laser  resonators  by  using  only  integral  equations  which  express  the 
current  (or  field)  on  one  mirror  entirely  in  terms  of  the  current  (or 
field)  on  the  other  mirror.  In  fact,  since  both  mirrors  normally  affect 
the  current  distributions,  one  should  generally  use  a pair  of  coupled 
equations  which  involve  both  self-induction  and  mutual-induction  terms. 
As  a result,  analyses  which  do  express  the  current  on  one  mirror  en- 
tirely in  terms  of  the  current  on  the  second  mirror  have  two  potentially 
important  deficiencies.  First,  such  analyses  completely  neglect  modes 
for  which  the  mutual-induction  and  self-induction  terms  are  either 
proportional  to  each  other  or  have  comparable  magnitudes.  Second,  some 
of  the  distributions  normally  regarded  as  resonator  modes  may  not  be 
modes  at  all.  This  statement  is  based  on  the  fact  that  some  of  these 
distributions  may  not  be  consistent  with  the  assumption  that  the  self- 
induction  terms  in  the  coupled  integral  equations  are  negligible. 

Another  important  point  in  this  analysis  is  the  fact  that  the 
fields  associated  with  open  resonators  may  be  discontinuous.  These 
discontinuities,  which  may  occur  at  the  mirror  edges  or  at  material 
interfaces  in  segmented  mirrors,  are  treated  simply  by  allowing  for  the 
presence  of  charges  and  currents  at  the  points  of  discontinuity. 
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Finally,  by  comparinq  the  equations  obtained  at  different  stages  of  this 
analysis  to  those  obtained  by  other  authors,  one  can  gain  considerable 
insight  into  the  actual  range  of  applicability  of  other  analyses.  For 
example,  the  integral  equations  obtained  in  this  paper  for  paraxial 
resonators  with  perfectly  conducting  mirrors  are  identical  to  those 
obtained  by  other  authors  under  what  appear  to  be  less  stringent  con- 
ditions. As  this  analysis  is  strongly  based  on  the  principles  of 
electromagnetic  theory,  this  suggests  that  the  equations  obtained  by 
those  authors  are  not  as  general  as  they  are  normally  considered. 

In  addition  to  the  insight  obtained  from  the  equations  and 
derivations  themselves,  one  can  gain  further  understanding  of  the  be- 
havior of  resonator  fields  by  actually  calculating  mode  distributions 
and  losses.  The  variational  method  presented  in  this  paper  is  one 
excellent  technique  for  performing  these  calculations.  Using  this 
technique,  which  is  implemented  by  apolying  the  Rayleigh-Ritz  procedure 
in  conjunction  with  the  expansion  functions  developed  in  chapter  V,  the 
analysis  of  laser  resonator  modes  is  reduced  to  a matrix  problem  that 
can  be  solved  using  standard  numerical  techniques. 

This  reduction  of  the  analysis  to  a matrix  problem,  which  can  be 
solved  using  standard  numerical  techniques,  has  two  important  advantages 
over  the  standard  iteration/orthogonal ization  methods  of  solution. 

First,  the  variational  method  is  usually  considerably  faster  than  the 
normal  iteration  approach,  which  actually  corresponds  to  solving  the 
integral  equation  by  the  method  of  successive  approximations.  In 
addition,  the  variational  approach  yields  information  concerning  several 
modes  at  once,  while  the  iteration  approach  requires  that  a complete 
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iteration/orthogonal ization  sequence  be  performed  for  each  mode. 

Finally,  the  method  of  solution  presented  in  this  paper  has  at  least  one 
important  advantage  over  other  matrix  approaches.  This  advantage,  which 
is  a result  of  the  way  in  which  the  expansion  functions  are  obtained,  is 
that  modal  calculations  can  be  performed  using  a relatively  small  number 
of  expansion  functions. 

As  indicated  earlier,  the  numerical  work  presented  in  this  paper 
has  the  rather  limited  purpose  of  supporting  the  basic  theory  and  method 
of  solution  covered  in  the  previous  chapters.  With  that  in  mind,  the 
numerical  work  was  performed  using  a collection  of  computer  routines 
which  were  easily  accessible  to  this  author.  No  extensive  effort  was 
made  to  write  or  find  routines  which  were  very  aeneral  or  extremely 
efficient.  As  a result,  this  collection  needs  inprovement  in  two 
important  areas.  First  and  foremost,  the  subroutine  used  to  compute  the 
confluent  hypergeometric  functions  should  be  modified  so  that  these 
functions  can  be  accurately  evaluated  for  intermediate  values  of  the 
spatial  variable  z.  This  improvement  is  essential  if  this  method  of 
solution  is  to  find  widespread  application.  Second,  the  integration 
routines  in  this  collection  should  be  modified  in  a manner  which  will 
reduce  the  effect  of  the  rapid  variation  associated  with  the  kernels  of 
these  integral  equations.  Such  a modification  would  allow  accurate 
computation  of  the  matrix  elements  with  a considerable  time  savings. 

This  is  especially  true  for  resonators  with  large  equivalent  Fresnel 
numbers. 

Another  aspect  requiring  improvement  is  the  analysis  of  the 
conditions  under  which  the  self-induction  terms  in  the  integral  equations 
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are  negligible.  The  estimates  given  in  this  paper  will  probably  suffice 
for  the  lower  loss  modes  of  paraxial  resonators;  however,  they  may  begin 
to  break  down  as  k"  increases  or  as  the  mirror  radii  become  large.  As  a 
result,  more  accurate  and  more  general  estimates  of  these  conditions  are 
definitely  needed. 

Following  these  improvements  in  the  analysis  and  computer 
routines,  this  analysis  should  be  extended  and  applied  to  study  the 
modes  of  a variety  of  laser  resonators.  For  instance,  this  analysis 
should  definitely  be  extended  to  determine  the  characteristics  of  Class 
II  and  Class  III  modes  of  conventional  stable  and  unstable  laser  reso- 
nators. It  should  also  be  used  to  study  toroidal  and  hole-coupled 
resonators  to  verify  that  the  expansion  functions  obtained  for  these 
cases  have  the  desired  characteristics.  When  the  applicability  of  these 
functions  has  been  verified,  and  when  the  above  improvements  have  been 
made,  the  basic  theory,  method  of  solution,  and  computational  procedure 
presented  in  this  report  should  be  extremely  valuable  tools  in  the  modal 
analysis  of  paraxial  resonators  with  perfectly  conducting  mirrors. 
However,  even  with  these  improvements,  this  work  will  not  be  the  last 
word  in  resonator  theory.  On  the  other  hand,  the  integral  equations 
derived  in  this  paper  can  and  should  form  the  basis  for  studies  of  more 
complex  resonators  for  which  one  or  more  of  the  limiting  assumptions  do 
not  apply. 
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APPENDIX  A 


Green's  Function  Derivation 


The  objective  of  this  appendix  is  to  derive  an  outgoing  wave 
Green's  function  to  be  used  in  calculating  laser  resonator  modes.  This 
Green's  function  is  a solution  of  Eq.  (20),  which  is  repeated  below 

7'2  S + k.2  r,  = - ! sir  - r')  (Al) 

J 

where  k.  = k'.  + k'J  with  k'.  > 0,  T = a a + a a + a a , 6(r  - r')  is 
JJJ  J xxyyzz 

the  Dirac  delta  function,  and  the  Laplacian  is  understood  to  operate 
on  the  rectangular  components  of  Substituting  the  form  5 = !(})  into 
Eq.  (Al),  one  obtains  the  following  equation  for  <t>. 

V'2  ,|,  + k.2  4,  = - 6(r  - r')  (A2) 

vJ 

To  solve  Eq.  (A2),  it  is  helpful  to  express  both  <i>  and  the 
delta  function  as  Fourier  integrals  of  the  form 


/ 


dF  e"'' 


where  R = r - r ‘ , and  the  integrals  are  evaluated  over  all  F.  Sub- 
stituting these  expressions  into  Eq.  (A2),  evaluatino  v'2 
solving  for  and  substituting  the  result  into  Eq.  (A3)  leads  to  the 
following  expression  for  <t>. 


^ = 


g-iF.1T 
|F|2  - kj2 


(A5) 


Integrating  over  angles  and  denoting  |F|  by  k and  |1T1  by  R, 
this  expression  is  simplified  to  yield  Eg.  (A6). 


^ " (2ir)2iR 


/ 


k dk 

(k  + k.)  (k  - k.) 


(A6) 


By  writing  the  right  side  as  the  sum  of  two  integrals  and  letting 

1 kR 

k ^ -k  in  the  integral  involving  e , Eq.  (A6)  may  be  rewritten  in 
the  form  shown  below. 


1 f k dk 

(2tt)2iR  j (k  + kj)  ■(■k“ncj) 


(A7) 


To  obtain  the  outgoing  wave  solution  for  kj'  < 0,  one  evaluates 
the  integral  in  Eq.  (A7)  using  contour  integration  and  the  theory  of 
residues.  The  closed  contour  in  the  complex  k-plane,  which  is  traversed 
in  a counterclockwise  direction,  includes  the  real  axis  and  a semi- 
circle of  infinite  radius  (centered  at  the  origin)  in  the  upper  half- 
plane as  shown  in  figure  10.  Since  the  integral  over  the  semicircle 
vanishes,  (()  is  given  by  Eq.  (A8), 


with  the  contour  and  pole  diagram  shov/n  in  Figure  10. 


(A8) 
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Fiqure  10.  Closed  Contour  in  the 
Complex  k-Plane 

To  avoid  systematically  excluding  outgoing  wave  solutions  for 
k'l  ^0,  one  must  analytically  continue  this  solution  into  the  upper 

J 

half  of  the  complex  k^.-plane.  To  do  this,  one  simply  distorts  the 
contour  to  insure  that  the  path  of  inteoration  includes  the  outgoing 
wave  pole  as  this  pole  crosses  the  real  axis.  That  is,  the  contour 
is  distorted  so  that  it  does  not  cross  any  poles  as  -k.  crosses  the 

V 

real  axis  so  that  -kj  takes  on  negative  values.  This  contour,  which 
analytically  continues  <p  into  the  upper  half  of  the  k. -plane,  is 

J 

shown  in  figure  11 . 
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Figure  11.  Contour  Defining  the  Analytic 
Continuation  of 


i 
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Thus  for  all  values  of  kl|,  the  only  pole  enclosed  by  the  con- 
tours  is  at  k = -k..  Accordingly,  the  contour  integrals  can  be 

sJ 

written  as 


/ 


k dk  _ o . 

- k.J-  - 


(A9) 


where  L',  the  residue  of  the  integrand  at  k = -k.,  is  given  by  Eq.  (AID) 

J 


L'  - 


lim  ('<  * '‘il  ® 


+ikR 


(k  + k.)  (k  - kj) 


1 -ik<  R 
je  3 


(AlO) 


102 


1 


Thus,  the  function  ()i  has  the  form 


-ik.R 


4ttR 


(An) 


for  kl  > 0 and  -«  < k"  < +» 

J J 

Finally,  using  Eo.  (All)  and  the  relation  6 = T(j>,  one  obtains 
the  following  expression  for  the  outgoing  ^ave  Rreen's  function. 


5 = T 


4irR 


(A12) 
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This  appendix  has  two  objectives.  The  first  is  to  show  that 
the  expressions  for  the  fields  for  both  open  resonators  and  segmented 
boundary  resonators  satisfy  Maxwell's  equations.  The  second  is  to  show 
that  these  fields  can  be  divided  into  m partial  fields  such  that 

1.  The  m^"  partial  field  involves  only  the  field  over 
the  mirror  surface  or  the  segment  of  a 
closed  boundary;  and 

2.  Each  partial  field  satisfies  Maxwell's  equations. 

The  formulation  and  manipulations  needed  to  accomplish  these  objectives 
are  covered  below. 


Formulation 

The  field  expressions  for  the  two  cases  described  above 
(Eqs.  (28)  and  (29),  and  Eqs.  (32)  and  (33))can  be  summarized  by  the 
following  equations. 


• v'xS  + (nxVxE„)  • g + (n-E^)v' 
' m ' m' 


(nxH„)  • v'x'^  + (nxv'xHJ*  5 + (n*H„)V  • 5 dS' 
m m'  m m 


/ 

'““  m / 


V (ji  E„  • ds„ 
m m 


where  S = !(().  For  closed  resonators,  the  sums  in  these  equations  are 
taken  over  all  continuous  segments  of  the  closed  interface  or_ boundary. 

For  open  resonators,  the  sums  are  taken  over  the  surfaces  of  the 
resonator  mirrors. 

To  obtain  these  two  equations,  it  has  been  assumed  that  both 
the  electric  and  magnetic  fields  are  square  integrable.  Thus,  the  points 
of  finite  discontinuity  have  simply  been  removed  from  the  above  integrals, 
and  the  remaining  surface  and  contour  integrals  are  considered  to 
extend  to  within  an  arbitrarily  small  distance  5 from  the  points  of 
discontinuity.  As  they  do  not  actually  include  the  points  of  disconti- 
nuity, the  integrands  (and  their  derivatives)  of  these  integrals  are 
continuous.  As  a result,  integral  theorems  such  as  Stokes'  theorem  can 
be  applied  to  the  surface  and  contour  integrals  in  the  above  two 
equations. 

To  accomplish  the  objectives  of  this  appendix,  the  following 
approach  is  used.  First,  the  fields  F(r)  and  FT(F)  in  Eqs.  (Bl)  and 
(B2)  are  written  as  the  sum  of  m partial  fields. 


E(F)  = S E„(F) 
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H(^  = E 

m 


(B4) 


where  r„(F)  and  FT  (F)  are  the  fields  result! no  from  the  m^^  seament  or 
m m 

mirror  surface.  Then  Eqs.  (B3)  and  (B4)  are  substituted  into  Eqs.  (Bl) 
and  (B2)  to  yield  the  pair  of  equations  shown  below. 


(nxE|^)  • v'x^  + (nxv'xF^)  • S + (n*rj^)V 


5 


dS 


m 


(B5) 


m m % , 

^m 


(nxH„)  • V'xf5  + (nxVxH„)  • + (n-H„)V 

m ' m m 


(B6) 


Then  it  will  be  shown  that  each  of  the  oartial  fields  E (F)  and  FT(F) 

m m 

satisfies  Maxwell's  eouations.  From  this  demonstration,  it  follows 
directly  that  the  total  fields  in  Eqs.  (Bl)  and  (B2)  satisfy  Maxwell's 
equations . 

To  show  that  these  fields  satisfy  Maxwell's  equations  in  all 
space,  it  is  necessary  to  assume  that  they  satisfy  Maxwell's  equations 
over  the  indicated  surfaces  and  contours.  Of  course,  any  assumption 
to  the  contrary  immediately  invalidates  the  assertion  that  the  modal 
fields  are  electromaanetic  in  nature.  Thus,  Eqs.  (B7)  and  (B8) 
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can  be  used  to  simplify  the  field  expressions  as  shown  below. 


■/ 


[(nxE*^)  • V X ^ 


io)u  (fi  X H„) 
m 


5 + (h*E^)v'*  5 


m 


looe  -L 


7'  (j.  • ds„ 

m m 


(1 


m 


Hn,(r) 

m 


■/I 


(nxH  ) 
m 


v'  X 6 + iue  (nxE  ) • ^ + (n*H  )v‘ 
m m 


m 


+ 


1 

iuu 


(I 


m 


Then,  subject  to  the  condition  that  these  fields  must  be  square  inte- 
grable  (which  was  assumed  earlier  in  this  aopendix),  the  objectives  of 
this  appendix  will  be  accomplished  if  and  aiven  above  are 

solenoidal  (divergence! ess)  and  satisfy  Eqs.  (B7)  and  (B8). 

Finally,  to  further  simplify  these  equations  and  to  aid  in 
later  manipulations,  the  following  vector  identities  are  listed. 

(7  X R)  • S’  = 7(}i  X S’  (I 

S’  • (7  X Pi)=  S’  X 7<J>  (I 

7 X [S’(r')  X 7'  (j>]  = S’(r‘)  - S’(f')7'2,*,  (I 


• 5 = ?T(|) 

A X W(ti  = (v  X ^V(j)  - V X (AV((i) 

V • (A'xf)=B*VxA‘-S'-vxS' 


Using  Eqs.  (B12),  (B14)  and  (B15),  Eqs.  (B9)  and  (BIO)  are  rewritten  as 
shown  below. 


■I 


(nxr  ) X v'(ji 

m 


- io)u  (nxHj^)  $ + v'(|)  (n-Ej^)  dSJ^ 


— fL  • ds„ 

Wt  I 


= / 
•'c  I 


(nxH^)  X V (|i  + iojG  (nxE^)  (p  + v'(j)  dSJ^ 


■? — V ‘f’ 

luy  / m m 


Manipulations 


To  show  that  the  fields  in  Eqs.  (B18)  and  (B19)  have  the 
desired  properties,  a procedure  similar  to  that  used  by  Stratton 
(ref.  33,  pp.  469-470)  will  be  employed.  One  first  takes  the  curl 
of  Eq.  (B18)  with  the  subscript  dropped  to  obtain 


V X [nxHr')  X v'<)>]  - icoyv  x [nxH(r' )](}.[  dS' 


where  7(|)  = - v'(|)  and  the  identity  v x V(t)  = 0 have  been  used  to  eliminate 
the  last  two  terms.  Using  Eg.  (B13)  and  the  fact  that  (t>  satisfies  the 
scalar  Helmholtz  equation  (Eg.  (A2)),  this  expression  is  simplified  to 
yield 


1. 1" 


V X rCr)  = - iwu  / |(nxff)  x 7'<|)  + iue  (oxF)  41 


dS' 


(nxE)  • v'v'(^dS' 


(B21) 


Using  Eqs.  (B7)  and  (B16),  and  the  fact  that 

(nxF)  • vV(p  »=  n • {ExvV((>} 

this  last  integral  can  be  simolified  as  shown  below. 


(nxE)  . v'v'4)dS' 


(n-H)  v'^  dS' 


[v'x(Fv'4))]  dS' 


(B22) 


(B23) 


The  last  integral  in  Eg.  (B23)  is  now  converted  to  a line  integral  by 
usino  Stokes'  theorem  applied  to  dyadics  (Collin,  ref.  34,  p.  569). 
The  result  is  shown  below. 


J n • [v'x(EV^)]  dS' 

S' 


(B24) 
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Then  substitutinq  Eqs.  (B23)  and  (B24)  into  Eq.  (B21),  one  obtains  the 
following  result. 


7 X r(F)  = - iojji  / (nxH)  X 7'iji  + iwe  (nxD'i' 


/ 7'())E  • ds 
C 


(B25) 


Comparing  the  right  side  of  Eq.  (B25)  with  Eq.  (B19),  one  can  see  that 
Eq.  (B7)  is  satisfied.  A similar  procedure,  which  is  not  shown  here, 
can  be  used  to  show  that  the  fields  in  Eqs.  (B18)  and  (B19)  satisfy 
Eq.  (B8). 

The  next  step  is  to  show  that  the  fields  in  Eqs.  (B18)  and 
(B19)  are  solenoidal.  One  begins  by  taking  the  divergence  of  Eq.  (B19) 
and  applying  Eq.  (ul7)  to  obtain 


V • H(F) 


■/I' 


(ijev  . [(nxE(r')(j)]  - v'24,(n»iT)  dS' 


-J- 

lOJW  J 


2(j)r  • ds 


Applying  the  Helmholtz  equation,  this  form  is  simplified  to  yield 


iwe  (hxF)  • 7'4i  + k^  (n*FT)  dS' 


Jsi  . 

iiou  Jr 


(B27) 


no 


Using  Stokes'  theorem,  the  line  integral  in  Eg.  (B27)  is  expressed  in 
terms  of  two  surface  integrals  as  shown  below. 


ds  = 


v'(|,  X E + (}i(v'xE) 


ndS' 


Using  Eq.  (B7)  and  the  fact  that 


(B28) 


(v'(t)xr)  • n = - (nxF)  • v'<t> 


(B29) 


the  line  integral  is  expressed  as 


(nx?)  • v'(j>  - ituy  (n*iT)  <(> 


dS' 


(B30) 


Substituting  Eq.  (B30)  into  Eq.  (B27)  and  using  = oi^ue,  one  can 
show  that  V • H"  = 0.  A similar  orocedure  can  be  used  to  show  that 


I 


APPENDIX  C 

Elimination  of  the  Self-Induction  Integral 

On  page  37  in  this  report,  it  is  shown  that  for  resonators 
containing  two  perfectly  conducting  mirrors  for  which  a/|R|  « 1,  the 
currents  on  mirrors  #1  and  #2  are  governed  by  the  following  pair  of 
equations , 


where 


>< 

< 

II 

^12‘^x2 

(Cl) 

'^22'^x2 

(C2) 

Vxt  = 2(-l)^/o,,(F')  ,^(r|r') 


r=r. 


'^qq'^xq 


= 2(-l)' 


■'xqf'-q) 


1 + ik  r - r 


r - r; 


<>(r'|rq) 


- ""q  • "q^ 


(C3) 


(CA) 


^(rqlr-)  = 


-ikkq  - r^l 


4tt  r_  - r' 

' q q' 


(C5) 


and  cos  (r  - r',  n ) is  the  cosine  of  the  angle  between  the  unit 

q q q 


normal  at  r and  the  vector  r„  - r'  It  was  further  shown  that  when 

q q q 

the  self-induction  integrals  (or  terms)  are  negligible,  that  is,  when 

IVxql  ""  I'^xql’  ^0^  q = 1.  2 (C6) 
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The  purpose  of  this  appendix  is  to  estimate  the  conditions  under  v/hich 
the  self-induction  integrals  are  negligible. 

To  obtain  this  estimate,  the  followina  approach  is  used. 

First,  it  is  assumed  that  the  resonators  being  considered  contain 
circular  mirrors  of  radius  a.  Then  the  inequality  in  Eg.  (C6)  is 
specialized  to  apply  to  the  center  of  mirror  #2  for  modes  which  have  no 
azimuthal  variation.  Finally,  a simplified  version  of  the  expansion 
functions  (see  chapter  V)  is  substituted  for  resulting 

integral  is  approximately  evaluated  for  mirrors  such  that  a/|R|  « 1. 
Thus,  the  estimate  itself  will  take  the  form  of  an  inequality  involving 
the  various  resonator  and  mode  parameters.  However,  since  this  inequal- 
ity is  merely  an  estimate  of  the  conditions  under  which  the  self- 

induction  integrals  are  nenliaible,  the  inequality  is  denoted  by 
■? 

the  symbol  <<. 

It  is  apparent  that  this  procedure  involves  many  approximations 
and  assumptions  concerning  the  resonators  considered  and  their  asso- 
ciated fields.  To  show  that  this  procedure  at  least  leads  to  the  right 
order  of  magnitude  of  the  conditions  under  which  the  self-induction 
terms  can  be  nenlected,  a second  aooroach  is  used  to  check  the  condi- 
tions for  one  particular  case.  This  second  approach,  which  employs 
the  method  of  steepest  descents  to  approximate  an  inteoral,  beoins 
with  Eq.  (C33). 
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w 


Using  the  first  approach  described  and  denotinq  the  current 
distribution  for  the  mode  by  u-(r),  one  specializes  the  inequality 

J 

in  Eq.  (C6)  to  obtain  Eq.  (C9). 


-ik.r'  h + ik.r' 

e I jcos  r^dr^ 


? 

<< 


Uj(0) 


(C9) 


The  expansion  functions  which  were  obtained  in  chapter  V,  have  the  form 


ik 


Uj(r)  = 9j(r)  e 


t ^ Vn2  - l'  r2 


(CIO) 


where  Re(k.)  = kl  > o,  and  q.(r)  is  assumed  to  be  slowly  varying  in 
^ ^ J ^ I 

-ik;r  ± - l'  r2 


comparison  to  the  exponentials  e ^ and  e 
Substituting  this  form  1 
inequality  shown  below. 


Substituting  this  form  for  u.(r)  into  Eq.  (C9),  one  obtains  the 

J 


9j(o)  I gj(r^)  exp  jikj  [-r'  i ^ j 


COS 


(Cll) 


j 


where  6„  = \o^  - 1 
9 

To  simplify  this  expression  further,  one  needs  to  approximate 
the  term  cos  (a^,rp.  To  do  this,  one  uses  figure  12  to  see  that  the 
vector  F2  can  be  expressed  in  the  form 
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where  R is  the  radius  of  curvature  of  mirror  #2  and  the  distance  z 
has  been  approximated  by  z = r«^/2R,  which  applies  if  jIt 


MIRROR  #2 


Figure  12.  The  Reometry  of  a Larne 

Radius  of  Curvature  Mirror 


Then  usino  Eq.  (C13) , 


®z  ’ ’"2  ®z  * ^2 
cos  (a  ,r')  = - — « 

2 ^ I M i I • 


in  conjunction  with  Eq.  (C12),  the  estimated  condition  is  simolified 
to  yield 
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a 

” m\j ["'2 ■ 


+ r-a' 

- 2L  ^2  . 


(1  + ikjr^)dr'  (C14) 


Since  the  integrand  contains  terms  of  the  form 
+k';r  +k';e.r2/2L 

e and  e ^ , the  probability  that  Eq.  (C14)  will  hold  (for 

both  + and  - signs)  should  increase  as  kj  decreases.  With  this  in  mind, 

and  to  simplify  the  analysis  as  well,  it  is  assumed  that  k'.'  is  suffi- 

+k!i'r  +k!i'8„r2/2L 

ciently  small  so  that  the  functions  e J and  e 4 ^ are  slowly 

varying  across  the  entire  mirror  surface.  With  this  assumption,  the 

magnitude  of  the  integral  in  Eq.  (C14)  will  be  largely  determined  by 

-iklr  ±ik^e  r^/ZL 

the  oscillatory  functions  e ^ and  e ^ .In  fact,  these  two 
oscillatory  functions  should  combine  to  produce  a significant  contribu- 
tion to  the  integral  in  regions  where 

1.  The  exponentials  oscillate  nearly  in  phase,  or 

2.  The  exponentials  are  nearly  conjugates 
over  some  significant  distance. 

There  are  two  types  of  regions  over  which  these  two  conditions 
appear  likely  to  occur.  The  first  type  consists  of  regions  surrounding 
points  where  the  two  phase  terms  are  equal;  that  is,  where 

k!e 

with  <()i(r)  a kjr  and  <l>2(r')  = 2L  course, 

‘t'-lCi")  = at  r = 0 and  at  r = 2L/8g. 

Although  the  phase  terms  are  equal  at  these  two  points,  their 
rates  of  change  are  not  equal  there.  In  fact,  at  r = o, 
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i 


dr 


kl  and 

V 


d(j)2(r) 

dr 


0 


while  at  r = — 


a 


9 


d(j)^(r) 

dr 


kl  and 


d(|)2(r) 

dr 


(C15) 


(C16) 


At  the  optical  or  infrared  wavelengths  of  interest  for  this  analysis, 
the  difference  in  these  two  rates  will  be  sufficiently  large  to  prevent 
the  exponentials  from  being  in  phase  (or  conjugates)  over  any  signifi- 
cant distance.  As  a result,  the  contribution  to  the  integral  from 
these  two  regions  (especially  the  one  near  the  point  r = o)  will  be 
small  in  comparison  to  the  contribution  from  a stationary  phase  type 
region  where  <Jii(r)  and  ij!2{r)  are  changing  at  nearly  the  same  rate. 

The  only  region  of  this  second  type  is  located  in  the  vicinity  of 
r = L/Bg.  To  estimate  the  width  of  this  stationary  phase  type  region, 
it  is  assumed  that  the  exponentials  remain  nearly  in  phase  or  nearly 
conjugates  over  the  distance  for  which. 


A(|>2(r)  - A(^^(r) 


(C17) 


where  the  changes  in  the  <|)^  are  computed  relative  to  their  values  at 
r = L/Bg.  Substituting  the  forms  for  and  <{i2(>')  into  Eq.  (C17), 
one  obtains  Eq.  (Cl 8)  below. 


^ r^ 
2L 


(C18) 


I 

i 

i 


n 
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This  equation  has  real  solutions  of  the  form 


(C19) 


where  S = X6g/2L.  Thus,  the  so-called  stationary  phase  region  will  be 
taken  as  the  region  extending  from  r = ^ (1  - \/?)  to  r = ^ (1  + 

g g 


To  estimate  the  magnitude  of  the  integral  (in  Eq.  (C14))  when 
the  mirror  radius  a is  sufficiently  large  to  include  all  or  part  of 
this  region,  the  following  procedure  is  used. 

1.  The  assumed  slowly  varying  functions,  g^(r)  and  the 
exponentials  involving  kV,  are  evaluated  at  the  midpoint 

V 

of  the  integration  interval  and  factored  out  of  the 
integral . 

2.  The  remaining  exponentials  are  taken  as  being  exactly  in 
phase  or  exactly  conjugates  over  the  integration  interval. 

3.  It  is  assumed  that  k^  » kV. 

J J 

Then  applying  these  steps,  taking  j |«l  g,-(o)  |>  and  assuming 

J Pg  J 

a > ^ (1  + -\/s^ , the  proposed  inequality  becomes 

^g 


k!.'L  k';L 

+ 

1 J e L§ 1 

2lRl 


/ 


^ (1  + V^) 

g 


iki 


I-  0 - 


(1  + ikjr^)  dr^ 


(C20) 
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JB. 


For  the  case  where  the  exponentials  are  in  phase  (upper  sign), 
Eq.  (C20)  takes  the  form  shown  below. 


3k';i 

J 

? ^ 9 


r 


-2iklr' 
e ^ 


r (1  -vH) 


-j/ 


i-n  *\Ts) 


e 1“  ^g'^2  ridr; 


By  determining  the  bounds  on  these  two  integrals,  one  can  show  that 
Eq.  (C21)  will  hold  if 


3k';L 

• 1 ^ 
^ I\R[  ® 


[2-  3gj 


If  the  exponentials  are  taken  as  conjugates  (lower  sign)  over  the 
stationary  phase  region,  Eq.  (C20)  takes  the  form 


Wl 


^ 0 + Vs) 


(1  + ikjrpdr’ 


|-(1  - ^^S) 
g 
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This  inequality  will  hold  if 


(C24) 


One  can  use  this  same  procedure  to  estimate  the  conditions 
under  which  the  self-induction  integrals  can  be  neglected  for 

T ® ^ (1  1 however,  for  this  case,  a is 

g g 

taken  as  the  upper  limit  of  integration  in  Eq.  (C20).  The  result 
is  two  rather  complicated  conditions  which  are  not  given  in  this 
report. 

This  same  procedure  is  also  used  to  estimate  the  integral  in 

Eq.  (Cl 4)  for  the  case  where  a < ^ (1  - Vs) . For  this  case,  it  is 

assumed  that  the  integration  region  extends  over  the  comparatively 
small  distance  (see  the  discussion  following  Eq.(C16))  corresponding 
to  the  condition, 

- A<^^(r)l  = ± f (C25) 

where  the  changes  in  (j)i(r)  and  computed  relative  to  their 

values  at  r = 0.  Application  of  this  condition  leads  to  the  equation 


(C26) 


120 


which  has  solutions  of  the  form 


r = | 1 ±>/l  +S 

9 


Since  r < ^ (1  - \[s) , one  chooses  both  minus  signs;  thus,  the 
region  of  integration  has  a width  w given  by  Eq.  (C28). 

w = (l  - >/l  - S%  (C28) 


For  the  case  when  S « 1 , w « x/4. 

To  at  least  partially  include  the  fact  that  the  exponentials 


+kVr  +kVr28  /2L 


e and  e 


J g 


increase  with  increasing  radius,  this  region  of 


width  w will  be  placed  at  the  outer  edge  of  the  mirror  of  radius  a. 
Thus,  for  this  case,  the  integral  in  Eq.  (C14)  will  be  approximated 
using  the  procedure  described  earlier  and  limits  of  integration  of 
r = a-w  and  r = a.  Then  evaluating  the  exponentials  involving 

k'!  at  r = a - j and  factoring  them  and  g-(r)  out  of  the  integral,  the 

J J 

proposed  inequality  takes  the  form 


j ' 


ik'. 

0 r '2 
2 2l  g 2 


(1  + ikj  rpdr' 


The  inequality  in  Eq.  (C29)  will  be  satisfied  if 
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+k!.‘(a  + ^ a2)  . 

1 


r 

w)  <<  1 


(C30) 


For  the  case  where  S « 1,  Eq.  (C30)  reduces  to  Eq.  (C31), 


+k"(a  + ^a2) 


li)  + <;  1 


(C31) 


which  will  be  satisfied  if 


(C32) 


For  future  reference,  the  estimates  of  the  conditions  under 
which  the  self-induction  integrals  may  be  neglected  are  summarized 
below.  First,  for  the  case  where  <<  1,  a < ^ (1  - V^) » 

X « — , this  estimate  corresponds  to  the  condition  shown  in  Eq.  (C32), 

“g 

For  the  resonators  considered  in  this  paper  (those  for  which  the  wave- 
length is  small  in  comparison  to  all  resonator  dimensions),  this 
condition  will  hold  for  all  but  the  lossiest  modes.  For  the  case  where 

« 1 and  j (1  this  estimate  corresponds  to  the  pair 

of  conditions  shown  below  and  referenced  for  your  convenience. 


1 „ ^'"g 


(C22) 
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2B 

L_  p ^ 

i:  FT  ® 


These  conditions  will  hold  only  for  cases  where  the  mirror  radius  of 
curvature  is  very  large  in  comparison  to  the  mirror  separation. 

As  indicated  earlier  in  this  appendix,  various  assumptions  and 
approximations  have  been  made  to  obtain  these  estimates  of  the  condi- 
tions where  the  self-induction  terms  may  be  neglected.  Many  of  tnese 
assumptions  were  made  to  simplify  the  integral  in  Eq.  (C14).  To  at 
least  show  that  this  simplification  process  has  led  to  the  right 
order  of  magnitude  for  these  estimates  (Eqs.  C22,  C24,  and  C32),  the 
integral 


/-IK. 


* a Vi!" 


(1  + ikjr')dr' 


will  be  approximated  for  the  case  where  < 3 using  the  method  of 

g 

steepest  descents  (Erdelyi,  ref.  41,  vol.  2,  pp.  24-27).  To  remain 
consistent  with  the  calculations  performed  earlier  in  this  appendix, 
it  is  assumed  that  k'.  » k'.'.  With  this  assumption,  I may  be  written 

V J 

in  the  form  characteristic  of  the  method  of  steepest  descents i that 


r kif(z) 

Iq  = J h(z)  e ^ dz 


(C34) 


123 


(C3 


h(z)  = gj(z)  (1  + ikjz)  e 
f(z)  = i ^2l 

and  z = x,+  iy. 

To  locate  the  saddle  point,  one  determines  the  point  at  which 

= 0 (C3 

For  this  problem,  the  saddle  point  is  located  at  z = L/Sg  and,  the 
steepest  descents  contour,  which  is  specified  by  the  equation 
Im  f(z)  = Imf(L/Sg).  corresponds  to  the  line  y = x - L/Sg  in  the 
complex  z-plane.  As  part  of  the  overall  procedure  to  estimate  1^, 
this  steepest  descents  contour  (C^^)  has  been  chosen  as  one  segment 
of  the  closed  contour  shown  below  in  figure  13,  where  C-j  and  are 
lines  parallel  to  the  imaginary  axis.  begins  at  the  point 
z = a and  ends  at  z = a(l  + i)  while  begins  at  z = -i  L/s^  and 
ends  at  z = 0.  The  segment  extends  along  the  real  axis  from 


X = 0 to  X = +a. 


Figure  13.  Closed  Contour  Chosen  to  Include  C 


Sd 

By  applying  the  theory  of  residues  to  the  closed  contour 
shown  in  figure  13,  one  obtains  Eq.  (C38) 


where 


+ 12  + 


^Sd^ 


(C3b 


I 


h(z) 


+k;f(z) 

e 


dz 


(C3y 


and  so  forth.  However,  by  writing  out  the  expressions  for  and 

and  using  the  fact  that  >>  k'.',  one  can  show  that  these  ii.tegrals 
involve  exponentials  which  decay  very  rapidly  with  increasing  distance 


from  the  real  axis.  As  a result,  [1^1  and  II2I  should  be  small  in 
comparison  to  Therefore,  they  will  be  neglected  for  the 

remainder  of  this  analysis. 

Finally,  to  simplify  the  analysis,  only  the  first  term  in  the 
asymptotic  expansion  of  (Erdelyi,  ref.  41,  vol . 2,  p.  26)  will 
be  used.  One  thus  obtains  the  following  expression. 


+iLK'7er 


I ~ /HE!  q /k  \ e ^ J 9 ^ U ■ Bg 

4d~Vk-eg 


(C40) 


Taking  the  bound  of  the  right  side  of  Eq.  (C40) , the  following  condi 
tion  is  obtained. 


1 + 2ir 


(C41) 


Substituting  Eq.  (C41)  into  Eq.  (C14)  and  taking 


9j(o) 

’j  (t) ' 


(C42) 


the  proposed  inequality  for  the  case  where  L/B^  < a takes  the  form 


Vs  3[i*2.^]  <? 


(C43) 
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Comparing  this  expression  with  Eq.  (C24) , one  can  see  that  the 

1 


two  conditions  differ  by  a factor  -j=p 

V2 


Although  this  small  difference 
does  not  prove  the  validity  of  the  inequalities  in  this  appendix,  it 


does  indicate  that  they  have  the  right  order  of  magnitude. 


APPENDIX  D 


Detailed  Calculation  of  Expansion  Functions 


This  appendix  presents  the  detailed  calculations  of  the  expansion 
functions  used  in  conjunction  with  the  Rayleigh-Ritz  procedure  discussed 
in  the  body  of  this  report.  These  calculations  are  performed  for  reso- 
nators containing  toroidal  as  well  as  spherical  mirrors.  The  spherical 
mirror  cases  are  considered  first. 


Spherical  Resonator  Calculations 

Spherical  Mirrors  of  Rectangular  Projection.  As  indicated  in 
chapter  V (page  56),  the  expansion  functions  used  in  computing  laser 
resonator  modes  are  approximate  solutions  to  an  integral  equation  for 
a symmetric  laser  resonator.  The  particular  symmetric  resonator 
chosen  consists  of  two  mirrors  identical  to  the  one  on  which  the 
current  is  being  analyzed.  The  mirror  separation  is  the  same  as  in 
the  original  problem.  Thus,  to  determine  the  expansion  functions  for 
the  current  on  mirror  #2  with  radius  of  curvature  R and  transverse 
dimensions  2a  and  2b,  one  begins  with  Eq.  (71)  which  is  repeated  below. 


'^x2^'"2^ 


-^1 


''xl'-'i* 


where  R21  is  the  distance  between  two  points  on  mirrors  #1  and  #2, 

021  is  the  angle  between  R21  and  the  optic  axis,  and  Sj  is  the  surface 


of  mirror  #1 . 


However,  as  this  equation  is  to  be  applied  to  a symmetric  resonator. 


where  n'  = ivq  and  q is  an  integer.  One  can  see  this  by  choosing  the 
origin  midway  between  the  two  mirrors  and  realizing  that  the  fields 
must  be  even  or  odd  with  respect  to  z.  Then,  writing  Eq.  (Dl)  in  terms 
of  rectangular  coordinates,  using  Eq.  (D2)  to  eliminate  and 

dropping  the  prime,  one  obtains  the  following  result, 

• t f f -ikR2i 

Jx2(^2’^2^  " ■ ^ ® J J '^x2^^r-^l^  cosa2^dx^dy^  (03) 

where  the  projection  of  R21  in  the  x-z  plane,  denoted  R^-j , is  shown  in 
figure  14. 


Figure  14.  Geometry  of  a Spherical  Resonator 
with  Rectangular  Mirrors 
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The  distance  R21  is  given  by 


R|l  = (x^  - X2)2  + (y^  - y2)^+  (z^  - 22)^ 


(D4) 


To  simplify  the  procedure  for  obtaining  the  expansion  functions, 
it  is  assumed  that* 


« 1 


(D5) 


(Do) 


k"  <<  k' 


(D7) 


d_  « L 
m 


(D8) 


where  d^  = a^  + b^.  With  these  assumptions,  cosa„,  « 1,  and 
m 21 


-ik'R 


igk' 

21_^-ik'L  ^ 2L 
~ e e 


(xf  + x|  + yf  + yp  (x^X2  + y^y2)  (D9) 


where  g = 1 - L/R  and  the  first  two  terms  of  the  binomial  expansion  have 
been  used  to  approximate  R21 . Finally,  it  is  assumed  that 


+k"L 

e 


(DIO) 


*For  resonators  for  which  Eqs.  (D5)  through  (D8)  are  not  valid,  the 
effect  of  these  approximations  will  be  simply  to  require  the  use  of  a 
relatively  large  number  of  expansion  functions  to  represent  the  modes. 

tin  general.  <J2  - (xf  t 


across  the  entire  mirror  surface.  Using  these  approximations,  Eq.  (D3) 
can  be  manipulated  to  yield 


g+itJ ' - ikL 


(Dll) 


(D12) 


(D13) 


where  5,  = x . , n,  = T ' “x  ” VF  "y  ’ 
Jj^^(x.y)  = u(x)v(y). 


Since  Eqs.  (Dll)  and  (D12)  are  uncoupled  equations  of  identical 
form,  the  remaining  calculations  are  performed  considering  only  the 
x-variation.  To  simplify  these  calculations,  one  lets 


K(52Ui) 


(DU) 


so  that  Eq.  (Dll)  can  be  written  in  the  form  shown  below. 


u(C2)  = / u(Ci)  K(C2Ut)  d5^ 

-H.. 


(Dlb) 
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In  the  next  step,  one  determines  the  operator  M (see  pages 
57  - 62)  associated  with  the  eigenvalue  problem  for  the  expansion 
functions.  To  do  this,  one  begins  with  Eqs.  (D16)  and  (D17)  below. 


d^u(C2) 


c12K(52Ui) 


dS] 


(D16) 


d5iK(C2l5i) 


d^u(?^) 

d^^ 


u(^l) 


dcf 


(D17) 


where  Eq.  (D16)  was  obtained  by  computing  d^/dc|  of  Eq.  (Dl5),  and 
Eq.  (D17)  was  obtained  by  integrating  the  term  K(C2U-|)  d^u(5-|)/dc^  by 
parts  twice.  The  function  Ru(C-|)  in  Eq.  (D17)  is  defined  by  Eq.  (018) 
below. 


Ru(s^) 


u(5t) 


dK(52Ui) 


d5i 


K(52Ui) 


du(e^) 

dCn 


(018) 


Then,  evaluating  the  derivatives,  substituting  the  expression 
for  d2K(C2l5p/d^|  Jnto  Eq.  (016),  and  using  Eq.  (015)  to  simplify 


terms  of  the  form  0(5^)  K(C2l5i)  -Qne,^c^^  obtain  the 


integrodifferential  equation  shown  below. 


1 . 
— + + ig 

dc|  ^ 


u(5o)  = 


2^  - ■ 


r ^ 

/ u(?,)  K(C2Ui)  Sfd^i 

-H.. 


r ^ 

+ 2gC2  / C-|U(5i)  K(?2Ui)  d5^  (D19) 

-H.. 


A similar  procedure  involving  d^K(52l5i)/<^5]  Eq.  (D17)  yields 
Eq.  (D20). 

H 

-{C|  + ig)  u(52)  = K(C2Ui; 


) 


Si-  t 

dSf 


ufC^)  dc^ 


;> 

■ 2«2’'x  y 

-H.. 


5-|U(Ci)  K(52Ui)  dc^  + Yj^Ru(5i)  (D20) 


One  then  adds  Eq.  (D19)  to  Eq.  (D20)  to  obtain 


where 


Mu(C2) 


K(C2l5i)Mu(5-,)dC^  - Yj^Ru(Ci) 


(021) 


M = ^Mg2  - 1)  (022) 

is  the  desired  operator.  Using  this  operator,  the  eigenvalue  problem 
is  formulated  as  the  following  differential  equation. 
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+ (g^  - 1)5^  + u(5)  = 0 (D23) 

) 

To  solve  this  equation,  one  first  considers  the  case  where 

g2  > 1 and  makes  the  substitutions  d^  = 2 yj(g^  - 1)  and  z = ^ 5d 
to  obtain  Eq.  (D24), 

^^^+-4+(v  + i)  u(z)  = 0 (D24) 

dz^  ^ ^ 

. 

where  v + 1/2  = - is^/d^.  This  equation  is  Weber's  differential  equation 
(Whittaker  and  Watson,  ref.  39,  p.  347),  which,  if  v is  not  an  integer, 
has  a general  solution  of  the  form 

u(z)  = AD^(z)  + BD^(-z)  (D25) 

where  the  Dy(z)  are  parabolic  cylinder  functions  (Lebedev,  ref.  40, 
chapter  10) . 

By  letting  B = ± A,  these  solutions  are  specialized  to  apply 
to  either  odd  or  even  modes.  Making  this  substitution  and  writing  the 
D^(z)  in  terms  of  the  confluent  hypergeometric  functions  of  the  first 

A 

kind  (Erdelyi,  ref.  41,  vol . 2,  p.  123),  one  obtains  the  following 
expressions. 

gygp  modes 

Ae"^  z$  (ij^»  modes  (D<^6) 

At  this  point,  the  only  remaining  step  is  to  specify  the 
parameter  v by  following  the  procedure  outlined  on  page  62.  That 
procedure  corresponds  to  the  requirement  that  as  g -►  1 , the  expansion 
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functions  must  reduce  to  the  plane  parallel  solutions  of  Vainshtein 


shown  below. 


u(5)  = 


cos  S5  for  even  modes 
sin  S5  for  odd  modes 


(D27) 


where  s = mTr/(M  + (1  + i )g)  > 6 = -c(l/2)/-(/i' « 0.824,  M = \Z8irN',  and  N 
is  the  resonator  Fresnel  number.  For  even  modes,  m = 1,3,5 — and  for 
odd  modes,  m = 2,4,6 — . 

This  procedure  for  specifying  the  parameter  v can  be  simplified 
considerably  by  realizing  that  as  g 1 , d^  -►  0,  and  [vj  Letting 

a = - v/2  or  (1-  v)/2  and  y = z^/2,  this  procedure  also  corresponds  to 
investigating  the  behavior  of  the  $(a,Y,z)  as  |a|  ®.  Several  in- 

vestigations of  this  type  nave  been  performed. 

To  obtain  the  particular  form  used  for  one  such  investigation, 
one  sets  < = y/2  - a,  and  requires  ay  to  be  bounded  in  absolute  value. 
Then,  since  y = 1/2  or  3/2,  <y  is  also  bounded  in  absolute  value.  With 
these  restrictions,  it  has  been  shown  that  (Erdelyi,  ref.  41,  vol . 1, 
p.  280) 


1 - Y rj 

Y.  y)  t(y)  (<y)  ^ ^ i(2v^Ky) 


(D28) 


One  then  uses  the  fact  that  since  < = (v/2  + 1/4)  = - is^/2d^, 

<y  = s^c^/4.  This  expression  for  <y  is  then  substituted  into  Eq.  (028) 
to  obtain  Eq.  (029) . 


Y id^c^ 

e 4 J 
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Then,  for  the  even  mode  solutions  of  Eq.  (D26),  one  lets 


Y = 1/2, 

u(c)  -*■  A r(4)  (sc) 

g 1 ■'* 

and  for  odd  modes,  y = 3/2 

u(5)*^  Msr  J^(S5) 

Substituting  the  well-known  forms  (Whittaker  and  Watson,  ref.  39, 
p.  364) 


’ ■'/S’ 


J_V  (S5)  = “S  « 


into  Eqs.  (D30)  and  (D31),  one  finds  that  for  even  modes. 


u(c)  -►  A cos  sc 
g ^ 1 


and  for  odd  modes, 


u(c) 


g 


A . 

— ^sin  sc 


Thus,  for  the  desired  reduction  to  occur,  it  must  be  true  that 
s = mTr/(M  + (1  +i)3).  Therefore,  v has  the  value  given  below. 

1 

I 

I 


1 

2 

1 + 2i 

mWd  1 

I 

M + TT  + i)6  J 
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To  complete  the  procedure  for  obtaining  the  expansion 


functions,  one  returns  to  the  differential  equation,  considers  the 


case  < 1 , and  makes  the  substitutions  h^  = 2-^'l  - g^'  and  y = h5. 


The  result  is  Weber's  equation  as  shown  below. 


j(.  4)  - 4j  „(y)  = 0 


For  all  values  of  v,  the  general  solution  of  this  equation  may  be 


written  in  the  form  (Whittaker  and  Watson,  ref.  39,  p.  348), 


u(y)  = AD^(y)  + B0_^  _-i(iy) 


However,  for  real  values  of  y,  D_^  _1  (iy)  increases  exponentially  as 
[y|  increases.  Since  this  directly  contradicts  the  known  behavior  of 


the  modes  of  stable  resonators,  one  requires  that  B = 0.  Finally,  the 


requirement  that  the  modes  be  either  even  or  odd  yields  v = 0,2, 4, 6, — 


for  even  modes  and  v = 1,3,5, — for  odd  modes.  Thus,  for  g^  < 1 , the 


expansion  functions  are  given  by  Eq.  (D39). 


u(y)  = 


(.n  I 

\ 2’  2’  2/ 


for  even  modes,  n = 0,2,4, — 


f.  ^ 


for  odd  modes,  n = 1,3,5, — (039) 


Then  combining  the  forms  in  Eqs.  (026)  and  (039),  the  expansion 


functions  for  one  transverse  dimension  are  given  by. 


u(0  = 


j, 


1 h2£2 


for  even  modes 


2 2 

-H- 

56  $ 


for  odd  modes 


(040) 


J 
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where  = 2 \ l-  g^  and  5 = V x.  For  stable  resonators, 

V = 0,2,4, — (even  modes)  or  v = 1,3,5, — (odd  modes).  For  unstable 
resonators,  v is  given  by  Eq.  (036)  where  m = 1,3,5, — (even  modes) 
or  m = 2,4,6, — (odd  modes). 

Spherical  Mirrors  of  Circular  Projection.  To  obtain  the 
expansion  functions  for  spherical  mirrors  of  circular  projection,  one 
begins  with  the  integral  equation  for  a symmetric  resonator  in  circular 
cylindrical  coordinates  as  shown  below, 

. -2tt  a -ikR.,2 

Jxi(pi » ®i)  = - ® / ^^2  j Jxl^‘^2’  ®2^  ^ R,p 

■•'o  0 

C0Sai2P2^P^  (D41J 

where  n'  = irq  as  before.  Then,  in  addition  to  assuming  that  Eqs.  (05) 
through  (010)  hold,  one  assumes  solutions  of  the  form 

0x2(p.  e)  = ^(P)  (D42) 

and  applies  the  identity  (Erdelyi,  ref.  41,  vol . 2,  p.  7), 

2tt 

1"  2.J„(z)  -f  d»  (D43) 

0 

where  J„(z)  is  a Bessel  function  of  the  first  kind.  The  result  is  the 
n 

following  integral  equation  for  the  radial  mode  function  i^p(c)» 

r®  - (Ci  + Co) 


(04 


where  = -^k'/L  p^ , and  g = 1 - L/R.  Finally, 

one  converts  this  equation  to  an  integral  equation  with  a symmetric 
kernel  by  making  the  substitution  u^(c)  = v^(c)/vT-  The  result  is 

shown  below  in  Eq.  (D45), 


"a 

= ’'n/  ''n‘«2>  '*«2 


(D45) 


where 


_ i£  (^2  + ^2) 

K„(C,I«2)  Jn<«l«2) 


(D46) 


To  obtain  the  eigenvalue  problem  for  the  expansion  functions, 
one  begins  with  the  following  pair  of  equations. 


' 0 ' 


d\(52l5,) 


(D47) 


"a 

^n| 


d2vn(C2) 


d52  = 


■ 1 Hr 

•2^  , ^^2 

^ dd 


where 


+ v„Rv„(52> 


dv  (So) 

Rv„(52>  = K„(52l5,)  - v„(«2> 


(D4d) 


dKn(52Ui) 


(D49) 
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The  indicated  derivatives  of  •^p(C2U-|)  tlien  evaluated 
and  a procedure  identical  to  the  one  used  for  rectangular  mirrors 
followed  to  obtain  Eqs.  (D50)  and  (D51)  shown  below. 


dsf 


(t  ■ "^)l 

+ g^cf  + 2ig  + ^ ^ 


'1 


0 

»a  3 


- 2i9  y„f  (5,52)^J;(5,52>  E(52l«l>  '*«2 


{21g  + Sf}v„(5,) 


- Y 


n / K„(52l5,) 

0 


a a & ■ "')i 
- — + g^5o  + ' 

dc|  ^2 


Vn(52)dC2 


+ 2ig 


(51^2)^  J;(5iC2)  E(C2Ui)  dS2  - YpRv^(52) 


where 


E(C2Ui)  = e 


if  (c|  * cf) 


is 


(D50) 


(Dol) 


(062) 


Finally,  one  adds  Eqs.  (D50)  and  (D51)  to  obtain  Eq.  (Db3). 


(053) 


One  can  see  that  this  equation  has  the  same  form  as  Eq.  (113)  in 
chapter  V.  However,  Eq.  (053)  will  not  have  the  characteristic  de- 
scribed in  condition  #1  following  Eq.  (122)  unless  the  term  in  braces 
in  Eq.  (049)  vanishes  at  ^2  ~ 0*  Since  the  functions  ^^^(52)  liave 
not  yet  been  selected,  it  is  initially  assumed  that  this  term 
vanishes  at  the  origin.  The  validity  of  this  assumption  is  demon- 
strated later  in  this  appendix  (beginning  with  Eq.  (077)). 

Then,  subject  to  the  above  assumption  and  in  accordance  with 
Eq.  (125),  the  operator  M is  chosen  as  shown  below. 


M 


(054) 


Thus,  the  eigenvalue  problem  for  the  expansion  functions  is  expressed 
as  the  following  differential  equation. 

/ 


d"v„(c) 

d4^ 


(g2  - 1)  52 


(i  ■ •■) 

+ ' ' + 


v,(c)  = 0 


(055) 
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To  solve  this  equation,  one  first  considers  stable  resonators 


(g^  < 1)  and  makes  the  substitutions. 


w„(5)  = v„(5)  ^ 


Z = 


(D56) 

(D57) 


J \ 

where  a'  = - 9^-  Then,  letting  s^  = - (2v  + l)a',  one  obtait 

Eq.  (D58). 


1 , . „„U)=0 


(D58) 


Equation  (D58)  is  Whittaker's  equation  (Whittaker  and  Watson,  ref.  39, 
p.  337),  which  has  a general  solution  of  the  form 


w „ (z)  = A W Iz)  + 8 W „ (-z) 

K,P  ''  ' <,  p -K,  p ' " 


(059) 


where  k = - 1/4  (2v  + 1),  p = n/2,  and  W (z)  is  a Whittaker  function 

• P 

of  the  second  kind.  Now,  the  Whittaker  functions  are  related  to  the 
confluent  hypergeometric  functions  of  the  second  kind,  'i{o.,y,z),  by 
the  following  equation  (Lebedev,  ref.  40,  p.  274), 


. z Y 

W (z)  = e ^ -{-(a,  y,  z) 

^ -5- 


(060) 


where  a = (n  + 1)/2-k  and  y = n + 1.  When  n is  zero  or  a positive 
integer,  which  it  is  for  this  problem,  'i{a,  n + 1 , z)  is  represented 
by  the  series  (Lebedev,  ref.  40,  p.  264), 


/ i\n  ^ 

4'{a,n+1,z)  = (n+k)'ik!  ’!'(“'*‘k)-<l^(1+k)-i()(n+k+1  )+1n(z) 


*rW  S 


^ (-1  )*^(n-k-1 ) ! (a-n) 


(D61) 


where  (a)j^  = r(a+k)/r(a),  and  i|)(x)  is  the  logarithmic  derivative  of  the 
gamma  function  (Erdelyi,  ref.  41,  vol . 1,  p.  15). 

From  this  series,  one  can  see  that,  unless  a = - m where  m 
is  zero  or  a positive  integer,  'V{a,  n + 1 , z)  has  a singularity  at 
z = 0.  However,  if  a = -m,  the  series  is  indeterminant,  and 
v(a,  n + 1,  z)  must  be  evaluated  using  a limiting  process.  For  this 
case,  one  finds  that  (Erdelyi,  ref.  41,  vol.  1,  p.  268) 


'l'(-m,  n + 1 , z)  = (-l)"’m!  Fj^"(z) 

where  f^^'^(z)  is  a generalized  Laguerre  polynomial.  Thus,  if  the 
origin  (z  = 0)  is  included,  the  solution  to  Eq.  (D58)  becomes 

• A'e  c"  * ' F^" 


t B'  e e"  * ’ F^"  (-a'£2) 


(D62) 


(D63) 


However,  since  stable  resonators  do  not  produce  current  distri- 
butions for  which  the  current  grows  exponentially  with  increasing 
miroDt.radius  (Bergstein,  ref.  2,  p.  500),  B'  m^^  b^^t 
zero.  Thus,  the  expansion  functions  are  given  by  Eq.  (064). 
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w„„,(c)  = A*e  “2~^n  + 1 (064) 


For  unstable  resonators  (g^  > 1),  it  is  advantageous  to  convert 
the  differential  equation  back  to  an  equation  for  To  do  this, 

one  substitutes  = v^(5)/vT  in  Eq.  (D55)  to  obtain. 


d\(z)  ^ 1 du„U) 

dz^  ^ dz 


U„(2)  = 0 


(D66) 


where  z = n?,  t = - (s/2fi)^,  and  = 1/2  ^g^  - V.  This  equation 

A 

has  a general  solution  of  the  form  (Erdelyi,  ref.  41,  vol . 2,  p.  1Z6), 


u 


n,T 


AM  (iz^)  + B W ^ (+iz2) 
ix,'^  ix,'^ 


(D66) 


where  the  M and  W are  Whittaker  functions  of  the  first  and  second 

K,U  K,U 

kind. 


As  with  rectangular  mirror  resonators,  one  must  now  insure 
that  these  solutions, 

1.  Have  no  singularities,  and 

2.  Reduce  to  the  proper  forms  as  g 1 (n  0). 

Dealing  with  the  singularities  first,  one  uses  the  fact  that 

Wix  n/2^^^  ® singularity  at  the  origin  unless  ix  = - (m  + (n  + l)/2), 

where  m = 0,1,2,-—.  For  this  case  (Lebedev,  ref.  40,  p.  274), 

iz^  n + 1 

M ^ (iz^)  = e (fz^)'  *^  ■ '$((11  + n •■+•1 ,-  n + t,  '^iz2)-  * (067) 
ix , j 
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W ^ (+iz2)  = e ‘‘  (+iz^)  4'(m  + n + 1 , n + 1,  +iz^)  (D68) 

i T j "2’ 


and 


u 


n.T 


+ B' 


A'e  $(8,  6,  iz^) 
- ( 

e fCs,  6,  +iz^) 


(D69) 


where  8 = m + n + 1 and  6 = n + 1 . However,  as  -*■  0,  z -►  0 and 
't(s’.  iz^)  -*•  1 and  J,  +iz2)  1 (see  appendix  H).  Thus,  as 

-*•  0,  the  solution  u^(z)  approaches  a constant  value,  which  does  not 
correspond  to  one  of  the  solutions  (for  the  plane  parallel  case) 
shown  below. 


(V  p/a  \ 

rrjhr-ijm) 

where  is  the  m root  of  the  n’^^  order  Bessel  function  of  the  first 
kind,  8 = 0.824,  M = yjs-nH,  and  N is  the  resonator  Fresnel  number.  As 
a result,  the  case  where  1t=  - (m  + (n  +i)/2)  will  not  be  considered 
further.  Therefore,  to  avoid  the  singularity  at  the  origin,  one 
must  set  B = 0 in  Eq.  (D66)  to  obtain 

u^  ^ (z)  = AM_  ^ (iz2)  1 (D71) 

’ i T , •2‘ 

where  the  yet  to  be  determined  values  of  t correspond  to  those  for 


which  the  u (z)  reduce  to  the  functions  in  Eq.  (070)  as  Q 0. 


From  the  relation  t = - (s/2fi)2,  one  can  see  that  as  a -*■  0, 


lixl  «.  Thus,  as  with  the  rectangular  mirror  case,  it  is  necessary 
to  investigate  the  behavior  of  the  solution  as  a parameter  becomes 
large.  To  do  this,  one  first  writes  u (z)  in  terms  of  $(a,  y,  iz^) 
as  shown  below  (see  appendix  H). 

n + 1 _ iz^ 

(-i)  ^ ^ M ^ (iz2)  = Ae’  z"  - ix,  n + 1,  iz^  j (D72) 

ix,  ^ ' 


Then,  using  Eqs.  (D71)  and  (D72)  in  conjunction  with  Eq.  (D28) , where 
a=  Cn  + l)/2+i  (s/2n)2,  Y = n + 1 , < = y/2  - a,  and  <y  = Kiz^  = (sc/2)^, 
one  can  show  that  if  |Kyl  is  bounded,  then 


n^ 

u„  (5)  ^ A j i (ix)  ^ n!  ) 0^(5?) 
n,x  g ^ T * in 


(D73) 


Thus,  the  u (c)  reduce  to  the  desired  form  if 
n ,T 


fk^  c ""nm 

' at  1 t iVJli 


(D74) 


where  $ and  M were  defined  in  the  discussion  following  Eq.  (D27). 
Eq.  (D74)  corresponds  to  a value  of  x such  that 


X = 


nm 

M + (1  + i)6 


_[ 


(D75) 


For  resonators  containing  mirrors  with  central  coupling  apertures 
this  expansion  set  must  be  altered  somewhat.  The  alteration,  which 
results  from  the  fact  that  the  singularity  at  the  origin  is  no  longer 
present,  simply  corresponds  to  retaining  the  terms 


6 


, 


t ' 

f:  ! 

S I 


(iz2) 


(D76) 


in  Eq.  (D66)  with  t given  by  Eq.  (D75). 

In  the  discussion  just  following  Eq.  (D53),  it  was  assumed  that 


= 0 
^2  = 0 


(D77) 


where  l^p(C2l5i)  is  given  by  Eq.  (D46)  and  the  v^(S2)  a>'e  the  expansion 
functions  for  resonators  with  spherical  mirrors  of  circular  projection. 
Using  Eqs.  (D64)  and  (D71),  one  can  show  that  for  stable  resonators. 


^n„  U2)  = A'e^C2""?  C 


(D78) 


and  for  unstable  resonators. 


''n.,  («2) 


_ 1 

B'52^M  n (1^2^52^ 
if , ■2' 


(D79) 


Further,  by  using  a relation  (see  appendix  H)  between  n/2^^^ 
and  $(a,  y,  z),  Eq.  (D79)  can  be  rewritten  in  the  form  shown  below. 


^ .in2^2 

''n,T  ^^2^  " B'(ifi2)^C2^e  ^ - ix,  n + 1,  (D80) 


Since  the  functions  (a'5|)  and  $((n+l)/2  - ix,  n+1 , 
are  entire  functions  of  52*  the  expansion  functions  of  both  stable  and 
unstable  resonators  can  be  summarized  by  the  form 
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Vn(C2)  = ^2  ^(^2^ 


(D81) 


where  ^{^2)  is  an  entire  function  of  ^2^.  The  kernel,  ^^(^21^1^’ 


written  in  a similar  form  as  shown  below, 


'^n(«2l«1>  ' «2 


(D82) 


where 


— - (£^  + £^1 
G„(«2l5,)  = 1,^  e ^ ^ 


(D83) 


is  an  entire  function  of  ^2- 

To  show  that  Eq.  (D77)  holds,  one  differentiates  '^^^(^2^ 
Kn(52Ui)  using  the  forms  in  Eqs.  (D81)  and  (082)  to  obtain  the 
following  pair  of  equations. 


^^n^^2^  _ 1 ■ I , . 2 

F(52)+52 


(084) 


(085) 


Substituting  Eqs.  (084)  and  (085)  into  the  term  in  brackets  in 
Eq.  (077),  one  obtains  the  following  result. 


n ■ £ 


Uo) 


n c.  • I 


d5 


= Sc 


GpU2Ui) 


2 

dF(5o) 


dSc 


- F(Co) 


dGnU2Ui) 


dc. 


(D8i 


Since  the  functions  F(C2)  and  G^(52Ui)  are  entire  functions  of  ^2.  the 
term  in  braces  in  Eq.  (086)  is  finite  at  ?2  ~ G-  Thus,  Eq.  (077)  holds 
as  assumed. 


Toroidal  Resonator  Calculations 

As  with  spherical  mirror  resonators,  one  begins  this  procedure 
for  obtaining  expansion  functions  with  the  integral  equation  for  a 
symmetric  resonator  as  shown  below, 

b Zt[  -ikR-|2 

Jj^l(Pl,  6i)  = - ^ e'^^^  J J J^i(p2,  62)  ® r"Y" 

a 0 

cosa^2p2‘^‘^2^®2 

where  R^2  is  given  by 

2 

R-| 2 ~ P^  P2  ” 2piP2COs(02  “ S'! ) ■*■  (^2  ~ 


and  a and  b are  the  inner  and  outer  mirror  radii  as  shown  in  figure  lo. 


Figure  15.  Intersection  of  a Toroidal  Mirror 
With  the  x-z  Plane 


The  distance  is  equal  to 

(Zg  - z^)  = L - - A2  (D89) 

where  the  distance  A^. , which  is  depicted  by  A in  figure  15,  is  given 
by 


Pe)' 


2R 


(D90) 


where  R is  the  radius  of  curvature  of  the  individual  arcs  shown  in 
figure  15  and  is  the  distance  the  axis  of  each  of  these  arcs  is 
displaced  from  the  optic  axis.  Substituting  Eqs.  (D89)  and  (D90) 
into  Eq.  (D88)  and  neglecting  terms  of  second  order  in  A^. , one  obtains 
the  following  expression  for  R^2' 
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I ■"  99%^ 


R^2  « L ( 1 


[(p^  - Pg)^  + (P2  - 


LR 


P-|^  + P2^  2p-|P2 


1 

II 


cos  (e,  - ej 


X2 '■''2  1 

Then,  in  addition  to  assuming  that  Eqs.  (D5)  through  (DIO) 
hold,  one  assumes  solutions  of  the  form 


(D91) 


0^1 (p, 9)  = u^(p)  e 


+ine 


(D92) 


and  applies  Eq.  (D43)  to  obtain  the  following  equation  for  u^(p), 
b 


^n  T J '^n(r  ‘'l^z)* 


j(Pl-Pe)"  + (P2-Pe^"j 


K f 1^  ^l) 


P2^P2  (D93) 


, _ ^n+1  ,-ikL+in' 


where  i e 


Equation  (D93)  is  simplified  considerably  by  making  the  sub- 
stitution result  is  shown  in  Eq.  (094), 


Hb 

%(«!>* %(52>  ® 

*Ti. 


4 ^ 


+ i(g  - l)5g(£i  + C2)  ^ 

2 ® ^2^52  (D94) 
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■9! 


where  = \ k'/L  a,  e , and 

g = 1 - L/R.  Finally,  to  simplify  the  forthcoming  computations, 
the  following  definitions 


E(52Ui)  = e 

are  employed  to  yield  Eq.  (D97) 


- l)5e(?i  * 52) 


= j'  %< 


“n(5l)  = \j  %<«2>  “«2 

“h. 


(D96) 

(D96) 


(D97) 


The  manipulations  used  to  obtain  the  expansion  functions  for 
these  toroidal  resonators  are  somewhat  more  complex  than  the  manipula- 
tions already  performed  for  spherical  mirror  resonators.  For  this 
case,  one  begins  with  the  following  two  equations. 


fUalfll  + L.  fliniill 

dcf  h '^'1 


u„(52> 


, 1 d 

dEf 


K„(52l5,)  152 


(D98) 


(D99) 
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where  Eq.  (D98)  was  obtained  by  operating  on  Eq.  (D97)  with 
d^/d^i^  + 1/5^  d/d5^ , and  Eq.  (D99)  was  obtained  by  integrating  the 
integrand  of  the  left  side  of  Eq.  (D99)  by  parts  twice.  The  function 
SUn(S2)  (D99)  is  defined  below  in  Eq.  (DlOO). 

IH. 


"u  du  ■ 


dK„(C2U,) 


(DlOO) 


H. 


One  then  calculates  the  indicated  derivatives,  substitutes 
the  results  into  Eqs.  (D98)  and  (D99),  and  uses  Eq.  (D97)  for 
simplification  purposes.  The  two  simplified  equations  are  then  added 
to  yield  the  rather  complex  equation  shown  below. 


‘ / •'n<52l5l>  ^ “r<«2>  “^2  ‘ ^^(52) 


g)  I (C^  - C2) 


where 


^ ° ak  * i * (g^->  )5^  - 2g(g-l  )5g5  - (0102) 


To  simplify  the  forthcoming  eigenvalue  problem  for  the  expan- 
sion functions,  it  is  desirable  to  eliminate  the  integral  involving 
'^n^^2^1^'  this,  one  applies  the  same  sequence  of  operations 

(the  sequence  described  between  Eqs.  (DlOO)  and  (DlOl))  to  the  pair 
of  equations  below, 
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n 


(D103) 


g 

dCi 


where 


The  result  is 

' 

al;’*  '55, 


u„  (5,)  » Y 


"i 


%<«2) 


du„(52> 


(62I5,)  -35 


^1 


u„(52> 


dKn(52l5,) 


d5c 


'*'2  V“n'^2> 


T“n(52>  = |“n‘5l>  ^‘«2l«1> 


Eq.  (D106). 


= ''n 


K„(52l5,) 


^ * '552  - ^ 


u„(52>  J52 


Un(?2)  ^^1  ‘ ^2^  ^^2 


" Vn 


(D104) 


(DlOb) 


(D106) 
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One  then  multiplies  Eq.  (D106)  by  2i5g(l-g)  and  adds  the  result 
to  Eq.  (DlOl)  to  obtain 


b 


21V„(1  -g)  J"  K„(52U,)  2 4 * ^ u„(52)  dS,  (D107) 


where 


iC^d  - g)  n2 


(D108) 


RUn(C2)  = - 2i5g(l  - g)  Tu^($2) 

Then,  using  the  operator  M in  the  manner  discussed  in  chapter  V 
(pages  57-62),  the  eigenvalue  problem  for  the  expansion  functions 
corresponds  to  the  following  differential  equation. 

-JJJ— Mt*  2’5g(i  - g)  -3j- 


(D109) 


+ (g2  . 1)^2  + e . n ^ y (5)  = 0 


(DUO) 


As  a first  step  in  solving  this  equation,  one  lets 

-iCgd  - g)c 

u(5)  = w(5)  e and  substitutes  this  form  into  Eq.  (DUO). 

This  substitution  leads  to  the  following  equation  for  w(c). 
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|(g2-i)c2  - fr+  ♦ (g-D^sd  «U)  = o (oin) 

Then  taking  > 1 , and  letting  z = = 1/2  \g^  - l',  and 

, (s2+(g-l)252| 

T = - >,  one  obtains  Eq.  (D112). 

j ^ ° (D112) 

This  equation,  which  is  identical  to  Eq.  (D65),  has  a general  solution 
of  the  form 


w 


n.T 


„(iz^)  + BW 
H 
2 


(D113) 


For  these  toroidal  resonators,  the  values  of  x are  chosen  so 
that  the  solutions  reduce  to  the  spherical  mirror  expansion  functions 
as  Cg  -*■  0.  This  choice,  which  automatically  incorporates  the  proper 
reduction  to  the  plane  parallel  solutions,  simply  amounts  to  the 
requirement  that  (see  Eq.  (075)) 


- ( ""nm  \ ]_ 

^ + (1  + i)e  ) ^2 

Thus,  for  unstable  resonators,  the  expansion  set  consists  of 
functions  u^(5)  given  by 


^n,.  = 


+i5„(g  - 1)5  1 

A'e  ® n (’“  ^ ) 

ix.  2 

+U«(g  - 1)5  1 

B'e  ® 7^  n 

^ ir  — 

IT,  2 


(0114) 


(0115) 
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where  t is  given  by  Eq.  (D114). 

For  < 1 , one  returns  to  Eq.  (Dill)  and  makes  the  sub- 
stitutions w(c)  = v(5)/5  and  z = to  obtain  Eq.  (0116), 


d^v(z)  1 
dz2  ■ T 


1 + inL^  + 2y_M 

Z2  Z 


V(2) 


(D116) 


where  a'  = >/l  - g^'and  s^  + (g  - 1 )^  5^  = - a' (2v  + 1 ) . As  pre- 
viously indicated  (see  Eq.  (D58)),  the  general  solution  of  this 
equation  is 


V ^ (z) 

ic,p 


AW 


K,P 


(2)  + 


BW  (-z) 
-K.P 


(D117) 


where  < = - v/2  - 1/4  and  p = n/2.  By  requiring  these  solutions  to 
reduce  to  those  for  the  spherical  case  as  5^  -*■  0,  one  obtains  the 
following  expression 


w_ 


n,m 


(5)  = e 


+ne(3 


1)5  - 
e 


a'C' 


.n  + 1 rO 


m 


(D118) 


for  the  expansion  functions  when  g^  < 1 . 
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Analytical  Results  for  Simple  Systems 

The  purpose  of  this  appendix  is  to  present  the  results  obtained 
from  the  analysis  of  three  relatively  simple  resonant  systems.  These 
analytical  results,  which  supplement  the  numerical  results  presented  in 
chapter  VI,  were  obtained  to  support  the  basic  theory  and  method  of 
solution  discussed  in  the  body  of  this  report.  The  three  systems  con- 
sidered are  all  symmetric  resonators  with  perfectly  conducting  mirrors. 
They  are  an  infinite  radial  waveguide,  a stable  resonator  with  rec- 
tangular mirrors,  and  a stable  resonator  with  circular  mirrors. 

For  each  of  these  resonators,  the  analysis  is  performed  in  the 
following  manner.  First,  the  applicable  integral  equation  is  con- 
verted to  an  integral  equation  involving  a surface  integral  over  a 
single  resonator  mirror.  Then  using  the  expansion  functions  presented 
in  chapter  V,  and  replacing  the  finite  limits  of  integration  with 
infinite  limits,  the  integral  is  evaluated  to  yield  an  oscillation 
condition  for  the  resonator  being  considered.  This  oscillation  condi- 
tion is  then  compared  to  similar  conditions  obtained  by  other  authors. 

The  Infinite  Radial  Waveguide 

The  system  considered  in  this  section  consists  of  two  plane 
mirrors  that  are  infinite  in  extent.  The  region  between  the  two 
mirrors,  which  are  separated  by  a distance  L,  contains  a homogeneous 
medium  characterized  by  parameters  u.  e,  and  o,  where  a > 0.  Since 
the  mirrors  in  this  system  are  infinite  in  extent,  finite  fields 
will  be  present  only  if  any  gain  resulting  from  the  polarization  of 
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the  medium  is  less  than  the  loss  resulting  from  the  conductivity  a. 
Thus,  the  propagation  constant  k is  written  in  the  form  k = k'  + ik", 
with  k"  < 0. 


I 

I 

i 

.1 

> 1 

‘ < , 

J 


To  take  advantage  of  the  symmetry,  a sysjtem  of  coordinates  is 
established  so  that  the  origin  is  located  midway  between  the  two  mirrors. 
With  this  system,  a point  on  the  i^^  mirror  has  position  vector 

± L/2  a^,  where  = p^-Sp.  Sp  is  a unit  vector  in  the  radial 
direction,  and  a^  is  a unit  vector  parallel  to  the  optic  axis.  Then, 
using  this  system  of  coordinates  and  Eqs.  (53)  and  (54)  as  the  basic 
forms,  one  substitutes  Eq.  (54)  into  Eq.  (53)  and  applies  the  fact 
that  the  self-induction  integrals  vanish  identically  (cos(n,  r-r')  = 0). 
The  result  is  the  following  expression  for  the  current,  u(F^),  on 
mirror  #1 , 


u(F^) 


'ikf  f - f ^ 

^1  j “("l^  j "°^“12  c°Sa2-,dS2dsj 


(El) 


where  R.j2  = IFpFj'l , R21  = , and  is  the  angle  between  R^^ 

and  the  optic  axis  of  the  guide. 

The  next  step  in  the  analysis  is  to  assume  that  the  current  at 
a point  with  coordinate  on  mirror  #1  is  primarily  a result  of  the 
current  over  a small  region  surrounding  the  point  with  coordinate 
F-]  on  mirror  #2.  This  effective  region  is  assumed  to  be  sufficiently 
small  so  that  the  paraxial  approximation  and  the  condition  shown  in 
Eq.  (E2)  are  satisfied  over  the  entire  region. 
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! Illipipj  1,^ 


-ik 


e 


ik  I — 
^-ikL  ■ 2L  '^1 
e e 


^2' 


(E2) 


As  will  be  seen  later,  these  assumptions,  which  certainly  apply  for 
lossy  systems,  correspond  to  paraxial  modes  within  the  radial  guide. 
Applying  the  paraxial  approximation,  the  distances  R..  are 

' J 

approximated  by  in  the  amplitude  terms  to  yield 

* J 


Substitution  of  Eq.  (E2)  into  Eq.  -(E4)  yields  the  following  expression 
for  the  kernel  K(F^|Fp. 


K(r^|rj)  = e‘ 


ik  ^ 
2ikL  / ■ 2L 


“/ 


jlPl  - P2l‘ 


+ Ip]  - P2I 


’ \ 

\ dS,  (E5) 


To  evaluate  this  integral,  one  makes  the  substitution 


q = V - p. 


(E6) 


where  v = 1/2  (p-j  + p])  and  p = (p^  - p] ) to  obtain  Eq.  (E7) 

^ IF  + T Pp  + IF  - 7 F!^!  _ 


K(v|F)  = e / e 

q 


2ikL  'I"'- I*- 


(E7) 


160 


Writing  out  the  terms  in  the  exponent  and  integrating  over  angles 
yields 

. li  d2  7 . jk  2 

K{v|p)  = e / e ^ q dq  (E8) 

0 


Recalling  that  k = k'  + ik"  with  k"  < 0,  the  integral  in  Eq.  (E8)  is 
easily  evaluated  by  making  a change  of  variables.  Substituting  the 
result  into  Eq.  (E3),  one  obtains  an  integral  equation  in  terms  of  an 
integral  over  a single  mirror. 


(E9) 


To  reduce  this  equation  to  one  involving  only  the  radial 
coordinates,  one  assumes  that  the  modes  have  the  form 


u(5",)  = u„(pj) 


(ElO) 


and  performs  the  azimuthal  integration.  The  result  is  shown  below. 


U„(p,)  = 


^-2ikL  k .n+1 
e 2T  ’ = 


ik  2 
^^1 
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Finally,  substituting  the  form  for  the  plane  parallel  resonator 
expansion  functions  given  in  Eq.  (129)  (J^Ck^p),  where  k^,  the  radial 
component  of  the  wave  vector,  is  chosen  to  be  real  so  the  expansion 
functions  remain  finite  as  p ->  «>)  into  Eq.  (Ell),  and  evaluating  the 
result  (Erdelyi,  ref.  41,  p.  50),  one  obtains 

ikp^  ^ 

Jn(kpP)  |l  - e"”  |=  0 (E12) 

This  condition,  which  must  hold  for  all  p on  the  mirrors,  will 
be  satisfied  if 


k 


P 

~Z1T 


(E13) 


where  q is  an  integer. 

The  exact  results  for  this  problem  (specialized  to  the  lossless 
case)  are  reported  by  Harrington  (ref.  32,  p.  209).  For  real  k^,.and 
with  p = 0 included,  these  results  correspond  to  waveguide  current  modes 
of  the  form, 

Un(p.e)  = J^(kpP)  e-"'"®  (E14) 

where  k'  is  specified  by 


k' 


(E15) 


At  optical  or  infrared  wavelengths,  k'  is  a large  number. 
Thus,  if  kp  is  not  large,  q must  be  a large  integer.  Under  these 
conditions,  the  square  root  is  accurately  approximated  by  using  the 
first  two  terms  of  its  binomial  expansion.  That  is. 
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Lk^ 

k'  ~ I3.+  —p. 
~ L 2Trq 


which  is  very  nearly  equal  to 


(E16) 


k^ 

(E17) 

Then,  specializing  Eq.  (E13)  to  the  case  where  k"  -►  0 and 
comparing  the  result  to  Eq.  (E17),  one  can  see  that  for  paraxial  modes 
(modes  for  which  kp/k'<<l),  the  procedure  for  approximately  solving  the 
integral  equation  yields  essentially  the  same  results  obtained  by  other 
authors. 


Stable  Resonator  Calculations 

The  systems  considered  in  this  section  are  symmetric,  stable 
resonators  (o<g<l)  satisfying  both  the  paraxial  approximation  and  the 
condition  shown  in  Eq.  (E2).  For  these  resonators,  which  are  formed 
by  two  perfectly  conducting  mirrors  of  either  rectangular  or  circular 
cross  section,  the  spatial  dependence  of  the  modes  separates  to  yield 
two  independent  governing  equations.  Each  of  these  equations  is  written 
in  operator  notation  below. 


Ui  = yK^2'^21^1 


(E18) 


where  u-|  represents  the  current  distribution  on  mirror  #1,  and  the 
K. . involve  the  coordinates  of  both  mirrors. 

' J 

As  the  resonators  being  considered  are  symmetric,  the  integral 
operators  K^2  1^21  identical  in  form.  Therefore,  the  eigen- 

vectors of  the  operator  K-, 


Jl  = ± YtK2-|U^ 


(E19) 
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are  also  eigenvectors  of  the  operator  K^2*^21’  ~ ^ 

result,  the  following  analysis  of  these  stable  resonators  is  based 
on  relatively  simple  equations  of  the  general  form  shown  in  Eq.  (E19). 

Rectangular  Mirror  Resonators.  For  resonators  with  rectangu- 
lar mirrors,  the  set  of  equations  corresponding  to  Eq.  (E19)  is  shown 
below  (Eqs.  (Dll)  through  (D13)), 


'i 

2tt 


- + d)  + iC,5o 

u„(Ci)  e 2 ' 2 ' 2 d5, 


-H. 


(E20) 


^m^'^2^  = ^m  V 


, Py  - + n|)  + in-.no 

• / e ^ ^ ^ ^ dn 


H 


(E21) 


- o+itTrq  - kL} 

Y Y 6 

'n’m 


(E22) 


where  c = ^ k/L  x,  n = ^k/T  y»  = yj  k/ll  a,  = yfkJZ  b,  and  q is 
integer. 

To  obtain  the  oscillation’  condition,  the  general  form  of  the 
expansion  functions  for  this  resonator  (see  appendix  D)  is  substi- 
tuted into  Eq.  (E20).  As  that  form  is  negligible  for  large  [x],  the 
finite  limits  of  integration  are  replaced  by  infinite  limits  to 
obtain  the  equation  shown  below. 


an 


= Y 


r— » 

n j °n^^' 


Cl)  e 


■ ^^1  ^2^  ^ ^^1^2 


(Eii3) 
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is  substituted  into  Eq.  (E23)  and  the  order  of  integration  is  inter- 
changed. Then,  the  integral  over  5-j  is  evaluated  by  completing  the 
square  in  the  exponent,  and  the  result  is  simplified  to  yield  the 
following  equation. 


An  identical  procedure  involving  the  y-variation  yields  a 
similar  result  for  By  substituting  these  expressions  into 
Eq.  (E22),  the  desired  oscillation  condition  is  obtained. 
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g+i{irq  - kL}  _ ^-(n  + m + 1) 


(E27) 


Finally,  assuming  that  k is  real,  this  equation  is  manipulated  to 
yield. 


f I cos-'  gl 

mnq  2L  | ^ tt  ^ j 


(E23) 


which  is  identical  to  the  result  reported  by  Pressley  (ref.  43,  p.  433). 

Circular  Mirror  Resonators.  The  appropriate  set  of  equations 
for  these  resonators  is  (see  Eq.  (D44)) 


(E29) 


where 


_ ii  (52  + ^2) 

K„(e,|S2)-e  2 J„(5,52)5, 


(E30) 


and 


nir 


+i{Trq  + -s — kL} 


Yn  = 1 e 


(E31) 


As  before,  the  general  form  of  the  expansion  functions  (see 
Eq.  (064)) 


- — 

%(«2>  = 52"  ^ 


(E32) 


with  a'  = yj]  - g2,  is  substituted  into  the  integral  equation.  Then 
using  infinite  limits  of  integration,  the  integral  is  converted  to  a 
known  form  (Erdelyi,  ref.  41,  vol . 2,  p.  43),  and  the  result  is 
simplified  to  yield  the  following  expression. 
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APPENDIX  F 

Additional  Numerical  Results 


The  purpose  of  this  appendix  is  to  present  the  circular 
mirror  results  not  included  in  chapter  VI  as  well  as  all  of  the 
rectangular  mirror  results.  The  circular  mirror  results  are  pre- 
sented first,  followed  by  those  for  rectangular  mirror  resonators. 

Circular  Mirror  Resonator  Results 

The  only  results  for  circular  mirror  resonators  not  covered 
in  chapter  VI  are  six  mode  plots  in  the  series  for  a plane  parallel 
resonator  with  N = 10.  These  six  plots  include  tne  third  and  fourth 
azimuthal ly  symmetric  (n  = 0)  modes  and  the  first  four  n = 1 modes. 

As  was  the  case  for  the  two  modes  covered  in  chapter  VI,  there  is 
excellent  agreement  between  the  distributions  obtained  by  this  author 
and  those  obtained  by  Fox  and  Li  (ref.  25,  p.  464,  figs.  5 and  6). 

This  agreement  is  especially  good  for  the  relative  phase  distributions 
(in  degrees),  where  the  only  disagreements  of  any  significance  occur 
in  regions  where  the  field  magnitude  is  quite  small. 

These  mode  distributions  (denoted  by  +)  as  well  as  the 
comparative  distributions  from  Fox  and  Li  (denoted  by  *)  are  given 
in  figures  16  through  21  on  the  following  pages. 
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Figure  20.  Relative  Magnitude  and  Phase  Distributions 
for  the  Third  Lowest  Loss,  n = 1 Mode 
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As  with  circular  mirror  resonators,  there  are  two  types  of  data 
to  be  presented.  The  eigenvalue  data,  which  correspond  to  a single 
pass  through  the  resonator,  are  covered  first.  These  data  are  followed 
by  a series  of  mode  plots  for  rectangular  mirror  resonators  with  g = l.d 
Eigenvalue  Data.  The  first  cases  for  which  eigenvalue  data 
were  obtained  were  two  plane  parallel  resonators  with  N = 0 and 
N = 8/ir.  The  percentage  power  loss  for  the  first  several  even  symmetric 
modes  of  these  two  resonators  is  presented  below  in  tables  X and  XI. 

The  comparative  data  were  taken  from  Sanderson  and  Streifer  (ref.  9, 
p.  135,  tables  IV  and  V) . 

Table  X 

Percentage  Power  Loss  for  Rectangular  Mirror  Resonator 
with  g = 1.0,  N = 10  for  Even  Symmetric  Modes 


Table  XI 

Percentage  Power  Loss  for  Rectangular  Mirror  Resonator 
with  g = 1.0,  N = 8/tt  for  Even  Symmetric  Modes 


'.J'. 


As  reported  in  chapter  VI,  the  difference  for  any  of  the  first  three 
modes  is  less  than  2.5%. 

The  next  cases  considered  were  symmetric  rectangular  mirror 
resonators  with  g = 1.2.  The  percentage  power  loss  for* the  first  three 
odd  symmetric  modes  for  three  values  of  is  displayed  in  table  XII. 
For  this  series,  the  comparative  data  were  taken  from  Sanderson  and 
Streifer  (ref.  42,  p.  2131,  fig.  7). 

Table  XII 


Percentage  Power  Loss  for  Rectangular  Mirror  Resonator 
with  g = 1.2  for  Odd  Symmetric  Modes 


DOUGHTY 

SANDERSON  & STREIFER 

f'e 

0.542 

0.54 

1.0 

0.681 

0.68 

0.913 

0.92 

0.562 

0.57 

l.t) 

0.677 

0.68 

0.793 

0.78 

0.60 

2.0 

0.668 

0.68 

0.773 

0.72 

Finally,  eigenvalue  data  for  a wida  range  of  equivalent  Fresnel 
numbers  for  symmetric  resonators  with  g * 1.8  are  presented  in  tables 
XIII,  XIV,  and  XV.  The  data  listed  correspond  to  the  percentage  power 
loss  for  several  even  symmetric  modes  (Sanderson  and  Streifer,  ref.  42, 
p.  2132,  fig.  10). 
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Table  XIII 


Percentage  Power  Loss  With  g = 1.8  for  Even  Symmetric 
Modes  at  Three  Non-integral  Values  of  2Ng 


DOUGHTY 

SANDERSON  & STREIFER 

2Ne 

1 

0.639 

0.64 

4.5 

J 

0.864 

0.90 

r 

0.923 

0.92 

I 

0.997 

1.0 

1 0.717 

0.75 

7.5 

i 0.909 

0.86 

0.924 

0.94 

0.947 

0.95 

0.794 

0.72 

28/ IT 

f 

0.794 

0.77 

0.927 

0.96 

0.979 

0.98 

1 

Table  XIV 

■ 

Percentage  Power  Loss 

with  g = 1.8  for  Even 

Symmetric 

r 

Modes  at 

Three 

Integral  Values  of  2N^ 

DOUGHTY 

SANDERSON  & STREIFER 

2Ne 

0.712 

0.72 

6.0 

0.866 

0.31 

[ 

y 

1 

0.933 

0.98 

0.968 

0.98 

\ 

0.647 

0.66 

9.0 

E 

0.916 

0.90 

1 

0.954 

0.91 

i 

0.959 

0.96 

i 

0.738 

0.74 

12.0 

0.858 

0.86 

1 

0.869 

0.91 

0.937 

0.94 

1 

■ 

1 

[ 
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Table  XV 

Percentage  Power  Loss  with  g = 1.8  for  Even  Symmetric 
Modes  at  Large  Values  of  2Ng 


DOUGHTY 

SANDERSON  & STREIFER 

2ttN^ 

e 

0.652 

0.65 

34.0 

0.899 

0.91 

0.931 

0.93 

0.709 

0.68 

36.0 

0.899 

0.91 

0.924 

0.92 

0.651 

0.65 

40.0 

0.861 

0.89 

0.905 

0.91 

Examination  of  the  results  in  these  three  tables  reveals  that  the 
differences  are  typically  4 percent  or  less.  However,  the  maximum 
difference  was  7 percent,  which  occurred  for  the  second  mode  with 

2Ng  = 6.0. 

Mode  Distributions.  To  further  check  the  analysis  presented 
in  this  paper,  a series  of  mode  plots  was  made  for  resonators  with 
g = 1.8.  These  mode  plots  (figs.  22  through  28)  include  the  relative 
intensity  and  phase  ( radians  ) for  the  first  two  even  symmetric  modes 
for  2nNg  = 34,  36,  and  40  and  the  relative  intensity  distributions  for 
2TrNg  = 18.  For  this  last  case,  relative  phase  data  were  not  included 
as  they  were  not  available  in  the  paper  from  which  the  comparative 
data  were  taken  (Sanderson  and  Streifer,  ref.  42,  p.  2133,  figs.  13 
through  20) . 

After  comparing  the  results  for  these  cases,  several  general 
remarks  are  offered  in  the  next  few  pages.  First,  the  intensity 
distributions  obtained  in  this  analysis  are  rather  slowly  varying. 
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Second,  the  relative  values  of  some  of  the  peaks  and  troughs  are 
somewhat  different.  Although  this  is  especially  true  at  points 
where  the  intensity  is  low,  it  is  certainly  not  uncomnon  at  other 
points.  These  two  facts  are  likely  the  result  of 

1.  The  tendency  of  a variational  method  to  average  the 
field  distributions  when  the  convergence  criterion  is 
based  only  on  the  eigenvalues,  and 

2.  Computational  errors  in  determining  the  expansion 
functions.  This  last  fact  also  causes  errors  in  the 
expansion  coefficients,  which  might  further  explain  any 
errors  or  differences. 

However,  despite  these  differences,  sharp  disagreement  does  not 
exist  for  any  of  these  intensity  distributions.  In  fact,  the  basic 
behavior  or  nature  of  the  intensity  distribution  is  correctly  pre- 
dicted in  every  case. 

Further  examination  of  these  plots  reveals  that,  except  for 
the  second  mode  with  2irNg  = 34,  there  is  basic  agreement  between  the 
phase  distributions  for  all  modes.  For  the  one  case,  the  disagreement 
occurs  over  a region  in  which  the  phase  is  changing  quite  rapidly  so 
that  the  difference  could  result  from 

1.  Inability  of  the  variational  method  to  follow  such 
rapid  changes  with  the  number  of  functions  used, 

2.  Computational  errors,  or 

3.  2it  ambiguities  in  the  phase  data. 

Interestingly,  for  regions  where  the  phase  is  not  changing  so  rapidly, 
elimination  of  these  2ir  ambiguities  results  in  excellent  agreement. 
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In  any  case,  improved  computer  routines  would  allow  the  use  of  more 
expansion  functions  and  would  likely  result  in  excellent  agreement  for 
these  modes. 

Since  these  and  other  results  in  this  paper  are  discussed  in 
the  last  section  of  chapter  VI,  no  discussion  is  included  with  the 
mode  plots  presented  in  this  appendix.  In  these  plots,  the  results 
obtained  by  Sanderson  and  Streifer  are  denoted  by  an  * while  those 
obtained  by  this  author  are  denoted  by  a +. 
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Figure  22.  Relative  Distributions  for  the  First 
Two  Modes  with  g = 1.8  and  2TiNg  * 18 
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Figure  25.  Relative  Distributions  for  the  Lowest 

Loss  Mode  with  g = 1 .8  and  2irN  = 36 
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Figure  26.  Relative  Distributions  for  the  Second 
Lowest  Loss  Mode  with  g = 1.8  and  2itN 
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Figure  27.  Relative  Distributions  for  the  Lowest 
Loss  Mode  with  g = 1.8  and  2ttN  = 40 


Figure  23.  Relative  Distributions  for  the  Second 
Lowest  Loss  Mode  with  g = 1.8  and  2Trt| 
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APPENDIX  G 

Mathematical  Approach  for  Solving 
the  Matrix  Problem 

In  chapter  IV,  the  problem  of  determining  the  modes  of 
paraxial  resonators  was  reduced  to  solving  a particular  matrix  equa- 
tion. The  equation,  which  involves  the  symmetric  matrices  K and  'a 
and  the  column  matrix  A,  is  shown  below. 

KA  = y'HA  (Gl) 

The  purpose  of  this  appendix  is  to  describe  the  procedure  for  solving 
Eq.  (Gl)  such  that  the  resulting  approximations  to  the  modes,  u(?), 

J 

obey  the  orthogonality  condition 

(Uj(c)|Uk(c))  = 

where 

(Uj(c)lUk(c))  = y Uj(c)  Uk(t)  dc  (G3) 

c 

However,  before  describing  the  procedure,  the  method 'used  to  obtain 
Eq.  (Gl)  is  briefly  reviewed. 

The  original  problem  of  solving  for  the  resonator  modes, 
u(c),  was  formulated  in  terms  of  the  integral  equation 

Ku(c)  = 7u(c)  (G4) 

where  K denotes  the  appropriate  integral  operator,  and  7=  . To 

solve  this  equation,  the  modes  were  expanded  in  terms  of  a set  of 
known  functions  {>((^(1:)}  with  coefficients  A^, 
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r' 
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1 

w • ; 

H 


u(i)  > E 


(65) 


and  the  Rayleigh-Ritz  procedure  was  applied.  That  procedure  led  to 
Eq.  (Gl)  with  the  elements  of  the  matrices  K and  'V  given  by 


(G6) 

(G7) 


The  first  step  in  solving  this  matrix  problem  is  to  reduce 
Eq.  (Gl)  to  an  eigenvalue  problem  involving  a symmetric  matrix  B.  To 
do  this,  one  decomposes  the  symmetric  matrix  'V  using  a Cholesky 
decomposition  (Wilkinson,  ref.  44,  p.  229).  This  decomposition  yields 
the  form, 

T 


>?  = L L' 


(68) 


T 


where  Lisa  lower  left  triangular  matrix  and  L is  the  transpose  of 
L. 

To  determine  the  matrix  L,  it  is  helpful  to  write  Eq.  (G8) 
in  the  expanded  form  shown  below. 


■^11 

^12- 

1 

C 

1 

1 

1 

O 

0 

1 

^11 

*■12  ■■  ^in 

^12 

1 

1 

1 

CM 

CM  - 

^2n 

1 

= 

^12 

1 

^22 

1 

1 

0 

0 

1 

*■22 

t 

/in 

1 

^2n- 

1 

f 

nn  ^ 

1 

. ^In 

1 

^2n” 

^nn  _ 

1 

0 

C 

C 

1 

1 

o 

(69) 
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Then,  perfonning  the  multiplications  involving  the  first  row  of  L, 
one  obtains  the  relation 


'‘'li  “ ^11  ^li 


(GIO) 


for  i = 1,2,— n.  These  equations  can  be  solved  for  the  elements  of 
the  first  row  of  L^.  To  compute  the  (n  - 1)  unknown  elements  of 
row  2,  one  uses  the  (n  - 1)  known  elements,  £-|^,  and  the  relation 


'‘'k2  ^ *"12  ^•ik 


This  process,  which  involves  the  general  relation 


(Gil) 


'*'kr  ~ ^5] 


(G12) 


is  continued  until  all  of  the  have  been  determined. 

Using  that  decomposition,  one  premultiplies  Eq.  (Gl)  by  L' 
and  makes  the  substitution 


A = ( J)’^F 


(G13) 


With  these  manipulations,  one  obtains  the  eigenvalue  problem. 


BF  * yF  (G14) 

where  B is  a symmetric  matrix  given  by  B = L~^  K 

To  solve  this  eigenvalue  problem,  B is  reduced  to  a tri- 
diagonal matrix  using  Householder's  method  (Acton,  ref.  4b,  p.  324). 
That  method,  which  preserves  symmetry  and  eigenvalues,  involves  the 
use  of  symmetric,  orthogonal  transformations  as  described  on  tne 
following  page. 


r\ 
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In  the  first  step  of  the  procedure  to  tridiagonal ize  the  matrix 

b 
b 

(615) 


11 

bi2  - - - 

^in 

12 

1 

^22 

1 

1 

1 

In 

1 

“^2n  --- 

*^nn 

one  finds  a symmetric,  orthogonal  matrix  P-j , such  that 

b. 


1 


'12 


Hn 


1 


(616) 


where  f?  = ^ b?..  The  matrix  P,  is  then  applied  in  the  manner  shown 
' i=2  1 

below. 
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“ 
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“ 

1 

0 — 

0 

1 

0 

— 

0 

‘^ll  ^ 

0 — 

0 

0 

0 

^1 

1 

1 

Pi 

H 

( 

1 

Pi 

= 

0 

1 

c 

0 

- 

_0 

- 

- 

.0 

- 

(617) 


where  C is  the  (n  - 1)^^  order  matrix 

C 


C = 


'11 


'12 

I 


^12  ““"^Un-l 


C22 


L^l,n-1  ^2,n-l-^n-l  ,n-l 


(618) 
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This  procedure  is  then  repeated  to  yield  a matrix 


(G19) 


n-1 

where  f ? = ^ C? . . 

i=2 

result  shown  below. 


The  matrix  is  then  applied  to  yield  the 


['1  w 

W b] 


“11  -- 


[']  W 

H U 


0 — 0 

fl  — 0 

0 f ) r T 


0 0 


(G20) 


where  I is  the  identity  matrix.  This  procedure  of  applying  symmetric, 
orthogonal  transformations  P^. , is  continued  until  a tridiagonal  matrix 
T is  obtained.  The  method  used  to  obtain  the  P^.  is  covered  in  Acton 
(ref.  45,  pp.  326-329). 

Once  the  matrix  T is  known,  the  eigenvalues  are  determined  using 
the  "LR"  algorithm  (Acton,  ref.  45,  p.  350).  This  algorithm,  which  is 
not  normally  used  for  a general  matrix,  was  chosen  because  the  tri- 
diagonal  form  is  preserved  at  every  step  of  the  iteration  procedure. 

This  preservation  of  form  makes  it  possible  to  accomplish  each  iteration 
using  only  (n  - 1)  multiplications  and  (n  - 1)  divisions. 
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To  apply  this  algorithm,  one  first  decomposes  the  matrix  T 


such  that 


T = LR  ( 

where  L and  R have  the  following  characteristics: 

1.  L is  a lower  left  triangular  matrix, 

2.  R is  an  upper  right  triangular  matrix,  and 

3.  The  elements  on  the  main  diagonal  of  L have  the 
value  one. 

Then,  the  matrix  T-j  is  computed  according  to 

T^  = RL  ( 

and  the  procedure  is  repeated.  After  the  (s  + 1)^*^  iteration,  one 


obtains 


such  that 


Ts  -M  = «s'-s 


+ 1 = ‘■s'*  (g: 

This  process  is  continued  until  the  set  {T^}  converges  to  an  upper 
right  triangular  matrix  with  the  eigenvalues  on  the  main  diagonal. 

In  the  next  step,  these  eigenvalues  are  used  to  compute  the 
eigenvectors  of  T by  employing  a procedure  which  utilizes  all  n 
equations  of  the  homogeneous  system 

(T  - Y|^I)  = 0 (G^ 

where  Wj^  is  a column  matrix  with  n elements.  This  inverse  iteration 
procedure  (Acton,  ref.  45,  p.  357)  is  based  on  the  iteration  scheme, 

(T  - v,I)  vj  , , - (Gi 
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where  (the  approximation  to  the  eigenvector  of  T)  is 
determined  using  Gaussian  elimination  with  interchanges.  As  the  eigen- 
values Y|(  have  already  been  determined,  this  iteration  procedure 
rapidly  converges  to  the  desired  eigenvectors.  These  solutions  are 
then  used  to  obtain  the  eigenvectors  (of  B)  by  applying  the  same 
similarity  transformation  used  to  obtain  T from  B.  The  Aj^,  which  are 
obtained  by  applying  the  relation  shown  below. 


A.  = (J)'^  F, 


(G27) 


are  then  substituted  into  Eq.  (G5)  to  obtain  the  resonator  modes. 
To  show  that  the  eigenfunctions  u.(c)  are  orthogonal  (see 

w 

Eq.  (G2)),  one  computes  the  scalar  product 


(Ui(5)lut(c))  = E A.  (ij^„(c)|ij;^(c))  A 


jq  "km 


(G28) 


which  can  be  written  in  matrix  form  as  shown  below. 

(Uj(c)|U|^(c))  = aJ'fAi^ 

One  then  applies  Eq.  (G27)  to  obtain  the  following  result 
(Uj(c)iu,^(c))  = FjJ  L‘^ 

Then  using  Eq.  (G8),  one  obtains 


(Uj(;)|U|,(t))  - Fj\ 


(G29) 


(G30) 


(G31) 


Since  the  eigenvectors  F.  of  a symmetric  matrix  B obey  the  relation, 

J 


fk  ' «jk 


the  modes  u.(c)  obey  Eq.  (G2)  as  asserted. 
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APPENDIX  H 


The  Confluent  Hypergeometric  Functions 

As  was  shown  in  chapter  V,  the  expansion  functions  for  the 
modes  of  both  rectangular  and  circular  mirror  resonators  are  directly 
related  to  the  confluent  hypergeometric  functions  of  the  first  and 
second  kind.  The  purpose  of  this  appendix  is  to  present  some  of  the 
more  important  characteristics  and  relationships  involving  those 
functions.  For  more  detailed  discussions,  the  reader  is  referred 
to  the  sources  of  this  appendix  (Lebedev,  ref.  40,  chaps.  9 and  10; 
Erdelyi,  ref.  41,  chaps.  6,  8,  and  10). 

As  the  name  implies,  these  functions  are  solutions  of  the 
confluent  hypergeometric  equation  which  is  shown  below  (Lebedev, 
ref.  40,  p.  262). 

^^-«u(z)=0  (HI) 

The  solution  of  the  first  kind,  denoted  $(a,  y,  z) , is  an  entire 
function  of  z,  but  it  does  have  simple  poles  at  the  points 
Y = 0,-1, -2,  — . For  values  of  y ^ 0,-1 ,2, — , this  function  is 
defined  by  the  series  (Lebedev,  ref.  40,  p.  260), 

E(a)|.Z 
k=0  '■^^k 

where  (t)^  = 1,  (t),^  = and  T = a or  y. 

Unfortunately,  the  solution  of  the  second  kind,  denoted 
f(a,  Y,  z),  is  somewhat  more  complex.  For  values  of  y such  that 
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Y t 0,±1,±2, — this  function  is  defined  in  terms  of  the 
$(a,  Y.  z)  as  shown  below  (Lebedev,  ref.  49,  p.  263), 


Y.  z) 


$(a,  Y^  z) 


+ + a - Y.  2 -Y,  Z)  (H3) 


with  [arg  i\  < tt.  However,  for  y ® n + 1,  where  n is  zero  or  a positive 
integer,  this  expression  is  indeterminant  and  f(a,  n + 1 , z)  must  be 
evaluated  using  a limiting  procedure.  The  procedure,  which  involves 
the  substitution  of  the  series  representations  for  $ and  r into 
Eq.  (H3)  as  well  as  the  application  of  L'Hospital's  rule,  leads  to 
the  following  equation  (Lebedev,  ref.  40,  p.  264), 


/ T^n+1  ^ (a)(.  _k 
'i'(a,n+l,z)  = 2^  (n+k)!k! 

k=0 


i|)(a+k)  -ij)(l+k)  -i(;(n+l+k)+ln(z) 


+ 


V (-’)^n-k-l)!  (g-n),  „ 

k=0 


(H4) 


where  i|;(x)  is  the  logarithmic  derivative  of  the  gamma  function 
(Erdelyi,  ref.  41,  p.  44)  and  jarg  z]  < u.  In  addition,  Eq.  (H4) 
can  be  used  in  conjunction  with  a recurrence  relation  (Lebedev,  ref.  40, 
p.  265)  to  define  4'(a,  y.  z)  for  all  values  of  y*  As  a result, 

'i'(a,  Yj  z)  is  not  only  an  analytic  function  of  z in  the  plane  cut 
along  {-  «,  0},  but  it  is  also  an  entire  function  of  a and  y. 

For  values  of  y such  that  y f 0,-1, -2, — and  for  values  of 
z such  that  jarg  zj  < u,  one  can  calculate  the  Wronskian  of  these 


two  solutions  to  obtain  the  result  shown  below  (Lebedev,  ref.  40, 


p.  265). 


W{*(o,Y,z);  f(a,Y,z)}  = - z"^e^ 


From  this  expression,  one  can  see  that  if  a ^ 0,-1, -2, — , then 
$(a,  Y.  z)  and  yCo,  Y.  z)  are  linearly  independent  solutions  of 
Eq.  (HI). 

In  addition  to  the  series  representations  shown  above,  * 
and  f are  also  represented  by  asymptotic  expansions  as  lz|  (or  one  of 
the  parameters)  becomes  large.  For  example,  if  jzj  is  large, 

Y f 0,-1, -2,  — , and  |arg  z[  - 6,  then  (Lebedev,  ref.  40,  p.  271), 


•(o.v.z)  ■ rT^  Z 1-')  k 

h=0  k! 

+ 0 ( izr"-'' A + IW-  y.  _-k 

\ / k=0  k!  ^ 

+ 0 (H6) 

where  the  plus  sign  is  chosen  if  Im  (z)  > 0 and  the  minus  sign  is  chosen 
if  Im  (z)  < 0.  Similarly,  one  can  show  that  for  the  identical  condi- 
tions on  z (Lebedev,  ref.  40,  p.  270), 


n (-1)  (a).  (l+a-Y)|^  I- 
4'(ot,Y,z)  = z"“  2^  z + 0 


(u,--) 


Another  asymptotic  relation  involving  $(a,  y.  z)  can  be 
obtained  by  letting  the  magnitude  of  o become  large.  For  that  case. 
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which  is  applied  several  times  in  this  report,  one  lets  k = ^ - a. 
Then  if  y and  tcz  are  bounded  and  if  jarg  z-arg  (Erdelyi, 

ref.  41 , p.  280) , 

1-Y  z 

®(a,Y,z)-*-  r(y)  (<z)  l(2i/<F”  ) + 0 (IkT^)  ( 

\<\^ 

In  addition  to  these  expressions,  there  are  many  important 
relationships  between  the  confluent  hypergeometric  functions  which 
have  played  a key  role  in  this  analysis.  Several  of  these  relation 
ships  are  listed  in  the  following  paragraph. 

The  first  two  relationships  involve  the  Whittaker  functions 
of  the  first  and  second  kind  (Erdelyi,  ref.  41,  p.  264), 

z 

M (z)  = e"^  «(%-k+u,2u+1 ,z)  ( 

f M 

z 

W (z)  = z^"^^  e"^ '?(%-K:+ii,2u+l  ,z)  ( 

where  jarg  z|  < tt.  Next,  an  important  relation  involving  the 
parabolic  cylinder  function  of  order  v,  D^(z)  is  shown  in  Eq.  (Hll) 
(Erdelyi,  ref.  41,  vol . 2,  p.  117). 


0 (z)  = 2^  e ^ 


r(%)  * / 1 

'n  ltv  * I - ^ I 


! 
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